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ABSTRACT 


r 

I 
I 

I  J 

The  procedure  for  solving  large-scale,  nonlinear  multiobjective 
|  optimization  problems  requires  the  use  of  a  computer  coded  algorithm 

which  is  both  computationally  efficient  and  mathematically  valid. 

I  The  multicriteria  methodology,  Goal  Programming,  has  proven  to  be  an 

j  effective  technique  for  the  solution  of  such  problems. 

The  problem  of  this  study  was  the  synthesis  of  the  goal  program- 
j  ming  methodology  into  an  efficient  computer  algorithm  which  could  be 

used  by  the  problem  solver  to  solve  a  general  class  of  nonlinear, 
j  multiobjective  (NLGP)  engineering  design  problems.  Three  such  codes 

were  developed  into  a  computer  package.  Each  code  is  capable  of  solv- 
^  ing  an  NLGP  problem.  Each  code  has  modifications  which  make  it  more 

|  attractive  than  the  others  for  application  to  particular  types  of  NLGP 

problems.  The  results  show  that  these  codes  are  effective  tools  for 

# 

i  solving  nonlinear  goal  programming  design  problems. 
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CHAPTER  I 
INTRODUCTION 


1. 1  Problem  Statement 

The  design  of  the  majority  of  large-scale,  complex  engineering  systems 
requires  the  specification  of  parameters  which  will  best  achieve  a  system 
of  multiple  and  conflicting  design  objectives.  These  objectives  may  be 
linear  and  nonlinear  in  nature.  Real  systems  can  often  contain  hundreds 
and  often  thousands  of  specific  design  objectives.  The  design  parameters 
which  optimize  the  performance  of  these  objectives  can  number  in  the  thou¬ 
sands.  Unfortunately,  traditional  methods  of  modeling  and  optimizing  such 
systems  are  restricted  to  problems  encompassing  a  single  design  goal.  Com¬ 
putational  techniques  for  these  methods  are  quite  often  further  restricted 
to  either  small  problems  or  uncomplicated  objectives,  or  both. 

A  satisfactory  technique  for  the  solution  of  large-scale  engineering 
design  problems  with  multiple,  conflicting  objectives  is  contained  within 
the  multicriteria  tool  of  goal  programming  [18].  This  computational  tech¬ 
nique  has  been  applied  to  a  wide  variety  of  problems.  The  computer  code 
presented  in  [18]  is  intended  for  very  small  problems,  while  those  in  [2] 
and  [14]  have  been  adapted  to  meet  the  design  requirements  of  very  specific 
problems.  There  is  a  need  for  an  improved  nonlinear  goal  programming  (NLGP) 
computer  code  which  can  be  used  without  modification  to  solve  a  wide  variety 
of  general  large-scale  nonlinear,  multiobjective  design  problems.  The  pur¬ 
pose  of  this  study  is  to  develop  such  a  computer  code. 

1.2  Specific  Objectives 

The  specific  objectives  of  this  study  are: 

1.  To  develop  a  computer  code  (or  codes)  which  can  solve  NLGP 
problems  with  up  to  thousands  of  objectives  and  thousands 
of  variables. 
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2.  To  test  and  validate  the  code(s)  for  computational  efficiency, 

3.  To  supply  a  User's  Manual  describing  in  detail  the  use  of  the 
code(s) . 

1.3  Summary  of  Results 

The  research  work  resulted  in  the  development  of  three  computer  codes. 
All  three  are  similar  in  technique.  Each  code  has  minor  modifications  from 
the  others  which  makes  it  more  suitable  for  application  to  particular  NLGP 
problems. 

The  codes  were  tested  on  a  variety  of  problems  varying  in  the  complexity 
of  the  nonlinear  objectives  and  problem  size  (the  number  of  variables  and 
objectives).  The  results  obtained  indicate  that  each  code  is  computationally 
correct  and  yielding  satisfactory  results. 

This  report  contains  information  about  the  development  of  the  computer 
codes  with  details  about  their  structure  and  performance.  A  User’s  Manual 
has  been  written  which  details  specifically  how  to  use  the  three  computer 
codes. 


1.4  Overview  of  Following  Chapters 

An  analysis  of  the  results  of  this  research  are  contained  in  the 
following  chapters. 

Chapter  II  contains  a  background  sketch  of  the  general  nonlinear  goal 
programming  model.  It  lists  advantages  and  limitations  of  the  technique 
in  solving  design  problems.  The  third  chapter  then  provides  a  description 
of  the  computer  codes,  their  limitations  and  modifications.  Some  examples 
are  presented  in  Chapter  IV,  along  with  the  empirical  results  from  the 
testing  of  the  codes.  Finally,  some  conclusions  about  this  study  and 
recommendations  for  further  research  are  provided  in  Chapter  V. 
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CHAPTER  II 

THE  MATHEMATICAL  PROGRAMMING  PROBLEM 

2.1  The  General  Nonlinear  Goal  Programming  Model 

The  primary  emphasis  of  this  research  has  been  directed  towards  the 
area  of  multicriteria  decision  analysis,  and  in  particular,  nonlinear 
goal  programming.  For  this  reason,  a  discussion  of  the  considerable  num¬ 
ber  of  other  mathematical  programming  optimization  techniques  will  not  be 
presented.  The  entire  focus  of  discussion  will  be  on  multicriteria  tech¬ 
niques. 

The  multicriteria  technique  called  Goal  Programming  was  first  presented 
by  Charnes  and  Cooper  [1]  in  the  1950's.  Ignizio  [18]  later  expanded  the 
method  to  solve  nonlinear  as  well  as  linear  models.  More  recently,  this 
technique  has  found  a  broad  range  of  applications  to  solving  design  problems 
in  many  areas  [3,13,14,15,16,17,19,20,21,22,23,24]. 

The  general  nonlinear  goal  programming  model  can  be  described  by  the 
following  [18]. 

A  set  of  mathematical  equations  are  established  by  the  analyst  or 
design  engineer  which  describe  the  problem  as  functions  of  the  design 
parameters.  These  design  equations  are  the  objectives,  some  or  all  of 
which  may  be  nonlinear  in  nature.  The  design  parameters  are  the  deci¬ 
sion  variables.  Associated  with  the  objectives  are  desired  aspiration 
levels.  They  represent  the  performance  of  the  system  at  the  optimal 


design  parameters. 
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Define  the  following: 

(1)  {f^(x),  f2(x),  ....  f^x)}  is  the  set  of  objectives 

(2)  x  =  (x^,  X2»  . ..,  Xj)  is  the  vector  of  decision  variables 

(3)  {b^,  b 2,  . ..,  b^}  is  the  set  of  aspiration  levels  for 

the  objectives. 

Then, 

> 

f^x)  =  bi,  i  =  1,  2,  .  . . ,  m 

< 

represent  the  objectives  or  goals  for  the  system. 

It  is  desired  to  find  the  values  of  the  decision  variables  such  that 
the  objectives  meet  the  specified  aspiration  levels.  Since  this  is  not 
always  possible  due  to  the  conflicting  nature  of  many  objectives,  devia¬ 
tion  variables  are  added  to  the  design  objectives.  These  variables 
measure  the  difference  between  the  desired  and  the  actual  aspiration 
levels. 

Then, 

ft(x)  +  ni  -  pi  =  bi,  i  -  1,  2,  ....  m 

are  the  design  objectives  where  n^  measures  the  negative  deviation  from 
the  aspired  level  for  objective  i,  and  p^  measures  the  positive  deviation 
for  objective  i. 

A  major  assumption  underlying  many  versions  of  goal  programming  is 
that  the  analyst  can  establish  preemptive  priorities  for  each  objective 
or  group  of  objectives.  This  implies  that  the  objectives  can  be  ordered 
such  that  achievement  of  the  objectives  at  any  one  priority  level  is 
immeasurably  preferred  to  the  achievement  of  the  objectives  of  any  lower 


level. 


Any  objectives  which  are  absolute  (i.e.,  they  must  be  satisfied)  are 
assigned  to  priority  level  (P^).  The  remaining  objectives  are  then  grouped 
and  assigned  priority  levels  (i.e,,  P^,  ...»  P^). 

The  goal  programming  algorithm  first  satisfies,  as  nearly  as  possible, 
the  objectives  with  the  highest  priority  level.  It  then  proceeds  to  satisfy 
the  objectives  of  the  next  priority  level,  as  nearly  as  is  possible,  without 
degrading  the  achievement  of  any  objective  in  a  higher  priority  level.  This 
process  continues  until  all  priority  levels  have  been  considered. 

The  satisfaction  of  these  various  priority  levels  is  measured  by  an 
achievement  vector,  a.  This  vector  has  one  component  for  each  priority 
level.  Each  component,  a^,  is  a  weighted  linear  function  of  the  deviation 
variables  which  are  contained  in  the  objectives  of  that  priority  level  m. 

Let  (g1(n,p),  g2(n,p),  ...,  g^n.p)}  represent  the  set  of  these 
linear  functions,  one  for  each  priority  level.  Each  function  can  be  repre¬ 
sented  as: 

8k<"»P)  =  ±t\  (wini  +  v^ 

k 

the  weighted  sum  of  negative  and  positive  deviation  variables  for  all  the 
objectives  contained  in  priority  level  P^,  where  and  w^^  are  the  weights. 

Since  both  negative  and  positive  deviation  variables  are  required  for 
each  objective,  minimization  of  these  variables  will  attempt  to  achieve 
the  aspiration  levels  of  the  objectives.  The  achievement  vector  a  is 
structured  as  an  ordered  set  such  that  the  preemptive  priority  structure 


is  maintained. 
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the  standard  form  of  the  nonlinear  goal  programming  model  is  as 
follows: 

Find  x  =  (x^,  x^,  ) 

to  minimize  a  =  (g^n.p),  g2(n,p),  ....  gk(n,p)} 
subject  to 

fi(x)  +  ni  -  pt  =  biS  i  =  1,  2,  m 

are  satisfied,  n  and  p  are  both  nonnegative, 

(1)  gk(n,p)  is  a  weighted  linear  function  of  the  deviation 

variables 

(2)  k  is  the  priority  level  associated  with  the  objective 

(3)  the  dimension  of  a  represents  the  number  of  preemptive 
priority  levels,  (k) , 

and 

(A)  the  number  of  preemptive  priorities  are  equal  to  or  less 
than  the  number  of  objectives. 

The  value  of  a  will  be  equal  to  the  zero  vector  if  all  the  objectives 
meet  their  aspiration  levels.  The  value  of  a^  will  be  a  positive  value  if 
one  or  more  objectives  in  priority  level  k  are  not  met.  The  goal  program¬ 
ming  solution  is  considered  optimal  if  it  is  the  lexicographic  minimum 
(i.e.,  the  value  of  a  is  lexicographically  the  same  as  or  less  than  the 
value  of  a  for  any  other  feasible  solution).  If  a^  is  non-zero,  an  abso¬ 
lute  objective  cannot  be  met,  and  the  solution  is  considered  unimplement- 
able. 

To  further  clarify  the  general  nonlinear  goal  programming  model,  an 


example  problem  is  presented  next. 
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2.2  A  Target  Allocation  Example 

The  example  problem  selected  to  demonstrate  the  formulation  of  an  NLGP 
model  is  a  target  allocation  problem  which  was  previously  formulated  and 
solved  by  Rice,  Bracken,  and  Pennington  (28)  using  a  single  objective  optim¬ 
ization  technique,  the  sequential  unconstrained  minimization  technique  [5,6]. 

This  problem  examines  the  allocation  of  weapon  systems  among  targets 
with  emphasis  being  placed  on  maximizing  the  amount  of  damage  inflicted  on 
enemy  targets.  The  weapon  systems  in  this  problem  are  limited  to  two  types 
of  attack  aircraft  originating  from  two  aircraft  carriers,  one  conventional 
and  one  nuclear.  Additionally,  the  number  of  each  aircraft  is  considered 
fixed.  This  creates  a  scenario  in  which  the  primary  objective  of  the  prob¬ 
lem  is  to  allocate  a  fixed  inventory  of  aircraft  among  carriers  so  as  to 
maximize  total  damage  on  the  enemy  targets. 

The  following  parameters  and  variables  are  introduced  for  this  problem: 


x.  . 


a 


jk 


y 


ijk 


the  number  of  aircraft  of  type  j  assigned  to 
carrier  i,  i  =  1,2  and  j  =  1,2 

the  fraction  of  target  k  (either  in  numbers 
of  specific  targets  or  target  area)  left 
undamaged  following  a  mission  by  an  aircraft 
of  type  j,  j  =  1,2  and  k  =  1,2,  ...,  5 

the  number  of  missions  flown  against  target  k 
by  an  aircraft  of  type  j  from  carrier  i, 
i  =  1,2,  j  =  1,2  and  k  =  1,2 . 5 


Additional  parameters  are  defined  for  the  objectives  with  regards  to 
the  fixed  supply  of  aircraft  and  carriers: 

bj  31  total  number  of  available  aircraft  of  type  j 

c,  ■  the  carrier  space  required  by  each  aircraft 

^  of  type  j 

d^  a  total  space  available  for  plane  storage  on 
carrier  i 

f, >k  ■  the  fraction  of  total  missions  available  in  the 
J  planning  period  used  in  flying  one  sortie  against 

target  k  by  an  aircraft  of  type  j  from  carrier  i. 


The  objectives  for  this  problem  are  formulated  as  follows: 


Aircraft  availability  constraints 
2 

I  x-m  I  bi  for  j  “  1.2 

i-1  J  J 


Carrier  space  constraints 
2 

T  c.x. .  <  d .  for  i  =  1,2 
. L .  j  13  —  1 
j=l  J 


Limitations  on  missions  due  to  aircraft  availability 
5 


fijk  yijk  -  0 


for  i  =  1,2 

j  -  1,2 


Damage  inflicted  on  targets 

The  total  damage  inflicted  is  the  sum  of  damages  to  all  targets. 

The  damage  to  any  individual  target  may  be  expressed  in  terms  of  the 

parameters,  a..  .  The  damage  inflicted  on  target  k  by  missions  of  all 
J  * 

aircraft  from  all  carriers  may  be  formulated  as: 


1 


2  2 
n  n 

1=1  j=l 


The  total  damage  inflicted  to  all  targets  is  the  sum  of  damages  to  each 
target  given  above.  Each  target  can  be  given  a  military  target  value, 
u^,  which  weights  that  target  in  terms  of  its  military  importance.  Mili¬ 
tary  importance  could  be  an  estimate  of  value  in  terras  of  specific  numbers 
(i.e.,  number  of  planes  destroyed),  a  target  area  or  fraction  of  area, 
dollars,  or  some  other  measure. 


The  damage  function  then  is: 


5  2  2  y 

l  \  a  ~  n  n  a  i3k) 

k=l  i=l  j=l  JK 


Also  note  that  the  variable  and  y^^  are  restricted  to  non-negative 

values . 

The  problem  to  be  solved  is:  given  the  parameters  for  the  problem, 

determine  the  values  of  the  variables  x .  .  and  y. such  that  the  total 

2-3 

damage  function  is  maximized  subject  to  the  aircraft,  carrier,  and  mission 
constraints. 

This  problem  will  now  be  converted  into  a  standard  NLGP  formulation. 
The  problem  variables  must  be  converted  into  a  standard  format  for  the 
NLGP  model.  These  model  variables  are  listed  in  Table  2.1.  The  parameters 
necessary  to  formulate  the  model  objectives  are  listed  in  Tables  2.2  and 
2.3  and  Figure  2.1.  Negative  and  positive  deviation  variables  have  been 


added  to  the  objectives. 

Aircraft  availability  objectives 


Goal  1:  xi  +  x3  +  n]_  -  P]_  =  27 
Goal  2:  x2  +  x4  +  n2  ~  ?2  =  102 


Carrier  space  objectives 


Goal  3:  2920  ^  +  1770  x2  +  n3  -  p3  =  112,300 
Goal  4:  2920  x3  +  1770  x^  +  n^  -  p4  =  147,100 

Limitations  on  missions  due  to  aircraft  availability 


Goal  5:  x^  -  0.06452  (x^  +  x^  +  x^)  -  0.06250  (x^  +  Xg) 

+  n$  -  p5  -  0 

Goal  6:  x^  -  0.05556  (x^q  +  x^2  +  x^4)  ~  0.05264  (x^  +  x^) 
+  n6  -  p6  -  0 
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Goal  7 : 


Goal  8: 


x^  -  0.06896  “  0.06452  (x^  +  x^g) 

-  0.06250  (x^  +  x^g)  +  rij  -  =  0 

x4  -  0.05882  x2Q  -  0.05556  (xn  + 

-  0.05264  (x22  +  x23)  +  ng  -  p8  =  0 


Damage  objective 


+ 

+ 

+ 


40 

(1  -  0.99978 

(x5 

+ 

x15) 

0.99953 

(x10 

+ 

/“s 

/ — s 

o 

CM 

X 

10 

(1  -  0.99978 

(x6 

+ 

x16> 

0.99953 

(xll 

+ 

x21>) 

50 

(1  -  0.99978 

(x7 

+ 

x17) 

0.99953 

(x12 

+ 

x22)) 

ng  *  Pg  =  0 


Non-negativity  objectives 

x,  n,  p  ^  0  (the  zero  vector), 


Note  that  only  those  terms  with  non-zero  values  appear  in  the  damage 

objective.  Scenario  1  data  were  used. 

All  that  remains  to  be  done  is  to  formulate  the  achievement  vector. 
This  problem  has  two  priority  levels.  Goals  one  through  eight  compose 
priority  level  one.  These  are  absolute  objectives  (i.e.,  limited 
resources)  which  cannot  be  violated.  Goals  one  through  four  are  resource 
objectives.  These  cannot  be  exceeded.  Therefore,  the  positive  deviation 
variables  for  these  objectives  (p^  -  p^)  are  minimized  at  priority  level 
one.  This  prohibits  overachievement  of  these  objectives.  Goals  five 
through  eight  relate  aircraft  to  missions,  where  the  total  number  sent  on 
any  mission  can  equal  to  or  less  than  the  total  supply  but  can  never 
exceed  it.  Hence,  the  negative  deviation  variables  (n^  -  ng)  are  mini¬ 
mized.  This  prohibits  underachievement  of  these  objectives. 
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Table  2.1  Target  Allocation  Problem  Variables 


Description 


Problem 

Variable 


NLGP  Model 
Variable 


Distribution  of  aircraft 
of  type  j  to  carrier  i 


Missions  flown 
from  carrier  1 
to  target  k 


Using 
aircraft 
of  type  1 


Using 
aircraft 
o  f  type  2 


Missions  flown 
from  carrier  2 
to  target  k 


Using 
aircraft 
of  type  1 


Using 
aircraft 
of  type  2 


y211 

X15 

y212 

x16 

y213 

X17 

y214 

X18 

y215 

X19 

y221 

X20  . 

y222 

X21 

y223 

CM 

CM 

X 

y224 

x23 

y225 

X24 

Table  2.2  Target  Allocation  Problem  Parameters  [28] 


Type 

Value  for  Model 

ajk 

an  =  0.99978 

a2  =  0.99953 

a12  =  0.99978 

a22  =  0.99953 

a13  =  0.99978 

a23  -  0.99953 

a  .  =  0.99935 

14 

a24  =  0.99960 

a15  -  0.99935 

a25  »  0.99960 

b1  =  27 

(A-6  Aircraft) 

b. 

J 

b2  =  102 

(A-7  Aircraft) 

c3  =  2,920  sq.  ft. 

(A-6  Aircraft) 

c . 

J 

c2  *  1,770  sq.  ft. 

(A-7  Aircraft) 

di 

d^  *  112,300  sq.  ft 
d2  ■  147,100  sq.  ft 

(CVA) 

(CVAN) 

f 

flu  =  0.06452 

f2n  =  0.06896 

ijk 

f112  =  0.06250 

f212  -  0.06452 

f113  *  0.06452 

f213  -  0.06250 

fliy  =  0.06250 

114 

f01/  =  0.06250 

f115  =  °-06452 

f215  =  °'06452 

f121  =  °*05556 

f221  «  0.05882 

f122  "  °* 05264 

f222  =  °'05556 

f123  -  0.05556 

f223  -  0.05264 

f124  «  0.05264 

f224  =  0.05264 

f125  =■  0.05556 

f225  =  0.05556 

;e:  f. ..  values  were  determined  by  distance  to  the  target 
^  and  required  maintenance  between  missions. 
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Table  2.3  Military  Target  Values  [28] 


Target  Set 

Scenario  1 

- 

Scenario  2 

Scenario  3 

1. 

Airfields 

40 

45 

0 

2. 

Fuel  Storage 

Depots 

10 

15 

0 

3. 

Close  Support 

50 

40 

45 

4. 

Short  Interdiction 

0 

0 

45 

5. 

Command  and 

Control 

0 

0 

10 

•mmm 


Figure  2.1  Sample  scenario  for  target 
allocation  problem  [28] 
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The  second  priority  level  contains  only  one  objective,  the  damage 
objective.  The  aspiration  level  for  Goal  9  is  set  equal  to  zero,  and 
since  this  objective  is  maximized,  underachievement  of  the  objective 
is  minimized  by  minimizing  the  negative  deviation  variable  (ng).  Any 
overachievement  will  then  be  indicated  by  a  non-zero  value  for  Pg,  the 
measure  of  total  damage  inflicted  on  all  targets. 

The  NLGP  model  for  this  problem  is: 

Find  the  values  of  x  =  (x^,  x? . 

so  as  to 

minimize  a  =  (Px  +  P2  +  P3  +  P4  +  n5  +  ng  +  n?  +  ng,  ng) 
subject  to 

Goal  1  through  Goal  9 
and  x,  n,  p  >_  0. 

The  discussion  of  the  solution  to  this  target  allocation  problem 
will  be  deferred  until  later  chapters. 

2.3  Advantages  and  Limitations 

There  is  still  no  single  best  method  for  solving  nonlinear  models, 
whether  they  contain  a  single  objective  or  multiple  objectives.  Each 
method  has  advantages  and  disadvantages  inherent  to  that  method.  Some 
advantages  and  limitations  of  the  NLGP  method  are  reviewed  here  [18]. 

The  greatest  advantage  of  the  NLGP  method  over  conventional  methods 
is  that  this  method  allows  the  analyst  or  design  engineer  to  model  the 
problem  realistically.  Most  problems  have  more  than  a  single  objective, 
and  generally,  some  of  these  objectives  will  be  conflicting  in  nature. 
The  objectives  of  the  system  do  not  have  to  be  massaged  to  fit  into  some 
prescribed  mold  so  that  a  solution  can  be  obtained. 
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A  major  advantage  of  the  NLGP  method  over  conventional  methods  is 
that  a  solution  to  the  problem  is  always  obtained.  This  is  due  to  the 
direct  search  technique  used  in  the  method  and  the  ordered  set  of  pre¬ 
emptive  priority  levels.  Rather  than  considering  infeasibility,  solu¬ 
tions  are  either  implementable  (i.e.,  a^  =  0)  or  un implement able 
(i.e.,  a^  >  0).  The  achievement  vector  indicates  the  degree  to  which 
the  desired  aspiration  levels  are  met. 

In  general  in  nonlinear  optimization,  there  is  no  way  to  guarantee 
that  the  solution  obtained  will  be  a  global  optimum.  The  NLGP  technique 
performs  no  better  or  worse  than  other  methods  in  this  respect.  In 
actual  practice,  this  technique  has  proved  to  be  a  robust  method  for  the 
solution  of  some  very  complex  problems. 

The  priority  structure  does  have  limitations,  however.  The  existence 
of  a  preemptive  priority  structure  may  not  always  be  a  valid  assumption. 
Absolute  objectives  (i.e.,  natural  laws,  limited  resources)  cannot  be 
violated  and  belong  in  the  highest  priority  level.  The  number  and  arrange¬ 
ment  of  objectives  in  lower  priority  levels  may  not  be  as  obvious  and  can 
depend  upon  the  analyst,  or  engineer’s  understanding  and  insight  into  the 
problem.  Quite  often  some  objectives  must  be  switched  between  priority 
levels  or  entire  levels  rearranged  to  obtain  satisfactory  results,  espe¬ 
cially  when  little  is  known  apriori  about  the  problem. 

The  use  of  weighting  factors  in  the  achievement  vector  is  an  advan¬ 
tage  because  these  factors  enable  the  analyst  or  design  engineer  to 
weight  objectives  only  within  a  single  priority  level.  The  proper 
weighting  of  objectives  is  not  a  clearcut  procedure,  however,  and  these 
weights  often  must  be  established  through  an  iterative  process. 


Finally,  aspiration  levels  for  some  objectives  may  not  always  be 
obvious.  Again,  proper  determination  of  these  levels  may  have  to  be 
accomplished  through  a  number  of  trials  and,  in  any  case,  evaluated  in 
the  final  result  via  a  sensitivity  analysis. 
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CHAPTER  III 

COMPUTATIONAL  TECHNIQUES 


3.1  Introduct ion 

A  variety  of  methods  can  be  employed  to  solve  nonlinear  models 
[5,6,7,8,9,10,11,12,27,30].  The  success  of  a  method  over  other  methods 
to  find  a  solution  may  be  relative  to  the  particular  problem  as  well  as 
to  initial  problem  parameters.  A  straightforward  approach  for  solving 
large-scale  nonlinear  models  is  one  employing  direct  search  techniques. 

The  method  selected  for  the  nonlinear  goal  programming  technique  is 
based  on  an  extension  of  the  search  method  of  Hooke  and  Jeeves  [12]. 

The  Hooke-Jeeves  method  is  an  accelerated  pattern  search  technique. 

It  is  based  on  the  assumption  that  any  set  of  moves  which  has  been  suc¬ 
cessful  in  optimizing  the  objectives  is  worth  repeating.  The  algorithm 
starts  with  an  initial  point  (i.e.,  an  n-dimensional  vector),  perturbs 
each  variable  one  at  a  time  and  tests  for  an  Improved  solution.  When  a 
search  has  been  successful,  future  increments  for  the  pattern  search 
moves  are  increased.  When  a  search  is  unsuccessful  in  finding  an  improved 
solution,  future  increments  are  decreased.  This  pattern  search  technique 
has  been  shown  to  be  a  simple  and  effective  method  for  solving  actual  prob¬ 
lems  with  large  numbers  of  variables,  large  number  of  objectives,  and  a 
variety  of  convergence  criteria  [12].  It  is  a  simple  method  which  makes 
it  computationally  efficient.  This  makes  it  attractive  for  use  with  the 
NLGP  algorithm. 
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3.2  The  Nonlinear  Goal  Programming  Computer  Code 

The  software  package  developed  under  this  contract  contains  three 
separate  computer  codes.  All  three  codes  can  be  used  to  solve  an  NLGP 
problem.  All  three  codes  are  very  similar  in  technique.  Each  code  has 
some  minor  modifications  which  makes  it  better  suited  to  solve  a  differ¬ 
ent  variation  of  NLGP  problem.  The  basic  computer  code  and  its  operating 
characteristics  will  be  discussed  in  detail  in  this  section.  The  remain¬ 
ing  two  codes  will  be  discussed  in  the  section  covering  modifications. 

All  program  versions  are  written  in  American  National  Standard  FORTRAN  IV. 

The  basic  computer  code  is  referred  to  as  the  nonlinear  goal  program¬ 
ming/modified  pattern  search  -  ridge  search  code,  NLGP/MPS-RS  code.  This 
is  Version  1  of  the  code.  A  complete  listing  of  this  code  is  contained 
in  Appendix  A. 

The  NLGP/MPS-RS:  Version  1  code  is  composed  of  a  small  main  program 
and  seven  subprograms.  Two  of  these  seven  subprograms  must  be  supplied  in 
part  by  the  user  for  each  new  NLGP  problem.  Figure  3.1  is  a  schematic  dia¬ 
gram  showing  the  operating  relationships  of  the  subprograms. 

The  MAIN  program  contains  data  initialization  statements,  type  declara¬ 
tion  statements,  and  other  specification  statements  that  are  necessary  for 
the  definition  of  the  program  parameters  and  variables.  The  MAIN  program 
calls  subroutine  DATAIN  and  subroutine  HJALG,  the  Hooke-Jeeves  pattern 
search  algorithm.  Subroutine  DATAIN  is  used  to  input  all  problem  data. 

Once  this  is  accomplished,  all  further  control  of  the  program  proceeds 
through  subroutine  HJALG.  Table  3.1  gives  a  brief  description  of  the 
modules  making  up  the  computer  code.  Module  length  is  given  as  compiled 
on  an  IBM  370/3033  System  using  the  FORTRAN  WATFIV  compiler. 


ACHFUN 


DECIDE 


DEWAR 


OBJFUN 


Figure  3.1  Program  Relationship  Diagram  - 
NLGP/MPS-RS  Code:  Versions  1 

and  ? 


Table  3.1  Program  Module  Information  NLGP/MPS-RS  Code: 
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Each  subprogram  of  the  code  will  be  discussed  in  terms  of  its  opera¬ 
tion  and  the  program  parameters  and  variables  which  it  uses.  The  complete 
FORTRAN  language  listing  of  each  subprogram  is  contained  in  Appendix  A.  A 
list  of  all  program  variables  and  parameters  is  given  in  Table  3.2  A  com¬ 
plete  guide  on  how  to  select  the  values  of  program  parameters  is  contained 
in  the  User's  Manual. 

The  computer  program  is  presently  dimensioned  to  handle  problems  with 
up  to  2500  objectives,  2500  variables,  and  10  priority  levels.  These 
dimensions  can  be  altered  to  increase  the  problem  size  which  can  be  solved. 


Subroutine  DATAIN 


Subroutine  DATAIN  is  used  to  input  all  the  necessary  data  for  the 

NLGP  problem.  Once  this  data  is  read  into  internal  storage,  an  echo  check 

of  this  data  is  printed  as  a  means  of  checking  it.  Table  3.3  contains  a 

listing  of  the  data  which  must  be  input  and  its  form. 

Variables  and  parameters  used  in  DATAIN: 

NOBJ 

NPRIOR 

NVAR 

NMAX 

EPS 

ACCEL 

REDUCE 

DELTA 

X 

LBD 

UBD 

RHS 

EPSY. 


Table  3.2  Computer  Code  Variables  and  Parameters 
NLGP/MPS-RS  Code:  Version  1 


Name 

NOBJ 

NPRIOR 

NVAR 

NMAX 

NCYCLE 

X 

XBASE 

A 

ABEST 

APOS,  ANEG 

EPS,  IPS 

LBD,  UBD 
OBJ 
RHS 
N,  P 
EPSY 

ACCEL 

REDUCE 

DELTA 

IMPROV 


Description 

Number  of  problem  objectives 

Number  of  priority  levels  in  the  achievement  vector 
Number  of  decision  variables 

Maximum  number  of  pattern  search  cycles  allowed  before 
termination;  set  by  user 

Counter  of  pattern  search  cycles 

Decision  variable  vector  for  perturbations  in  pattern 
search 

Decision  variable  base  point  vector  for  pattern  search 
Achievement  vector  at  a  test  point 

Achievement  vector  for  best  solution  during  pattern 
search 

Temporary  achievement  vectors  for  the  resolution  ridge 
search 

Perturbation  step-size  scalars  used  for  exploratory 
moves 

Lower  and  upper  bound  vectors  on  decision  variables 

Objective  function  vector 

Desired  aspiration  levels  for  objectives 

Negative  and  positive  deviation  vectors 

Test  vector  for  termination  of  algorithm  by  meeting 
specified  tolerances 

Acceleration  factor  for  pattern  moves 

Step-size  reduction  factor 

Minimum  allowable  step-size  for  exploratory  moves 

A  test  switch;  if  IMPROV  =*  1,  X  improves  the  solution 
for  ABEST,  if  IMPROV  =  0,  X  does  not  improve  this 
solution 


n 


4 
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Name 

IPRINT 

ITERM 

IPROG 

IACCEL 

SAVE 

I 

J 

K 


Table  3.2  Continued 
Description 

An  output  switch;  if  IPRINT  *  1,  complete  results  are 
printed  out  at  every  pattern  search  cycle,  if  IPRINT  =* 
0,  only  the  final  results  are  printed. 

A  termination  testing  indicator;  if  ITERM  *  1,  the 
solution  for  A  is  compared  to  ABEST  and  then  tested 
for  termination  against  EPSY,  if  ITERM  -  0,  only  a 
comparison  of  A  and  ABEST  is  made. 

A  counter  for  number  of  improvements,  if  any,  during 
the  exploratory  moves 

An  indicator  for  pattern  moves;  if  IACCEL  *  1,  the 
pattern  moves  was  successful,  if  IACCEL  =*  0,  no 
pattern  move  was  attempted,  if  IACCEL  =  -1,  the 
pattern  move  was  unsuccessful 

A  temporary  storage  variable  for  pattern  and  ridge 
moves 

Counter  for  all  do  loops  containing  the  objective 
functions,  OBJ 

Counter  for  all  do  loops  containing  the  priority 
levels  of  the  achievement  vector,  A  and  ABEST 

Counter  for  all  do  loops  containing  the  decision 
variables,  X  and  XBASE 
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Table  3.3  Data  Input  Guide  NLGP/MPS-RS 
Code :  Version  1 


Card 

Data 

Type 

Format 

1 

NOBJ 

Integer  *4 

15 

NPRIOR 

n 

15 

NVAR 

it 

15 

NMAX 

it 

15 

IPRINT 

Integer  *2 

15 

,  .NVAR 

8 

1) 

X 

Real  *4 

8F10.0 

,  NVAR 

J  1  8 

1) 

LBD 

Real  *4 

8F10.0 

.  .NVAR 

4  -  (  g  - 

1) 

UBD 

Real  *4 

8F10.0 

.NOBJ 

K  8 

1) 

RHS 

Real  *4 

8F10.0 

,  ,NPRIOR 

6  -  (  8 

-  1) 

EPSY 

Real  *4 

8F10.0 

7 


EPS 

ACCEL 

REDUCE 

DELTA 


Real  *4 
11 

»1 


F10.0 

F10.0 

F10.0 

F10.0 


11 
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Subroutine  HJALG 


Subroutine  HJALG  is  the  Hooke-Jeeves  pattern  search  Algorithm  as 
modified  to  solve  NLGP  problems.  Figure  3.2  is  a  schematic  diagram 
showing  the  flow  through  the  pattern  search  algorithm.  A  detailed  com¬ 
puter  logic  flow  diagram  is  contained  in  Figure  3.3. 

Variables  and  parameters  used  in  HJALG: 


NOBJ 

IPRINT 

NPRIOR 

ITERM 

NVAR 

IPROG 

UMAX 

IACCEL 

NCYCLE 

IPS 

X 

EPS 

XBASE 

LBD 

A 

UBD 

ABEST 

IMPROV 

SAVE 

ACCEL 

REDUCE 

DELTA 

Subroutine  RIDGE 


Subroutine  RIDGE  is  used  to  calculate  a  new  temporary  head  point  when 
the  exploratory  moves  of  the  pattern  search  fail  to  find  any  improved  points 
near  the  base  point,  XBASE.  RIDGE  is  a  heuristic  technique  which  uses  infor¬ 
mation  for  each  set  of  pairwise  perturbations,  positive  and  negative  around 
a  variable  x(k).  It  attempts  to  move  in  a  direction  which  is  oblique  to 
the  original  axes  of  exploration,  and  therefore,  move  away  from  the  present 
base  point  around  which  the  pattern  has  collapsed. 

The  basis  for  RIDGE  comes  from  Rosenbrock's  technique  for  rotation  of 
coordinates  [29]  and  Mugele's  "poor  man's  optimizer”  [26].  RIDGE  requires 
much  less  storage  of  information  than  Rosenbrock's  method.  It  is  construc¬ 
ted  specifically  to  use  information  from  the  achievement  vector  to  establish 
a  direction  for  continued  search.  It  is  called  everytime  the  step  size,  IPS, 
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Figure  3.3  Continued 


becomes  smaller  in  value  than  the  minimum  allowable  step  size,  DELTA. 


Figure  3.4  shows  the  flow  diagram  for  subroutine  RIDGE. 
Variables  and  parameters  used  in  RIDGE: 


X 

NPRIOR 

XBASE 

SAVE 

IPS 

NVAR 

EPS 

ITERM 

LBD 

IMPROV 

UBD 

ACCEL 

APOS 

ANEG 

A 

ABEST 

Subroutine  DECIDE 

Subroutine  DECIDE  has  two  main  purposes:  (1)  To  compare  the  achieve¬ 
ment  vector  at  a  test  point.  A,  against  the  achievement  vector  for  the  best 
point,  ABEST.  If  A  is  preferred  over  ABEST,  it  becomes  ABEST  and  the  cor¬ 
responding  x  vector,  X,  becomes  XBASE,  the  new  base  point.  If  A  Is  not 
preferred  over  ABEST,  the  program  continues  with  no  changes  to  the  variables, 
and  (2)  To  test  the  achievement  vector  for  termination  of  the  algorithm. 
Termination  can  occur  by  two  methods:  first,  all  elements  of  the  achieve¬ 
ment  vector  have  zero  values;  this  is  an  optimal  solution,  and  secondly, 
when  all  elements  of  the  achievement  vector  have  values  less  than  the 
specified  tolerance  vector,  EPSY.  This  is  the  test  by  tolerances.  Figure 
3.5  is  the  flow  diagram  for  subroutine  DECIDE. 

Variables  and  parameters  used  in  DECIDE: 


A 

ABEST 

NPRIOR 

IMPROV 

ITERM 

X 
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Subroutine  RIDGE 


5 

Figure  3.4  Flow  Diagram  for  Subroutine  Ridge 


Figure  3.4  Continued 


Subroutine  DECIDE 


Figure  3.5  Flow  Diagram  for  Subroutine  DECIDE 
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A  and  X 
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Subroutine  ACHFUN 

Subroutine  ACHFUN  is  the  module  which  contains  the  functions  of  weighted 
positive  and  negative  deviation  variables.  These  functions  must  be  supplied 
by  the  user  for  each  new  NLGP  problem.  An  example  will  be  presented  in  the 
next  chapter  which  will  show  these  functions,  one  for  each  priority  level  in 
the  achievement  vector.  Figure  3.6  shows  the  computer  flow  diagram  for  this 
subroutine. 

Variables  and  parameters  used  in  ACHFUN: 

A 

N 

P 

Subroutine  DEWAR 

Subroutine  DEWAR  calculates  the  values  of  the  negative  and  positive 
deviation  variables.  It  uses  the  values  of  the  objective  functions  already 
calculated  and  compares  these  values  to  the  desired  aspiration  levels,  RHS. 
Figure  3.7  shows  the  flow  diagram. 

Variables  and  parameters  used  in  DEWAR: 

NOBJ 

OBJ 

RHS 

N 

P 

Subroutine  OBJFUN 

Subroutine  OBJFUN  is  used  to  calculate  the  values  of  the  objective 
functions  at  the  test  point  x.  These  functions  must  be  user  supplied  for 
each  new  NLGP  problem  to  be  solved.  An  example  of  this  subroutine  will  be 
presented  in  the  next  chapter.  A  flow  diagram  of  this  subroutine  is  con¬ 
tained  in  Figure  3.8. 

Variables  and  parameters  used  in  OBJFUN: 

OBJ 

X 


2 


Figure 


3.6  Flow  Diagram  for  Subroutine  ACHFUN 
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Subroutine  DEWAR 


Figure  3.7  Flow  Diagram  for  Subroutine  DEWAR 


Subroutine  OBJFUN 


Figure  3.8  Flow  Diagram  for  Subroutine  OBJFUN 
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3.3  Modifications 

The  second  and  third  versions  of  the  NLGP/MPS-RS  code  are  modifica¬ 
tions  of  version  1.  The  modifications  which  make  them  different  from  the 
basic  code  are  presented  here.  Otherwise,  these  codes  are  exactly  the 
same  as  version  1.  Complete  FORTRAN  language  computer  listings  of  these 
versions  are  contained  in  Appendices  B  and  C. 

The  NLGP/MPS-RS:  Version  2  Code 

This  computer  code  is  exactly  like  version  1  except  that  the  perturba¬ 
tion  step  sizes,  EPS  and  IPS,  are  vectors  instead  of  scalars.  This  enables 
the  user  to  input  a  separate  step  size  for  each  variable.  This  can  be 
especially  important  when  the  range  of  values  which  the  variables  can  take 
on  vary  considerably. 

The  major  difference  in  the  code  from  version  1  is  in  subroutine  DAT AIN . 
Now  a  vector  of  values  for  EPS  must  be  read  into  storage.  Table  3.4  lists 
the  necessary  cards  for  input  for  version  2. 

Minor  modifications  also  appear  in  subroutines  HJALG  and  RIDGE.  All 

statements  previously  using  IPS  and  EPS  now  use  the  variable  forms,  IPS(k) 

and  EPS(k) .  The  resetting  of  the  step  size  vector  in  subroutine  RIDGE  is 

accomplished  using  a  do  loop  of  the  form:  DO  10  K*l,  NVAR 

IPS(k)  =  EPS(k) . 

Advantages  for  using  this  version  as  compared  to  the  other  versions' 
will  be  presented  in  the  next  section. 

The  NLGP/MPS-RS:  Version  3  Code 

Version  3  of  the  computer  code  is  the  same  as  version  2  with  an  addi¬ 
tional  subroutine,  subroutine  UVALUE.  This  subroutine  has  been  added  to 
reduce  the  computation  time  necessary  for  the  evaluation  of  the  objective 


functions. 


mmmmm 


'VrTT^V, 
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Table  3.4  Data  Input  Guide  NLGP/MPS-RS 
Code:  Version  2 


Card(s) 

Data 

lXE£ 

Format 

1 

NOBJ 

Integer  *4 

15 

NPRIOR 

II 

15 

NVAR 

II 

15 

NMAX 

ft 

15 

IPRINT 

Integer  *2 

15 

2  -  -  1) 

X 

Real  *4 

8F10.0 

3  -  (^f  -  1) 

EPS 

Real  *4 

8F10.0 

4  -  (™  -  1) 

LBD 

Real  *4 

8F10.0 

5  -  rf  - 1) 

UBD 

Real  *4 

8F10.0 

6  -  (N0B8J  -  1) 

RHS 

Real  *4 

8F10.0 

_  .NPRIOR  - . 

Real  *4 

7  (  8  1) 

EPSY 

8F10.0 

8 

ACCEL 

Real  *4 

F10.0 

REDUCE 

II 

F10.0 

DELTA 

II 

F10.0 

5 


. -—-I  -  II I . . ' . . .  II  I f,iifi II-  ur****~~*****r - 


-  .  1 
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Whenever  a  new  test  point,  X,  is  formed,  subroutine  OBJFUN  is  even¬ 
tually  called  to  evaluate  all  the  objectives  at  this  test  point  so  that 
a  new  achievement  vector  can  be  formed.  There  are  many  cases,  however, 
when  it  is  not  necessary  to  re-evaluate  all  the  objective  functions. 

There  are  two  reasons  for  this: 

(1)  Many  problems  have  objectives  which  are  sparsely  populated  by 
the  variables.  It  is  not  necessary  to  re-calculate  all  the 
objectives  when  the  variable  being  perturbed  only  appears  in 
a  few  of  the  objectives,  and 

(2)  Since  the  pattern  search  algorithm  perturbs  one  variable  at  a 
time  during  exploratory  moves,  it  is  not  necessary  to 
re-calculate  every  term  in  each  objective  function.  Rather, 
only  those  terms  in  which  the  perturbed  variable  appears  need 
to  be  recalculated,  and  the  net  effect  on  the  value  of  the 
objective  function  can  be  determined. 

Both  categories  occur  frequently  in  large-scale  nonlinear  problems.  The 
use  of  the  UVALUE  subroutine  eliminates  many  unnecessary  calculations . 

This  results  in  significant  reductions  in  computation  time,  especially 
in  large  problems.  The  subroutine  can  be  especially  useful  during  the 
exploratory  moves. 

Additionally,  this  subroutine  is  problem  dependent  and  must  be 
supplied  in  part  by  the  user  for  implementation  with  the  code.  A  flow 
chart  description  of  subroutine  UVALUE  is  given  here  for  the  target 
allocation  problem  described  in  Chapter  II.  A  complete  description  of 
the  use  and  set-up  of  this  subroutine  is  contained  in  the  User's  Manual. 
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The  major  modifications  contained  in  this  code  will  now  be  presented. 
The  new  parameters  and  variables  which  are  contained  in  version  3  are 
listed  in  Table  3.5.  There  use  in  the  program  will  be  indicated  through 
flow  charts  showing  the  revised  subroutines. 

Subroutine  DATAIN 

Subroutine  DATAIN  is  exactly  the  same  as  that  used  in  version  1  with 
one  addition.  The  matrix  of  coefficients,  AMAP<  ’ ,  L '  must  be  read  into 
storage.  It  is  currently  dimensioned  to  (100,100).  Larger  problems  could 
be  solved  by  increasing  this  matrix  in  subroutine  LVALUE  and  DATAIN. 

Table  3.6  shows  the  revised  data  input  guide  for  version  3. 

Subroutine  HJALG 

Subroutine  HJALG  in  version  3  differs  from  that  in  version  ±  in  that 
it  now  calls  subroutine  UVALUE  instead  of  subroutine  OBJFUN  in  many  places 
in  the  program,  as  indicated  by  the  switch,  JPI.  Subroutine  OBJFUN  is 
called  initially  to  calculate  the  values  of  the  objective  functions. 

During  the  exploratory  moves,  subroutine  UVALUE  is  used  to  calculate  the 
incremental  changes  in  the  objective  functions  for  the  variable  being 
perturbed.  The  additional  variables  and  parameters  used  in  subroutine 
HJALG  are: 

JPI 

LNUM 

UVALUE 

SEPS. 

The  complete  flow  diagram  for  this  subroutine  is  in  Figure  3.9. 

Subroutine  RIDGE 

Same  as  in  version  1. 

Subroutine  DECIDE 


Same  as  in  version  1. 
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Table  3. 

Name 

LNVM 

DIFF 

JPI 

AMAP 

SEPS 

A ,  B ,  C ,  D ,  E 


5  Additional  Parameters  and  Variables 
for  NLGP/MPS-RS  Code:  Version  3 


Description 

Counter  used  to  indicate  which  variable 
K  is  being  perturbed  in  the  exploratory 
moves 


A  vector  of  the  incremental  changes  in 
the  objective  functions 


A  switch  indicator:  if  JPI  =  0  subrou¬ 
tine  UVALUE  is  not  used,  if  JPI  =  1, 
subroutine  UVALUE  is  called  to  calcu¬ 
late  DIFF  for  the  objectives 


A  matrix  of  coefficients  for  the 
objective  functions 


An  adjusted  step  size  scalar 


Subroutine  UVALUE  internal  variables: 
problem  dependent  on  actual  use 
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Table  3.6  Data  Input  Guide  NLGP/MPS-RS 
Code:  Version  3 


Card(s) 

Data 

TyPe 

Format 

1 

NOBJ 

Integer  #4 

15 

NPRIOR 

II 

15 

NVAR 

II 

15 

NMAX 

II 

15 

IPRINT 

II 

15 

IU 

11 

15 

,  ,NVAR 

''8 

1) 

X 

Real  // 4 

8F10.0 

-  ,NVAR 

J  V  8 

1) 

LBD 

Real  *4 

8F10.0 

.  .NVAR 

4  -  (  8  - 

1) 

UBD 

Real  *4 

8F10.0 

.NOBJ 

5  K  8 

1) 

RHS 

Real  *4 

8F10.0 

,NPRIOR 
6  “  (  8 

-  1) 

EPSY 

Real  *4 

8F10.0 

7  - 

EPS 

Real  *4 

F10.0 

ACCEL 

II 

F10.0 

REDUCE 

II 

F10.0 

DELTA 

II 

F10.0 

(NVAR x  NOBJ 

8  g  "  1) 

AMAP 

Real  *4 

8F10.0 

i 


Figure  3.9  Flow  Diagram  for  Subroutine  HJALG: 
Version  3 


Figure  3.9  Continued 


Figure  3.9  Continued 
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Figure  3.9  Continued 
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Figure  3.9  Continued 


Figure  3 . 9  Continued 
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Subroutine  ACHFUN 

Same  as  in  version  1  (problem  dependent) . 

Subroutine  DEWAR 

This  subroutine  is  slightly  different,  in  that  it  now  tests  to  see 
if  subroutine  OBJFUN  is  UVALUE  is  to  be  called,  then  calls  the  appropriate 
one.  Figure  3.10  shows  the  flow  diagram  for  this  subroutine. 

Subroutine  OBJFUN 

Same  as  in  version  1  (problem  dependent). 

Subroutine  UVALUE 

This  subroutine  calculates  the  incremental  changes  in  the  objective 
functions.  The  variables  and  parameters  used  in  this  subroutine  are: 


NOBJ 

SEPS 

L 

A,B,C,D,E 

LNUM 

Y9 

DIFF 

AMAP 

Calls  to  other  subroutines:  None. 

The  flow  diagram  for  this  subroutine  is  contained  in  Figure  3.11.  It 
should  be  restated  that  this  subroutine  has  many  statements  which  are  prob¬ 
lem  dependent,  and  those  shown  in  the  flow  diagram  represent  a  schematic 
for  the  objectives  used  in  the  target  allocation  problem  of  Chapter  II. 

A  more  thorough  discussion  of  the  formulation  of  this  subroutine  is  given 
in  the  User's  Manual. 

Further  discussion  of  the  three  versions  and  their  respective  advan¬ 
tages  and  limitations  will  be  given  in  the  following  sections  of  Chapter 


III. 


Subroutine  UVALUE 


Figure  3.11  Flow  Diagram  for  Subroutine  UVALUE 


Figure  3.11  Continued 


3.4  Specifications 


Each  version  of  the  code  has  different  computer  storage  specifications, 

due  to  the  modifications  contained  in  versions  2  and  3.  All  three  codes  are 

presently  dimensioned  for  problems  of  the  following  maximum  size: 

Number  of  objectives  =  2500 

Number  of  variables  =  2500 

Number  of  priority  levels  =  10  . 

The  computer  core  storage  necessary  for  each  version  of  the  code  will 
be  dependent  upon  the  number  and  type  of  objective  functions  and  the  number 
of  and  term  in  the  achievement  vector.  If  every  variable  appears  in  every 
objective  function,  the  core  storage  necessary  for  subroutine  OBJFUN  could 
be  conservatively  estimated  to  be: 

(80  Bytes/card)  x  (1  card/10  variables)  x  (2500  variables) 
=  20,000  Bytes/objective 

(20,000  Bytes/objective)  x  (2500  objectives) 

=  50  million  Bytes. 

This  is  an  estimated  minimum  storage  requirement  i_f  all  variables  appear  in 
all  objectives.  Since  most  problems  do  not  have  this  form,  actual  storage 
will  be  much  smaller.  Table  3.7  gives  the  requirements  which  are  minimum 
values,  and  actual  requirements  will  vary  greatly  depending  on  the  problem. 

3.5  Advantages  and  Limitations 

There  are  advantages  to  using  each  code  which  depend  on  the  type  of 
problem  to  be  solved.  Version  1  is  the  basic  code,  the  simplest  in  form, 
and  the  easiest  to  use.  This  code  is  differentiated  from  the  other  two 
versions  in  that  it  uses  a  scalar  for  the  step  size,  EPS,  for  the  pertur¬ 


bations  in  the  exploratory  moves.  This  is  advantageous  when  the  decision 


I 

I 
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Table  3. 7  Storage  Requirements 


Bytes  (in 

decimal) 

Version 

Object  Code 

Array  Area 

1 

22,016 

80,260 

2 

22,016 

90,260 

3 

22,784 

130,308 

r 
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variables  for  the  problem  are  all  similar  in  nature  and  in  the  range  of 
values  which  they  can  have.  Using  the  scalar  step  size  saves  10,000 
bytes  of  core  storage.  This  area  can  be  used  to  redimension  the  number 
of  variables  to  sizes  larger  than  2500. 

Version  1  is  limited  in  this  same  respect  if  the  range  of  values 
for  the  variables  vary  greatly  or  if  their  very  nature  makes  it  difficult 
to  select  only  a  single  step  size  for  the  exploratory  moves.  This  dis¬ 
advantage  is  eliminated  by  the  use  of  version  2.  With  a  vector  of  step 
sizes,  each  variable  can  have  a  distinct  step  size  to  account  for  the 
type  and  range  of  values  peculiar  to  it.  This  results  in  an  additional 
cost  of  10,000  bytes  of  core  storage,  but  only  a  negligable  increase  in 
computation  time. 

The  third  version  of  the  code  is  the  most  complex,  and  therefore,  the 
most  difficult  for  the  user  to  implement.  All  codes  must  have  user  supplied 
subroutines  OBJFUN  and  ACHFUN,  but  version  3  also  requires  a  user  supplied 
subroutine  UVALUE.  This  is  the  most  difficult  subroutine  to  formulate.  Its 
advantages  of  saving  considerable  computer  time  offsets  the  disadvantages  of 
its  additional  storage  requirements  and  difficulty  of  use  only  when  the  problem 
to  be  solved  is  large  (i.e.,  in  the  range  of  1000  or  more  variables  and  500 
or  more  objectives)  and  when  only  a  few  variables  appear  in  each  objective. 

It  is  the  opinion  of  the  authors  that  the  use  of  version  3  be  attempted  only 
after  the  user  has  gained  experience  using  versions  1  and  2,  and  also  has  a 
thorough  understanding  of  the  problem  to  be  solved. 

The  direct  search  technique  utilized  in  each  version  is  superior  to 


I 


gradient  search  techniques  in  the  number  and  type  of  evaluations  which  must 
be  made  to  calculate  an  achievement  vector.  No  gradients  have  to  be  calcu¬ 
lated  or  estimated  through  expansion  of  series. 
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This  technique  is  limited  in  that  one  can  only  hope  to  obtain  a  local 
optimal  for  the  problem  solution,  especially  for  large  and  complex  problems. 
(Global  optimal  solutions  can  be  obtained,  however,  when  the  problem  is  of 
a  special  nature  and  limited  complexity.) 

Finally,  the  computer  codes  are  designed  to  find  solutions  which  are 
continuous  -  valued  in  nature.  Problems  which  require  solutions  with  0-1 
valued  variables  and  mixed  integer-continuous  problems  present  special  dif¬ 
ficulties.  No  guarantee  can  be  made  that  the  solutions  obtained  will  be 
integer  valued.  Some  attempt  can  be  made  to  control  the  variable  values 
by  restricting  the  lower  and  upper  bounds  and  restricting  the  step  size 
used  for  variable  perturbations. 

More  comments  on  these  considerations  will  be  made  in  the  final 
chapter. 

3 . 6  Summary 

Three  different  versions  of  the  NLGP/MPS-RS  computer  code  have  been 
developed.  Each  version  has  advantages  over  the  others'  which  makes  it 
superior  for  solving  a  particular  type  or  class  of  NLGP  problems. 

A  brief  description  has  been  given  for  each  code,  explaining  its 
general  form  and  the  types  of  variables  and  parameters  used  in  each.  Flow 
diagrams  depict  the  subroutines  composing  each  code,  and  an  overall  flow 
schematic  shows  the  relationships  of  these  subroutines  to  form  an  entire 
computer  program.  Listings  for  each  are  given  in  Appendices  A,  B,  and  C. 

Some  general  advantages  and  limitations  for  the  use  of  each  version 
have  been  described.  More  specific  information  detailing  the  implementa¬ 
tion  of  each  code  can  be  found  in  the  User's  Manual  supplied  along  with 
this  report.  Considerations  on  problem  solution  and  computational  aspects 
will  be  given  in  Chapter  IV. 
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CHAPTER  IV 

COMPUTATIONAL  RESULTS 

4.1  Introduction 

The  general  nonlinear  goal  programming  problem  formulation  was  pres¬ 
ented  in  Chapter  II  along  with  an  example.  The  three  variations  of  the 
computational  techniques  used  to  solve  such  an  NLGP  problem  were  described 
in  Chapter  III.  The  solution  to  the  example  problem  will  be  presented 
here  first,  and  then  the  remainder  of  this  chapter  will  discuss  computa¬ 
tional  results  in  general. 

4.2  A  Target  Allocation  Problem:  Solution 

The  target  allocation  problem  formulated  in  Chapter  II  is  a  small- 
scale  problem,  but  it  shows  a  type  of  problem  suitable  to  naval  applica¬ 
tion.  The  NLGP  model  for  this  problem  is: 

Find  the  values  of  x  =  (x1,x2>  •••,  x^)  so  as  to  minimize 
a  =  (px  +  p2  +  p3  +  p^  +  n5  +  nfe  +  n7  +  ng,  ng)  subject  to: 

Aircraft  availability  objectives 


Goal 

1: 

x^  +  xg  +  nL  -  pL  = 

27 

Goal 

2: 

X2  +  X4  +  n2  '  p2  = 

102 

Carrier  space  objectives 

Goal  3:  2920  +  1770  x2  +  n3  -  P3  -  112,300 

2920  xg  +  1700  x^  +  n^  -  p^  = 


Goal  4: 


147,100 
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Limitations  on  missions  due  to  aircraft  availability 

Goal  5:  x^  -  0.6452  (x^  +  x^  +  x^)  -  0.06250  (xg  +  xg) 

+  n5  -  p5  =  0 

Goal  6:  x2  -  0.05556  (x^q  +  x^2  +  x^^)  -  0.05264  (x^  +  x^g) 

+  n6  -  p6  =  0 

Goal  7:  Xg  -  0.06896  x^g  -  0.06452  (x^g  +  x^g) 

-  0.06250  (x^  +  x^g)  +  n7  ~  P7  =  0 

Goal  8:  x^  -  0.05882  “  0.0556  (x^  +  x^) 

-  0.05264  (x22  +  x2g)  +  ng  -  p8  =  0 


Damage  objective 


Goal  9: 

40 

(1  -  0.99978 

(x5 

+ 

x15) 

0.99953 

(xio 

+ 

X2oh 

+ 

10 

(1  -  0.99978 

(x6 

+ 

X16} 

0.99953 

(xll 

+ 

x2ih 

+ 

50 

(1  -  0.99978 

(x? 

+ 

X17) 

0.99953 

(X12 

+ 

X22?  ) 

and  x,  n,  p  0 

This  problem  has  the  following  specifications  for  the  computer  codes: 

NOBJ  =  9 

NVAR  =  24 

NPRIOR  =  2  . 

The  objectives  are  divided  into  two  groups  for  the  priority  levels.  Priority 
level  one  contains  the  first  eight  objectives.  These  are  resource  objectives 
which  cannot  be  violated.  The  second  priority  level  measures  the  damage 
inflicted  on  all  targets  by  the  aircraft  from  the  two  aircraft  carriers. 

The  reader  is  referred  to  Chapter  II  for  a  complete  description  of  all 


problem  parameters. 
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Table  4.1,  4.2,  and  4.3  list  all  the  input  information  necessary  for 
the  solution  of  this  problem.  Table  4.1  contains  the  program  parameters 
necessary  for  control  of  the  pattern  search.  Table  4.2  contains  data 
related  to  the  variables,  and  Table  4.3  contains  the  data  as  related  to 
the  objectives  and  achievement  vector. 

The  solution  to  this  problem  is  given  in  Table  4.4  The  achievement 
vector  element  a(l)  =  0.0  indicates  that  the  first  eight  objectives 
(i.e.,  constraints)  have  all  been  satisfied.  The  second  element  a(2) 
measures  the  total  damage  to  enemy  targets  in  terms  of  the  military  value 
placed  on  each  target  set.  The  maximum  attainable  damage  would  be  100 
units.  The  value  a(2)  =  67.55  is  the  under achievement  from  this  maxi¬ 
mum  aspiration  level.  Thus,  the  difference  (100  -  67.55)  =  32.45  units 
is  the  measure  of  total  damage  inflicted  on  all  enemy  targets. 

The  decision  variables  represent  the  best  solutions  for  x  to  attain 
the  values  of  a(l)  and  a(2).  The  reader  is  referred  to  Table  4.5  for 
the  conversion  of  the  solution  to  the  actual  problem  variables. 

The  solution  which  has  been  calculated  has  continuous  values  for  the 
decision  variables.  The  nature  of  the  problem  and  the  variables  demands 
that  they  be  integer  values.  All  non-integer  values  appear  to  be  the 
result  of  round-off  errors  in  the  calculations.  The  solution  at  integer 
values  rounded  to  the  next  closest  value  is: 

a(l)  =  0.0 

a(2)  -  67.71 

This  results  in  a  military  value  for  the  damage  of  32.29,  which  is  slightly 
1'  than  the  continuous  solution  and  differs  by  0.49  percent.  The  computer 
code  gives  an  acceptable,  if  not  exact,  solution  to  the  NLGP  problem. 


Table  4.1  Input  Data  for  Program  Parameters: 
Target  Allocation  Problem 


Number  of  Objectives  = 

9 

Number  of  Priorities  = 

2 

Number  of  Decision  Variables  = 
Maximum  Number  of  Allowable 

24 

Pattern  Search  Cycles  = 

50 

Initial  Step  Sizes  = 

2.0000 

Acceleration  Factor  * 

2.0000 

Step  Size  Reduction  Factor  = 

0.5000 

Minimum  Allowable  Step  Size  = 

0.0001 

VECTOR  OF  INITIAL  ESTIMATES  FOR  DECISION  VARIABLES - : 
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Table  4.3  Input  Data  for  Objectives  and 
Achievement  Vector:  Target 
Allocation  Problem 


Right-hand  side  values  for  objective  functions - RHS 


RHS 

(1) 

27.0000 

RHS 

(2) 

= 

102.0000 

RHS 

(3) 

= 

112300.000 

RHS 

(4) 

= 

147100.000 

RHS 

(5) 

= 

0.0 

RHS 

(6) 

= 

0.0 

RHS 

(7) 

= 

0.0 

RHS 

(8) 

= 

0.0 

RHS 

(9) 

= 

100.0000 

Vector  of  achievement  function  tolerances - EPSY 

EPSY  (1)  -  0.0000100 

EPSY  (2)  =*  0.0000100 


I 

] 

] 


n 


Table  4.4  Target  Allocation  Problem  Solution 


Decision  Variables 


X 

(  1) 

= 

14.999939 

X 

(  2) 

= 

33.999939 

X 

(  3) 

= 

9.999939 

X 

(  4) 

= 

63.999939 

X 

(  5) 

= 

126.999939 

X 

(  6) 

= 

0.999939 

X 

(  7) 

3 

10.999939 

X 

(  8) 

= 

6.999939 

X 

(  9) 

6.999939 

X 

(10) 

576.000000 

X 

(11) 

= 

0.999939 

X 

(12) 

= 

0.999939 

X 

(13) 

= 

0.999939 

X 

(14) 

as 

0.999939 

X 

(15) 

s 

20.000000 

X 

(16) 

=2 

0.999939 

X 

(17) 

= 

76.999939 

X 

(18) 

= 

20.999939 

X 

(19) 

8.999939 

X 

(20) 

* 

0.999939 

X 

(21) 

* 

0.999939 

X 

(22) 

= 

1176.00000 

X 

(23) 

= 

24.999939 

X 

(24) 

= 

8.999939 

Values 

A 

(  1) 

3 

0.0 

A 

(  2) 

3 

67.554916 

i 


Table  4.5  Target  Allocation  Problem  Solution 


Description 


Distribution  of  aircraft 
of  type  j  to  carrier  i 


Missions  flown 
from  carrier  1 
to  target  k 


Missions  flown 
from  carrier  2 
to  target  k 


Using 
aircraft 
of  type  1 


Using 
aircraft 
of  type  2 


Using 
aircraft 
of  type  1 


Using 
aircraft 
of  type  2 


Problem 

Variables 

Solution  for 
Variable 

X11 

14.999 

x12 

33.999 

X21 

39.999 

x22 

63.999 

ylll 

126.999 

y112 

0.999 

y113 

10.999 

y114 

6.999 

y115 

6.999 

y121 

576.000 

y122 

0.999 

y123 

0.999 

y124 

0.999 

y125 

0.999 

y211 

20.000 

y212 

0.999 

y213 

76.999 

y214 

20.999 

?215 

8.999 

y221 

0.999 

y222 

0.999 

y223 

1176.000 

y224 

24.999 

y225 

8.999 

75 


Anytime  a  problem  is  encountered  which  must  have  an  integer  solution 
for  some  of  the  variables,  the  non-integer  values  can  be  rounded  to  integer 
values  and  resolved  for  a  problem  solution.  This  method  cannot  guarantee 
that  a  better  integer  solution  does  not  exist.  This  method  should  result 
in  a  solution  which  is  close  in  value  to  the  best  continuous  solution  to 
the  problem.  It  also  yields  a  good  starting  point  for  the  examination  of 
other  integer  solutions. 

For  comparison.  Rice,  et  al  [28]  obtained  similar  results  for  the 

tota]  damage  inflicted  using  SUMT  [5].  Using  the  target  values  for 

scenario  1,  their  solution  had  33.67  percent  of  the  total  target  value 

damaged,  as  compared  to  32.45  percent  using  the  NLGP  method.  They 

obtained  assignments  of  9  aircraft  of  type  1  and  69  aircraft  of  type  2 

to  carrier  2.  Unfortunately,  this  solution  violates  the  carrier  space 

2 

availability  constraint  for  carrier  2  (148,410  ft  used  versus  147,100 
2 

ft  available).  This  makes  their  solution  infeasible.  The  NLGP  solu¬ 
tion,  although  slightly  smaller  in  value  for  the  damage  percentage,  is 
better  since  it  does  yield  a  feasible  solution.  No  priority  level  one 
objectives  (i.e.,  resouce  constraints)  were  violated.  The  authors 
feel  that  in  more  complex  and  large  problems,  the  superiority  of  the 
NLGP  method  over  single  objective  methods  will  be  more  pronounced. 

4.3  Empirical  Results 

The  relationship  of  problem  size  (i.e.,  the  number  of  objectives  and 
number  of  variables)  versus  computation  time  and  storage  requirements  has 
been  examined.  The  results  of  many  problems  which  were  solved  will  be 
presented  on  a  aggregate  level. 

Testing  of  the  three  computer  codes  proceeded  from  an  initial  phase 
where  problems  were  small  and  simple  in  nature  to  a  secondary  phase  to  the 
final  phase  where  problem  size  exceeded  100  variables  and  100  objectives. 
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In  the  initial  phase,  a  number  of  problems  with  less  than  10  objectives 
and  10  variables  were  solved.  The  second  phase  of  testing  involved  problems 
ranging  in  size  from  10  to  25  objectives  and  25  to  75  variables.  The  last 
phase  involved  problems  with  over  100  objectives  and  over  100  variables.  A 
summary  of  the  results  is  in  Table  4.6. 

The  average  compilation  time  increases  because  the  object  code  increases 
with  problem  size  and  complexity.  This  is  because  the  user  supplies  sub¬ 
routines  OBJFUN  and  ACHFUN  (and  UVALUE,  if  version  3  is  used) .  These  will 
vary  in  size  by  problem. 

The  average  execution  time  of  course  increases  with  increases  in  problem 
size.  It  is  interesting  that  the  difference  between  the  last  two  categories 
is  slight.  This  is  due  in  part  to  the  fact  that  the  problems  with  greater 
than  100  objectives  and  100  variables  had  primarily  linear  objectives 
(usually  100)  and  only  a  few  nonlinear  objectives  (10  to  20).  This  made 
execution  of  these  problems  much  faster  than  would  normally  be  expected  with 
the  majority  being  nonlinear  objectives. 

Computer  storage  versus  problem  size  is  a  difficult  relationship  to 
predict.  The  array  area  is  fairly  constant  as  would  be  expected,  since 
this  area  depends  on  the  dimensioning  within  the  program.  The  object  code 
sizes  varied  and  tended  to  increase  with  increases  in  problem  size.  This 
storage  requirement  will  be  highly  dependent  on  the  type  of  objectives 
fi.e.,  linear  versus  nonlinear),  complexity  of  the  objectives  (i.e., 
quadratic  terms,  exponential  terms,  logarithmic  terms),  the  number  of 
terms  in  each  objective,  and  the  number  of  objectives  in  the  problem. 

These  areas  will  dictate  storage  requirements,  and  also,  execution  time. 


Table  4.6  Computation  Results 
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4.4  Summary 

All  three  codes  have  been  tested  and  give  the  same  solutions  to  the 
same  problems.  Execution  times  of  the  programs  are  fairly  equivalent  for 
all  three  codes.  Version  3  had  no  particular  advantage  over  the  others 
for  the  problems  on  which  it  was  tested,  but  those  problems  were  deemed 
too  small  to  show  appreciable  advantages.  Version  3  should  give  substan¬ 
tial  savings  in  the  area  of  execution  time  over  versions  1  and  2  when  the 
problem  size  is  large  and  the  objectives  are  sparsely  populated  by  the 
variables . 

Version  1  is  the  least  complicated  of  the  three  codes  and  easiest  to 
use.  Version  2  is  the  next  in  line  in  terms  of  ease  of  access.  Version  3 
is  the  most  difficult  to  set  up  for  computation  due  to  subroutine  UVALUE. 
Its  advantages  have  already  been  mentioned. 
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CHAPTER  V 
CONCLUSIONS 


5.1  Conclusions 

The  specific  objectives  of  this  study  have  been  met.  To  meet  the 
demands  for  a  general  computer  tool  to  solve  a  broad  range  of  NLGP  prob¬ 
lems,  three  computer  codes  have  been  developed.  Each  has  advantages 
over  the  others  which  make  it  superior  for  solving  certain  types  of 
problems. 

A  detailed  User's  Manual  has  been  supplied  along  with  this  report 
and  copies  of  the  codes.  The  User's  Manual  is  the  guide  on  how  to  use 
each  of  the  codes.  It  gives  example  problems,  examples  of  the  necessary 
computer  input,  aids  for  selection  of  this  information,  examples  of  the 
output,  and  its  analysis. 

5.2  Recommendations 

The  present  NLGP  algorithm  and  computer  codes  have  limitations  which 
have  been  presented  in  previous  sections.  There  are  certain  areas  in 
which  future  research  efforts  should  be  focused: 

(1)  Examination  and  development  of  other  heuristic  techniques  to 
solve  0-1  nonlinear  gnal  programming  problems. 

(2)  Development  of  methods  to  solve  mixed-integer  NLGP  problems. 

(3)  Development  of  generalized  sensitivity  analyses  for  NLGP  and 

implementation  of  these  methods  into  computer  packages. 

(4)  Analysis  of  these  and  future  codes  using  actual  design  prob¬ 
lems  with  over  1,000  variables  and  over  1,000  objectives. 

This  will  greatly  enhance  future  validation  efforts. 

(5)  Including  and  testing  future  additional  program  modifications 


as  they  are  developed. 
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C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


*** 


*  *  * 


NONLINEAR  GCAL  PROGRAMMING — FATTEFK  El  ARCH  CODF 


C 

c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c - 


THIS  IS  THE  NIGF/MPS-BS  NAVY  CODE 

NONLINEAR  GOAL  PROGRAMMING/ MODIFIED  PATTERN  SEAFC3 
WITH  A  RESOLUTION  RIDGE  SEARCH  TECHNIQUE 
VERSION  1:  SCALAR  STEP  SIZE  (EPS) 

- JULY  15,  1980  - 


♦♦♦CODE  SPECIFICATIONS*** 
2500  OBJECTIVES 
2500  V API AELES 

10  PRIORITY  LEVELS 


CODE 


c 

c 

*  *  *V  A 

♦♦♦NAME*** 

5IAFLES  AND  PARAMETER 
♦♦♦DESCRIPTION*** 

c  *** 

c 

NO  EJ 

NUMBER  OF  CEJECTIV 

ES 

c 

NPKIOR 

NUKBEF  OF  PRIORITY 

IEVELS 

c 

N  VAR 

NUMBER  CF  CECISICN 

VAEIAELES 

c 

NM  AX 

MAXIMUM  NC.  OF  PAT 

TERN  SEARCH  CYCLES  ALLOWED 

c 

NCYCLE 

COUNTER  F CF.  FATTER 

N  SEARCH  CYCLES 

c 

X 

DECISION  VAETAELE 

V ICTCR  USED  FOE  PER1URE ATICN 

c 

XBASE 

BASE  PCINT  VECTCF 

FCE  PATTERN  SEAPCE 

c 

A 

ACHIEVEMENT  VECTOR 

AT  TEST  POINT 

c 

ABEST 

ACHIEVEMENT  VFCTCF 

FCE  PEST  SCIUTICN 

c 

APCS , A  N  EG 

TEMPORARY  AC  HI EV EM 

ENT  VECTOFS  FOR  RIDGE  SEACH 

c 

EPS, IPS 

SCALARS  CF  FERTUFE 

ATICN  STEP  SIZES  USED  FOR 

LBE,  DBD 

OBJ 

RHS 

N,E 

EPSY 

ACCEL 

REDUCE 

DELTA 

IMPROV 


IPFINT 

ITERM 

IPROG 

IACCEL 

SAVE 


■  -C  AFC  1 


EXPLORATORY  MOVES 

LOWER  AND  UPPER  ECIED  VECTORS  FOR  DECISION  VAF 
CEJECT3VE  FUNCTION  VECTOR 

RIGHT-HAND  SIDE  VAIUE  VECTOR  FCR  OEJECTIVE  FNS 

NEGATIVE  AND  POSITIVE  DEVIATION  VARIABLE  VECT. 

TERMINATION  TCIEFANCE  TEST  VECTOR 

ACCELERATION  FACTOR 

STEP  SIZE  REDUCTION  FACTOR 

MINIMUM  ALLOWABLE  STEP  SIZE 

TEST  SWITCH - IF  IMEFCV=1,  X  IMPROVES  THF 

SOLUTION:  IF  IMPEOV=C,  X  DOES  NOT  IMPROVE  THE 
SCIUTICN 

OUTPUT  SWITCH - IF  IFETNT=1,  RESULTS  ARE 

PRINTED  OUT  AT  EVERY  PATTERN  SEARCH  CYCLE; 

IF  I P  FI NT=0 ,  CNIY  TEE  FINAL  REStlLTS  ARE  OUTPUT 

TERMINATION  TEST  SWITCH - IF  ITERM=1,  THE  • 

SCIUTICN  IS  TESTFE  ECH  TERMINATION;  IF  ITEFM=0 
ONLY  A  COMPARISON  TEST  IS  PERFORMED 
COUNTER  EOF  IMPROVEMENTS  DURING  THE  EXPLORA¬ 
TORY  MOVES 

INDICATCF  FCR  ACCEL  IF  ATT  ON  MOVES: 

IACCEL  =  1,  ACCELERATION  MOVE  WAS  SUCCESSFUL 
IACCEL  =-l,  ACCEIEB  ATICN  MOVE  WAS  UNSUCCESSFUL 
IACCEL  =  0,  ACCE IE  F  ATICN  MOVE  NOT  ATTEMPTED 
A  TEMPORARY  SWITCHING  INDICATCR  FOR  PATTERN 
AND  RIDGE  KCVES 


***D  AT  A  INPUT  GUIDE*** 

FORMAT 


nnon  o  non  n  on 
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c 

NCEJ 

15 

c 

NPRIOR 

15 

c 

NVAR 

15 

c 

NH  AX 

15 

c 

IPFINT 

15 

c - 

- CARD  2 - INITIAL  PCINT 

c 

X  ( K)  ,  K=  1  ,  N  V  AR 

8E1C.0 

c - 

- CARD  3 - LOWER  BOUNDS 

c 

LBE(K),  K=  1 ,  SV  AR 

8F1C.C 

c - 

- CARD  4 - UPPER  BOUNDS 

c 

UBE(K),  K=  1 , NV  AR 

eEio.c 

c - 

- CARD  5 - RIGHT-HAND  SIDES 

c 

RHS(I),  1=  1 ,  NOEJ 

8F10.C 

c - 

- CARD  6 - TOLERANCES  FCP  TEFMI 

NATICK 

c 

EPSY  (J)  ,J=1, NPRIOR 

8  F  1 C  .  C 

c - 

- CARD  7 

c 

EPS - INITIAL  STEP  SIZE 

F1C.C 

c 

ACCEL - ACCELSRATIC  N  FACTCF 

F10.0 

c 

REDUCE - STEP  SIZE  REEDCTICN 

c 

FACTOR 

F10.0 

c 

c 

DEITA - STEP  SIZE  TEST  FACTOR 

F1C.C 

c 

c 

*** S  UB RC  DTI NES*** 

c 

SUBROUTINES  OBJFUN  AND  ACHFUN 

MUSI  EE  USEF.  SUPPLIED 

c 

FOR  EACH  NEW  N  LG  P  FRCBLEM 

c 

SUBROUTINE  OBJFUN  CONTAINS  THE  OBJECTIVE  FUNCTIONS 

C  SUBSOU  T IN  E  ACHFUN  CONTAINS  THE  ACHIEVEMENT 

C  FUNCTIONS,  ONE  FOR  EACH  PRIORITY  IEVII 

C 
C 

C  ***£1  AIN  PROGRAM*** 

C 

RE AL*4  IPS, LEE, N 

IN  TEGER*2  IMPROV, IPR I  NT ,ITERM ,IPRCG, I  ACC El 
COSMON/COHMO 1/N03J, NP EIOP , NV A R , NMA X, NCYCIE 
COMMON/CO  MMO  2/X  (2  500)  ,  XEASE  (2500) 

CQMMON/COM MO  3/ A  (10)  ,  ABEST  (10) 

COEMON/COMMO 4/ACCEL,  R EDUCE, DELTA, IMPROV  ,IPEI NT, I TEPH 
CO MMON/COMM05/LBD  (2500)  ,0BD(2500)  , EES, IIS 
CO  EM  ON/CO  M  MO  6/C  BJ  (2500)  ,RHS  (2500) 

COMMON/COMM07/N  (2  500)  ,P{2500) 

CO  MMON/CO  M MO 8/EP  S  Y (10) 

CO  EMON/CO  MMO  9/SAV  E 

* 

BEGIN  THE  FBOGEAM 
CALL  DATAIN 

ALL  NECESSARY  DATA  IS  INPUT  AND  PPINTED  OUT 
CALL  HJALG 

BEGIN  THE  PATTERN  SEARCH  ALGORITHM 
ALL  FURTHER  WORK  IS  COMPLETED  IF  THE  SUERCUTINES 
OP  THE  CODE 
STOP 
END 

***************************************4****4*** **************** ****** 

SUEROUTINE  DATAIN 

*******************************************  44*** 4******* ************** 


SUBROUTINE  DATAIN  READS  IN  ILL  THE  NECESSARY  LATA 
AND  PRINTS  IT  OUT  AS  A  KEAKS  CP  CHECKING  IT. 


c 

c 
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r 


I 


J 

J 


1*1  TEG  IE  *2  IMERCV,  IPP.INT,  IT  ERM,IPFOG,  I  ACC  EL 
R  E  AL  *  4  IP  S,L3D 

CO EMO S/CO MaO  1/S03J,  NPRIOR,NVAE,NMAX,  NCYCIE 
COMMON/COKMO  2/X  (2500)  ,XEASE(25C0) 

CO  MMO  N/COMM04/ ACC  EL  ,  REDUCE  ,DEITA,IEFFCV,IEEINT,ITEPM 
CO  KMOK/COMMO  5/LBD  (2500)  ,UEE(2500)  ,  EPS,IES 
COSMO  S/CO  MMOS/CBJ  (2500)  ,RHS{2500) 

CO  MM  ON/ CO  MMO 8/ EPS Y (10) 

C  BEAD  ALL  THE  NECESSARY  DATA 

READ  170, NCBJ, NPRIOR, NV AR . NK A X , IPE INT 
READ  175,  (X  (K)  ,K=  1,  NVAR) 

HEAD  175,  (LBD  (K)  ,K=1  ,KVA3) 

READ  175,  (UBD(K)  ,  K=  1  ,  NVAR) 

READ  175,  (RHS  (I)  ,1=1  , NCBJ) 

READ  175,  (EPSY  (J)  ,J=  1,NFRICR) 

READ  175, EPS, ACCEL, SEDOCE, DELTA 
C  PRINT  AN  ECHO  CHECK  OF  THE  DATA 

100  PRINT  180, N03J,NPRI0R, NVAR, NMAX,  EPS, ACCEL,  REDUCE, DELTA 
PRINT  181 
DO  11C  K=  1  ,N  VAR 

110  PRINT  182,  K,  X  (K)  ,  K,LED  (K)  ,K,t7ED  (K) 

PRINT  183 
DO  120  1=1, NCBJ 
120  PRINT  184,  I,  RHS  (I) 

PRINT  185 
DO  1 3 C  J=1 ,NPRIOR 
130  PRINT  186 ,J, EPSY  (J) 

PRINT  187 

IF  (I PRINT  .EQ.  1)  GO  TO  140 

PRINT  188 

RETORN 

140  PRINT  189 
170  FORMAT  (1615) 

175  FORMAT  (8F  1  0.  0) 

130  FORMAT  (•-» ,25X,' NUMBER  OF  OBJECTIVES  = »  ,  T57 , 15 ,/ , •  '  ,25X ,' NUMBER  0 
3  F  PRIORITIES  = • , T57 , 1 5 ,/  ,  •  •  ,  2 5 X , ' J 0  EE E E  CE  DECISION  VARIABLES  =  ', 
$T57,I5,/,'  • , 25X , ' MAXIMUM  NOMEEB  OF  ALLOWABLE’ ,/, '  ', 29X, ’PATTERN 
3SEARCB  CYCLES  = ’ , T57 , 1 5,/ ,  •  • , 25X,  '  INITIAL  STEP  SIZES  =',T59,F10. 
*6,/,’  • ,25X, 'ACCELERATION  FACTOR  = •  ,  159  ,  FI  0. 6, /, •  *,25X,’STEP  SIZE 
3 RE DOCTICN  FACTOR  =' , T59, F10. 6,/,'  ', 25X  ,’ MINIM OM  ALIOS ABLE  STEP  SI 
SZE  ='  ,T59,F10.6) 

191  FOFMAT('- ’ ,25X, 'VECTOR  OF  INITIAL  ESTIMATES  FOR  DECISION  VARIAELES 

3 - X’,/,'  ',25X,' LOWER  AND  OPPER  BOUNDS  ON  DECISION  VARIABLES - LB 

3D  AND  OBD') 

182  FORMAT  ('  •  ,  10X  ,  ’  X  (' ,  15,  ' )  =  '  ,  F 1 C.  5, 7  5  ,  •  LBD  (•  ,  15  ,  • )  =  ',F10.5, 

S7X,'UED  (•  ,15,')  =  '  ,F  1 0 . 5) 

183  FORHATC-' ,25X, 'RIGHT-HAND  SIDE  VALDES  ICR  OBJECTIVE  FUNCTIONS - R 

SHS  •) 

184  FORHATC  •  ,25X,'RHS  ('  ,15,' )  =  ’,F20.4) 

185  FORMAT  (’-•  ,25X,' VECTOR  OF  ACHIEVEMENT  FUNCTICS  TOLERANCES - FPSY') 

186  FORMAT  ('  •  ,2  5X  ,  •  E  PS  Y  ( '  ,1 2  ,  ')  =  »,F10-6) 

187  FORMAT ('- ' ,25X,' ALL  DATA  HAS  EEEN  INPUT - FEGIN  THE  PATTERN  SEARCH 

S') 

188  FORHATC  ',25X,’NCTE:  IPBINT  =  0,  SC  CNIY  THE  FIN  A I  SOLUTION  W  ILL 
3BE  PRINTED  OUT') 

189  FORHATC  ’,25X,'NCTE:  IPHINT  =  1,  SC  THE  SOLUTION  WILL  EE  PPINTED 
SOUT  AT  EVERY  PATTERN  SEAFCE  CYCLE') 

RETURN 


non  nn  nn  nnnn  on  nnnn 
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END 

C  **** ******* ******************* ****4********* *****************  ********* 

SUBROUTINE  HJALG 

C  *************  *■***■*****************+***  +  *******+****************  ******* 

c 

C  SUBROUTINE  HJALG  IS  THE  HOCKE-JEEVES  PATTERN 

C  SEARCH  ALGORITHM.  IT  PERFORMS  EXFLOE A10RY  AND 

C  PATTERN  MOVES  TO  FIND  THE  'BEST'  SCIDTICN  TC  THE 

C  NONLINEAR  GCAL  PROGRAMMING  PROBLEM. 

C 

C 

R E  AL  * 4  IP  S  ,  LBD 

INTEGERS  IMERCV,  IPBINT,  IT  EEB  ,  IPROG,  I  ACCEL 
COMMON/COMMO 1/NOB J, NPRICR , RVA B , N MA X , RCY C IE 
CCEMON/COMMO  2/X  (2500)  ,XFASE(25C0) 

COMMON/COM  MO  3/A  (10)  ,ABESI  (10) 

CO  CMC K/CO M MO 4/ ACC  EL, REED C, DELTA, IHFECV  ,  IPRINT, ITEPK 
CO  MMO  N/COMM05/L3D  (250  0)  ,  UBD  (2  5  00)  ,  E  FS,  I  IS 
COMMON/COMM0  9/SAV  E 
CALL  ACU? ON 

SUBROUTINE  ACHEON  FORMS  THE  ACHIEVEMENT  VECTOE  FOE 
THE  INITIAL  BASE  POINT.  IT  IS  TESTEE  ECR  ALL  ZERC 

VAIOES - IF  30,  TEE  INITIAL  POINT  IE  OPTIMAL. 

OTHERWISE,  IT  IS  SET  AS  THE  BEST  VALUE  AND  IS  PFINTEE. 

DC  100  J=1,N PRIOR 
IF  (A  (J)  .  NE.  0.0)  GO  TC  203 
100  CONTINUE 

ALL  PRIORITY  LEVELS  ARE  SATISFIED  AT  ZERC  VALUE. 

THE  INITIAL  POINT  IS  OPTIMAL - TERMINATE. 


200 

PRINT 

290 

DO  201 

J= 

1 ,  nprio: 

201 

PRINT 

282 

,  J,  A  ( J) 

PRINT 

291 

DO  202 

:  K= 

1 ,N  VAR 

202 

PRINT 

280 

,  K,  X  (K) 

GO  TO 

295 

203 

PRINT 

281 

DO  204 

J= 

1, NPRIO 

ABEST  (J) 

*  A(J) 

204 

PRINT 

282 

, J, A (J) 

**  INITIALIZATION  CF  PATTERN  SEAPCE  FARAMETEFS  ** 
NCYCLE  COUNTS  THE  NUMBER  CF  PATTEFK  SFAPCHES  USEE 
I ACCEL  IS  AN  INDICATOR  WHICH  TELLS  IE  THE  EXPICFATORY 
MOVES  ARE  MADE  A3CUT  AN  ACCELEF ATICN  POINT 
NCYCLE  =  0 
IACCEL  =  0 

SET  THE  INITIAL  BASE  PCINT  AND  THE  INITIAL  STEP 
SIZE  BEFORE  THE  EXPLORATORY  MOVES 
DO  205  K= 1 ,N  VAR 

205  XBASE  (!C)  =  X  (K) 

IPS  =  EPS 

START  THE  FATTERN  SEARCE 
**  EXPLORATORY  MOVES  ** 

206  IF  (NCYCLE  .GT.  NH AX)  GO  TC  260 
I  PROG  =  0 

IPROG  CCU NTS  THE  NUMBER  OF  IMPROVEMENTS  DURING  THE 
EXPLORATORY  MOVES.  IT  IS  USED  TO  DECIDE  IF  i  PATTERN 
MOVE  SHOULD  3E  PERFORMED. 

207  DO  210  K=  1  , N V A R 

C  PERTURBATIONS  IS  THE  POSITIVE  DIRECTION 


no  nnnn  nnono  non  nnnn  o  n  nnnon  n  no 
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TEST  STEP  SIZE  COMPONENT.  IF  LESS  THAN  THE  MINIMUM 
ALLOWABLE  STEP  SIZE,  A  FIDGE  SEARCH  WILL  BE  PERFORMED 
IF  (IPS  .LT.  DELTA)  GC  TO  255 
X  { K)  =  X(K)  +  IPS 
I?  (X  ( K)  .LT.  LBD  (K)  )  GC  TO  208 
I?  (X  (  K)  .  GT.  UED  (K)  )  GO  TO  206 
CALL  ACBFUN 

ACHFON  FORMS  THE  ACHIEVEMENT  VECTOR  AT  THE  TEST  POINT 
ITESM  =  0 
CALL  DECIDE 

DECIDE  COMPARES  TEE  ACHIEVEMENT  VECTCES  OF  THE 
BASE  POINT  AND  THE  TEST  POINT.  IF  X  (K)  IMPROVES  THE 
SOLUTION,  IT  BECOMES  A  TEMPORARY  EEA I  PCINT  AND 
EXPLORATORY  MOVES  ARE  CONTINUED  FRCM  IT.  IF  THFPF 
IS  NC  IMPROVEMENT,  THEN  EXPLORE  IN  TEE  NEGATIVE  DIRECTION 
I?  (IM  PROV  .EQ.  0)  GO  TO  208 
I  PROG  =  I  PROG  +  1 
GO  TO  210 

PERTURBATIONS  IN  THE  NEGATIVE  DIRECTION 

208  X  ( K)  =  X (  K)  -  2.  0 *1  PS 

I ?  (X  (  K)  .  LT.  LED  (K)  )  GO  TO  20S 
I?  (X  (K)  .GT.  UBD  (K)  )  GC  TC  20S 
CALL  ACHFUN 

ACHFUN  FORMS  THE  ACHIEVEMENT  VECTOR  AT  THE  TEST  PCINT 
ITERM  =  0 
CALL  DECIDE 

IF  X  (K )  IMPROVES  THE  SCIUTICN,  IT  BECOMES  THE 
TEMFCRAFY  HEAD  POINT  AND  THE  NEXT  VAFIAELE  IS 
PEFTUREED.  IF  THERE  IS  NO  IMPROVEMENT,  XBASE(K) 

REMAINS  THE  TEMPORARY  HEAD  POINT. 

IF  (I M  EROV  .  EC-  0)  GO  TC  209 
IP HOG  =  IPROG  ♦  1 
GC  TO  210 

209  X  (K)  =  X  ( K )  +  IPS 

210  CONTINUE 

NCYCLE  =  NCYCLE  *  1 

THE  EXPLORATORY  MOVE  DC  LCCE  HAS  CONSTRUCTED  A 
NEK  POINT.  TEST  TEE  ACHIEVEMENT  VECTOR 
FOR  ALL  ZERO  VALUES 
DO  211  J= 1 , NPRIOR 
I?  (A  ( J)  .NE.  0.000)  GC  TO  212 

211  CONTINUE 
GO  TO  200 

THE  ACHIEVEMENT  VECTOR  IS  NOT  ALL  ZEROES 
TEST  THE  POINT  FOR  IMPROVEMENT.  IF  IT  GIVES 
IMPROVEMENT,  A  PATTERN  MOVE  IS  PERFORMED.  IF  NOT 
THEN  STEP  SIZES  Will  BE  REDUCED  AND  JEW  EXPLOF ATCFY 
MOVES  WILL  EE  PEETCRMED. 

212  IP  (IPROG  .GT.  0)  GO  TC  214 
IF  (I ACCEL)  250,213,215 

THE  EXPLORATORY  MOVE  HAS  EEEN  UNSUCCESSFUL.  STEP 
SIZES  WILL  EE  REDUCED  AND  NEW  EXFLCF  ATOPY  MOVES 
WILL  BEGIN. 

213  PRINT  283 

*<■  REDUCE  STEP  SIZES  ** 

IPS  =  IPS  *  REDUCE 
GO  TO  206 

THE  EXPLORATORY  MOVE  HAS  EEEN  SUCCESSFUL. 

A  PATTERN  MOVE  KILL  BE  PEPFCRMIC 

214  PRINT  279 


nnnnn  no  non  no  non 
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C  **  ACCELERATION  ** 

215  DO  220  K=  1  ,  N  VAR 
SAVE  =  X (  K) 

X  { K)  =  X  ( K)  ♦  ACCEL  *  (X(K)  -  XEASI(K)) 

IF  {X  { K)  .LT.LEE(K)  )  X  (K)  =  IBD(K) 

IF  (X  (  F)  .LE.  UED  (K)  )  GO  TO  220 
X  (K)  =  tJB  D  (K) 

220  XBASE  (K)  =  SAVE 

C  TEST  THE  ACHIEVEMENT  VECTOR  AT  THE  ACCELERATION  POINT 

CAIL  ACHFUN 
ITERM  =  1 
CALL  DECIDE 

IF  (IM  P30V  .EQ.  1)  GO  TC  225 
I ACC EL  =  -1 
GO  TO  230 
225  I ACC EL  =  1 
GO  TO  240 

THE  ACCELERATION  PT.  GIVES  NO  I ! P  FC V  E  M  ENT. 

EXPLORATORY  MOVES  WILL  EE  P  ER  FORMED  ABOUT  THE 
ACCELERATION  PT.  TC  FIND  I FFRCVEMENT. 

230  PRINT  286 

231  IF  (IP SINT  .EQ.O)  GO  TO  206 
PRINT  297 , NCYCLE 
DO  23  2  K= 1 ,N VAR 

232  PRINT  280  ,  K,  XEASE  (X) 

PRINT  238 
DC  235  J=1,NPEIOR 

235  PRINT  282, J, ABEST  (J) 

GC  TO  206 

THE  ACCELERATION  POINT  GIVES  IMEFCVE  KENT. 

EXPLORE  ABOUT  IT  EOF  IMPROVEMENT. 

240  PRINT  285 

IF  (IP PINT  .EQ.  0)  GO  TC  206 
PRINT  287,  NCYCLE 
DO  24  1  K=  1  ,  NV AF 

241  PRINT  230,  K,  X  (K) 

PRINT  288 
DO  245  J= 1 ,NPSIOS 

245  PRINT  282  ,J, A  (J) 

GO  TO  206 

EXPLORATION  AROUND  THE  ACCEIE  FATI CN  ICINT 
DOES  NOT  YIELD  ANY  IMPROVEMENT 

RESET  THE  BASS  POINT  AND  GC  TC  EX  PLC  £  ATCRY  MOVFS. 

250  IF  (IACCEL  .EQ.  0)  GO  TC  255 

251  DO  252  K=  1  ,N  VAR 

252  X  (  K)  =  XBASE  (K) 

IACCEL  =  0 
GO  TO  206 

STEP  SIZE  REDUCTIONS  EQUALLED. 

A  RESOLUTION  RIDGE  SEARCH  WILL  EE  FEFECEBED 
255  PRINT  284 
259  CALL  RIDGE 

SUBROUTINE  RIDGE  WILL  TRY  TO  FIND  A  NEW  BASE  POINT 
IN  AN  OBLIQUE  DIRECTION  TC  THE  PRESENT  EXPLCR  ATCRY 
SEARCH  DIRECTION.  IF  SUCCESSFUL,  THE  PATTERN 
SEARCH  WILL  BEGIN  AT  THIS  NEW  PCI  NT .  IT  UNSUCCESSFUL, 
THE  ALGORITHM  WILL  EE  TERMINATED. 

IACCEL  =  0 
GO  TO  206 

C  THE  MAXIMUM  NUMBER  OF  PATTERN  SEAFCH  CYCLES  HAS  BFEN 


non 
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EXCEEDED.  THE  BEST  VALUES  OF  TEE  VAFIAPLES  AND  THE 
ACHIEVEMENT  VECTOR  AEE  ERINTED  CUT,  AND  THE  A  LGCRITH  E 
WILL  BE  TERMINATED. 

260  PRINT  239 

IF  {IACCEL  .LT.1)  GO  TO  265 
DO  261  K=  1  ,  N  V  A  B 

261  PRINT  280 ,  K,  X  { K) 

PRINT  298 
DO  262  J= 1 , NP3I0R 

262  PRINT  282,0,  A  (J) 

GO  TO  295 

265  DO  266  K=  1  ,NVAR 

266  PRINT  280,  K,  XBASE  (K) 

PRINT  288 

DO  267  J=1,N PRIOR 

267  PRINT  282  ,  J,  AEEST  (J) 

270  GO  TO  295 

279  FORMAT ('-  * ,25X,'THE  EXPLORATORY  MOVE  HAS  BEEN  SUCCESSFUL ',/, •  *,25 
SX, 'A  PATTERN  MOVE  KILL  BE  FEFFCBCED ') 

230  FORMAT  ( *  '  ,  3  4X  ,  •  X  (’  ,  15  ,  •  }  =  *,F13.6) 

231  FORMAT  ('  •  ,2  5X  , 'I  NITI A L  ACHIEVEMENT  VECTOR  VALUES') 

282  FORM  AT  ( '  ' ,34X, ' A  {',12,')  =  *,713.6) 

283  FORMAT ('-* ,25X, 'TEMPORARY  HEAD  POINT  DCES  NOT  IMPFCVE  THE  SOLUTION 

S',/,'  ' ,2  5  X,  ' THE  STEP  SI2ES  WIIL  EE  RFCCCED' ,/, *  ',25X,'NEW  EXFLOP 

SAT  CRY  MOVES  WILL  EE  PERFORMED  ABOUT  X  F  A  S  E ' ) 

234  FORMAT (»  ' ,2 5X ,' MAXIMUM  STEP  SIZE  REDUCTIONS  EXCEEDED  FCR  THIS  EXP 
SLORATION*  ,/, '  ' , 2  5X , *  A  RIDGE  SEARCH  Kill  EE  EEEFOFMEE* ) 

285  F0EMAT('-'  ,25X,'THE  ACCELERATION  POINT  YIELDS  AN  IMPROVED  SOIUTICN 
3',/,'  *,25X, 'EXPLCRATCRY  MCVES  Will  EEG3N  AT  THIS  NEW  PCINT') 

286  FOFMAT ( ' ,25X,'THE  ACCELERATION  PCINT  CCES  NOT  YIELD  AN  IMRCVED 

jSCLUTION* ,/, '  ',25X, 'EXPLCRATCRY  MCVES  KILL  BEGIN  AT  THIS  POINT') 

287  FORMATC- '  ,25X, 'THE  BEST  SCLOTICK  AT  PATTERN  SEARCH  NUMEER  ',13,' 
5FGILOWS ') 

288  FORMATC  ',25X,'THE  ACHIEVEMENT  VECTOR  VALUES  ARE:*) 

289  FOFMATC  '  ,25X,'  MAXIMUM  NUMEEF  OF  PATTERN  SEARCH  CYCLES  EXCEEDED', 
$/,'  '  ,25X , 'THE  ALGORITHM  IS  TERMINATED',/,'  »,25X,'THE  BEST  SCIUTI 
SON  TO  THIS  POINT  FOLLOWS') 

290  FORMAT  ('-' ,20X, 43  (1H*)  ,/,  *  • , 2 5X , • OP TIM AL  ACHIEVEMENT  VECTOR  VALUE 

SS',/,'  * ,  20X  ,  4  3  (1 H*)  ) 

291  FORMATC-' ,20X, 43  (1H*)  ,/,*  *  ,25X  ,  •  CETI  E  A I  DECISION  VARIAELE  VALUES 
S',/,'  '  ,20X,  43  (1H*)  ) 

295  STOP 
END 


SUBROUTINE  DEWAR 

C  ********************************************************* 

c 

C  SUBROUTINE  DEWAR  CALCULATES  THE  POSITIVE  AND 

C  NEGATIVE  DEVIATION  VARIABLES  FCR  EACH  CEJECTIVE 

C  FUNCTION  USING  THE  PREVIOUSLY  EVALUATED  OBJECTIVE 

C  FUNCTIONS  AND  THE  RIGHT-HAND  SICE  VAIUES. 

C 

c 

REAL* 4  N 

CO  MMON/COMMO  1/N03J,  NPRIO R  ,  KV A  R  ,  N CA X  ,  KC  Y CL E 
COCMO K/COSMO  2/X  (2500)  ,  XE AS  E  ( 2500) 

COHflON/COMH06/OBJ  (2500)  ,RHS  (2500) 

COCaON/COSM0  7/N  (2500)  ,P(2500) 

CALL  CBJFUN 
400  DO  430  1=  1  ,NOSJ 
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IF  (OBJ  (I)  -  2HS(I)) 410,420,420 
410  M  ( I)  =  RHS  (I)  -  OEJ  (I) 

P (I)  =  0. 0 
GO  TO  430 
420  N  (I)  =  0.  00 

P  (I)  =  OBJ  (I)  -RHS  (I) 

430  CONTINUE 
RETURN 
END 

C  *********** *********************************************************** 

SUBROUTINE  DECIDE 

C  ********************************************************************** 
C 

C  SOBHOOTINE  DECIDE  IS  USED  FOR  TVC  FOIFCSFS: 

C  (1)  TO  COMPARE  TEE  ACHIEVEMENT  VECTOR  VALUES  OF 

C  THE  BASE  POINT  AND  A  TEST  POINT  DURING  THE  EXPLORATORY 

C  MOVES  AND  SELECT  THE  BEST  ECINT  FCF  FUTURE  MOVES 

C  (2)  TO  TEST  AN  ACHIEVEMENT  VEC10F  FOP  TERMINATION 

C  OF  THE  ALGORITHM.  THIS  CAN  CCCOF  IN  T  SC  METHODS: 

C  *  *  **  TERMINATION  EY  MEANS  OF  ZERC  FOR  ALL 

C  PRIORITY  LEVELS,  AND 

C  *  *  **>£  EHM I  NATION  EY  MEETING  SPECIFIED  TOLERANCES 

C 

C  THE  PARAMETER  ITERM  IS  AN  INDICATOR  FCF  TERMINATION 

C  TESTING.  I?  ITERK  =  0,  CNIY  THE  COMPARISON  TEST  IS 

C  PEFFCRMFD.  IF  ITE3M  =  1,  THE  TERMINATION  TESTS  WILL 

C  ALSO  BE  PERFORMED. 

C 

C 

INTEGERS  IKPRCV,  IP3INT,  ITERM 
COMMON/COMMO 1/NOB J, NPRIOR , NVA E , S EA I,  NCYCLE 
CC  SMC  S/CO  M  MO  2/X  (2500)  , XBASE  (2500) 

COMMON/COMMO  3/A  (10)  ,  ABEST(  10) 

CO MMON/CO  MM04/ ACCEL ,REDUC ,BE1T A,I EFECV ,IFFINT,ITERM 
COMMON/COM  MO  8/E  PS  Y  (10) 

C  **  COMPARISON  OF  ACHIEVEMENT  VECTORS  FOE  EXPLORATORY  MOVES  ** 

600  DC  610  J= 1 , NP RIOH 

I?  (A  (J)  .GT.  ABSST(J)  +0.0001  )  GO  TO  63C 
IF  (A  ( J)  +0.0001  .LT.  A  BEST  (J)  )  GO  TC  62C 
IF  (J  .EQ.  NPRIOR)  GO  TO  630 
610  CONTINUE  ' 

C  THE  TEST  POINT  IMPROVES  THF  SOLUTION 

620  I M FRO V  =  1 
GO  TO  640 

C  THE  TEST  POINT  DOES  NOT  IMPROVE  THE  SOLUTION 

630  13  PRO  V  =  0 

C  CHECK  THE  INDICATOR  FOR  FURTHER  TERMINATION  TESTING 

640  IF  (ITERM  .  EQ  .  0)  G  0  TC  670 
CONTINUE 

C  **  TERMINATION  TESTS  ** 

C  (1 )TERMINATION  TEST  EY  MEANS  OF  ZERC  FOB  ALL  PRIORITIES 

650  DO  655  J= 1, NPRIOR 

IP  (A  ( J)  .  NE.  0.0)  GO  TC  660 
655  CONTINUE 

C  ALL  PRIORITY  LEVELS  ARE  ZERC THE  OPTIMA!  SOLUTION 

C  HAS  SEEN  FOUND - TERMINATE  THE  ALGORITHM 

GO  TO  680 

C  (2 ) TER  MIN  ATION  TEST  BY  TOLERANCES 

660  DO  665  J= 1 , NPRIOR 

If  (ABS  ( A  ( J)  -  AE2ST  (J)  )  .GT.  EPSY  (J)  )  GC  TO  670 
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665  CONTINUE 

C  All  TOLEF  ANC2S  ARE  SATISFIED - A  SCIUTICN  HAS  BEEN 

C  FOUND - THE  SEARCH  ALGORITHM  IS  CCME1ETE 

GO  TO  680 

C  A  TOLERANCE  IS  NOT  SATISFIED - A  EETLEE  SOLUTION 

C  HILL  BE  SOUGHT - RETURN  TC  THE  RATTEEN  SEARCH 

670  I?  (IM  EROV  .  EQ.  0)GO  TC  695 
DO  675  J=1 ,  NPRIOR 
675  ABEST  (J)  =  A  (J) 

GO  TO  695 
680  FRINT  690 
PRINT  691 
DC  682  J=  1  , NPRIOR 
632  PRIST  692,  J,  A  (J) 

PRINT  693 
DO  68  5  K=  1 ,  N  VAR 
635  PRINT  694  ,  R,  X  (K) 

STOP 

690  FORMAT  (•-' ,25X, *T HE  SEARCH  ALGORITHM  IS  COMPLETE - ALL  TOLERANCES 

SFOR  A  SOLUTION  ARE  SATISFIED*) 

691  FORMAT  (*-• ,20X, 43  (1 H*)  ,/, *  *  ,25X, *  OPTIMAL  ACHIEVEMENT  VECTOR  VALUE 
5  S' ,/, *  *  ,  20X  ,43  ( 1H*) ) 

692  EO  RM  AT  ( *  •  ,  34  X  ,  •  A  (•  ,  I  2  ,  •  )  =  *,F13.6) 

693  FORMAT  (*-' ,20X, 43  < 1 U*)  ,/, *  25X, 'OPTIMAL  DECISION  VARIABLE  VAIUES 

3*,/,'  *,20X,43  <1H*)  ) 

694  FOFMAT  (*  •  ,3  4X  ,  *  X  (*  ,  15,  '  )  =  *,F13.6) 

695  RETURN 
END 

C  **********************************************  ************************ 

SUBROUTINE  RIDGE 

C  ********************************************************************** 
C 

C  SUEROUTINE  PIDGE  IS  CALLED  WHEN  TEE  PATTERN  SEARCH 

C  CAN  NO  LONGER  FIND  PATTERN  CR  EXFLCF  ATOPY  MOVES 

C  SHICH  IMPROVE  THE  ACHIEVEMENT  VECTCE.  RIDGE 

C  EVALUATES  EXPLCP  ATCEY  POINTS  HB2CH  ARE  IN  OBLIQUE 

C  DIRECTIONS  TO  THE  USUAL  EXRICFATI CN  AXES.  RIDGE 

C  ATTEMPTS  TC  FIND  A  RESOLUTION  RIDGE,  IF  IT  EXISTS, 

C  AND  MOVE  THE  PATTERN  IS  THAT  DIFECTICN. 

C 

C 

DIMENSION  APCS  (10)  ,  ANEG{  10) 

REAL*  4  IP  S , L3D 

INTEGER*2  IMPRCV,  IPR INT, IT ERK 

CO  MMO  N/CO  3M0 1/NOB  J, NPRIOR, NVAR,NCAX, >CY  CIE 

CO  MMO  K/COM  MO  2/X  (2  50 0 )  ,  X E  AS  E  ( 2  5 00) 

CO MMO N/CO MMO 3/ A (10)  , ABEST  (10) 

C0M3ON/COMH0  4/ACCEL, R E CUC E , D EITA , IM FRO  V , IPRINT ,1  TERM 
CO MMO N/CO MMO 5/ LBD  (2500) ,UBD (2500)  ,  EES, I ES 
CO  KMC  N/CO  MMO  9/SAVE 

C  PERTURB  EACH  VARIABLE  CNE  AT  A  TIEE 

C  PERTURB  IN  THE  POSITIVE  DIFECTICN 

DO  800  K=  1  , N V A R 

800  X  { K)  =  XBASE  (K) 

801  DO  865  K*  1 , N V A R 

X  (K)  *  X  (  K)  ♦  IPS 

I?  (X  (  F)  .  LT.  LED  (K)  )  GO  TC  805 
IF  (X  ( K)  .GT.  UBD  (K)  )  GC  TC  810 
GO  TO  815 
805  X  (K)  =  LBD  (K) 


nnnnn  n  n  n  n  no  non  nn  o  nn 
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GC  TO  815 
810  X  (K)  =  OSD  (K) 

8 15  CALL  ACHFUN 

ACHFUN  FORMS  THE  ACHIEVEMENT  VECTOR  IN  THE  POSITIVE 
DIRECTION 
DO  820  J= 1 ,N?RI0R 
820  A?OS(J)  =  A  { J) 

PERTORB  IN  THE  NEGATIVE  DIRECTION 
X(K)  =  X { K)  -  2.0  *  IPS 
IF  (X  {  K)  .  LT.  LED  (K)  )  GC  TC  825 
I?  (X  (K)  .GT.  ODD  (K)  )  GC  TC  830 
GO  TO  935 
825  X  (K)  =  LBD  (X) 

GO  TO  835 
830  X  { K>  =  UBD  IK) 

835  CALL  ACHFUN 

ACHFUN  FORMS  THE  ACHIEVEMENT  VECTOR  IN  THE  NEGATIVE 
DIRECTION 
DO  840  J= 1 ,  NP5IGR 
840  AN  EG  ( JJ  =  A  { J) 

COMPARE  APCS(J)  ARE  ANEG(J).  SELECT  THE  MINIMUM 
VIA  A  STEPWISE  PREEMPTIVE  PROCEDURE.  THIS  DETERMINES 
A  DIRECTION  FOR  THE  RIDGE  POINT. 

DO  850  J=  1  ,  N PRIOR 

IF  (APCS  (J)  .LT.  A  NEG  ( J) )  GC  TO  855 
IF  ( APO  S  (J )  .GT.  A  NEG  ( J)  )  G  C  TC  360 
I?(J  .EC.  NP RICH)  GO  TO  855 
850  CONTINUE 

855  X  ( K)  =  XBASE  (K)  ♦  IPS 
GO  TO  861 

860  X  ( K)  =  XBASE  (K)  -  IPS 

RETURN  X(X)  TO  ITS  BASE  PCINT  VALUE,  XEASE,  3EFCRE 
THE  PERTURBATION  OF  THE  NEXT  VARIAEII 

861  SAVE  =  XBASE  (K) 

XBASE  (K)  =  X  [K) 

365  X  ( K)  =  SAVE 

THE  NEW  POINT,  XRIEGE,  HAS  BEEN  FORMED  IN  AN  OFIICUE 
DIRECTION  TO  THE  OLD  EXELCEATCEY  KCVES.  ITS 
ACHIEVEMENT  VECTOR  IS  COMPUTED  AND  CCMPAF.ED  TO  THE 
BEST  VALUE  AT  POINT  XBASE. 

DO  866  K=  1  , N V A R 
SAVE  =  XBASE  (K) 

XBASE  (K)  =  X  (K) 

866  X  (K)  =  SAVE 

CAIL  ACHFUN 
I  TERM  =  0 
CAIL  DECIDE 

IF  THE  RIDGE  POINT  GIVES  NC  1 3FFCVES  ENT ,  THE  PRESENT 
BASE  POINT  IS  THE  BEST  SOLUTION,  A.KE  THE  ALGORITHM 
IS  TERMINATED. 

IF  THE  RIDGE  POINT  GIVES  IMEBOVEMEKT,  AN  EVEN  EETTER 
POINT  WILL  EE  SOUGHT  ALCNG  THE  FESCLCTICN  RIDGE. 

I?  (IN PROV  .EQ.1)  GO  TO  869 
CONTINUE 

C  XBASE  IS  THE  BEST  SOLUTION 

PRINT  891 
PRINT  897 
DO  867  K=1,NVAR 

867  PRINT  892, K, XBASE  (K) 

PRINT  898 


UNCLASSIFIED 

2  nr  2 
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DO  86 8  J=1,NPRI0S 

868  PRINT  893  /  J,  AEEST  (J) 

STOP 

THE  RIDGE  POINT  GIVES  AN  IXFROVED  SOIOTION. 

TRY  TO  ACCELERATE  ALONG  THE  RESOLUTION  RIDGE 

869  PRINT  894 

870  DO  87  1  K=  1  ,NVAR 
SAVE  =  X{  K) 

X  ( K)  =  X(K)  +  ACCEL  *  (X  (K)  -  XBASE  (K) ) 

871  XBASE  (K)  =  SAVE 
CALL  ACHED N 
ITERM  =  0 

CALL  DECIDE 

I?  THE  ACCELERATION  POINT  YIELDS  AN  IMPROVED  SOLUTION, 

TRY  TO  ACCELERATE  FARTHER  ALONG  THE  FESCLUTION  PIEGE 
IF  THERE  IS  NO  IMPROVEMENT,  TEE  FIDGE  POINT  BECOMES 
THE  NEH  POINT  FOR  THE  RENEKED  EATTEFN  SEARCH. 

IF  (I  B  FROV  .Eg.  0)  GO  TO  872 
PRINT  896 
GO  TO  870 
8  72  PRINT  395 
PRINT  897 
DO  873  K=  1 ,N  VAR 
873  PRINT  892  ,  K,  XESSE  {K ) 

PRINT  898 
DO  875  J=  1  ,NPEIOR 
875  PRINT  893,  J,  A3EST  (J) 

RESET  THE  EASE  POINT  AND  STEP  SIZE  VECTORS  AND 
RETURN  TO  THE  EX PICRATCRY  MOVES . 

NO  IMPROVEMENT  FROM  ACCELERATION  POINT 
DO  877  K=  1  ,  N  V  A  B 
877  X  (K)  =  XBASE  (K) 

IPS  =  EPS 
DO  879  J=1,N PRIOR 
879  A  (J)  =  ABEST(J) 

GO  TO  899 

891  FORMAT ,25X,'THE  RIDGE  SEARCH  IS  COMPLETE',/,'  *,25X,'NO  IMPRCV 
SED  SOLUTION  CAN  BE  FOUND',/,*  ',25X,'THE  CURRENT  EASE  POINT  IS  THE 
S  BEST  SOLUTION  THAT  CAN  BE  FOUND  THROUGH  THE  PATTEFN  SEARCH*) 

892  FORMAT  ('  '  ,3  4X,'X  (»,  15,  *)  =  ',F13.6) 

893  FORMAI(*  »  ,34X,  *  A  (»  ,12  ,  * )  =  *,F13.6) 

894  FORMAT  {'-• ,25X,'TBE  FIDGE  POINT  GIVES  AN  IMPROVED  SOLUTION* ,/,  '  *, 
S 25X, ' TRY  TO  ACCELERATE  ALONG  THE  RESOLUTION  RIDGE*) 

895  FOFMATC-' ,25X,' ACCELERATION  ALONG  THE  FESCLUTION  RIDGE  YIELDS  NO 
5IME30VED  SOLUTION',/,'  *,25X,'THE  RIDGE  POINT  IS  THE  NEK  BASE  PCIN 
ST') 

896  FOFMATC-*  ,25X,'  ACCELERATION  ALONG  THE  EESOLUTION  RIDGE  YIELDS  AN 
SIMPROVED  SOLUTION',/,'  ',25X,'TBY  TO  ACCELERATE  FURTHER ' ) 

897  FORMAT  (•-' ,25X, 'DECISION  VARI  AELES  * ) 

898  FORMAT  ('-•  ,25X,»  ACHIEVEMENT  VALUES') 

899  RETURN 
END 

*********************  *****************  4***4  4 ******************** ****** 

SUEROUTINE  ACHFUN 

***************************************  **** 4  ***************  *********** 


REAL*  4  N 

COMaON/COHHO 1/NOB J, NPBIOB, NV1E, NHAX, NCYCLE 
CO  MM  ON/CO  MMO  2/X  (2  500)  , XBASE  (2500) 
COMHON/COMM0  3/A  (10)  ,  AB  ESI  (10) 


n  n 
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COBMOR/COHK07/N(2SOO)  #P  (25C0) 

CA1L  ESVVAB 

C  ACHIEVEMENT  FUNCTIONS  ARE  ENTERED  HEFE. 

RETURN 

END 

C  ************* ********************************************************* 

SUBROUTINE  OBJFUN 

********************************************************************** 

COEMON/COHBO 1/NOBJ, NFBIOS, NVAR,NMAX, KCYCIE 
CO MSON/CO M MO 2/X  (2500)  , XBASE  (2500) 

CO KHON/CO3H06/OBJ (2500)  ,RHS  (2500) 

C  OBJECTIVE  FUNCTIONS  ARE  ENTERED  HEBE 

RETURN 
END 


APPENDIX  B 


NLGP/MPS-RS  CODE:  VERSION  2 
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C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


**************** *********  *  ************************************* ******* 
NONLINEAR  GOAI  PROGRAMMING —  PA  HERN  SEARCH  CODE 
********************************************************************** 


THIS  IS  THE  NLGE/MPS-HS  NAVY  CCEE 

NCNLINEAB  GCAL  EROGR AMMING/MODIEIED  PATTERN  SEARCH  CODE 
WITH  A  BESCLOTICN  PICGE  SEARCH  TECHNIQUE 
VEFSICN  2:  VECTCB  STEP  SIZES  (EPS) 

-  JULY  IE,  1980  - 


♦♦♦CODE  SPECIFICATIONS*** 
2500  OBJECTIVES 
2500  VARIABLES 

10  FRICRITY  LEVELS 


***VAEI AELES  AND  PARAMETERS*** 
***NAME***  ***DE'CEIPTION*** 


NCBJ 

NEEIOF 

NVAB 

KMAX 

NCYCLE 

X 

XEASE 

A 

A  BEST 
APOS, ANEG 
EES, IPS 

LEE, OBD 

CBJ 

PHS 

N.P 

EFSY 

ACCEL 

PEDOCE 

DELTA 

IMPROV 


IPRINT 


ITERH 


IFBOG 
I ACC EL 


SAVE 


NOREEF  OP  OBJECTIVES 
NOEEER  CF  PRIORITY  LEVELS 
NOEEER  CF  DECISION  VABIAELES 

MAXIMUM  NO.  CE  PATTERN  SEARCH  CYCLES  ALLOWED 
CCCNTEF  FOR  PATTERN  SEARCH  CYCLES 
DECISION  VARIABLE  VECTOR  USED  FOR  PERTURBATION 
EASE  POINT  VECTOR  FCR  PATTERN  SEARCH 
ACHIEVEMENT  VECTOR  AT  TEST  POINT 
ACHIEVEMENT  VECTOR  FCR  BEST  SOLUTION 
TEE ECRABY  ACHIEVEMENT  VECTORS  FOR  RIDGE  SEACH 
PEETOPBATION  STEP  SIZE  VECTORS  USED  FOR  EC* 
EXELCRATORY  MCVES 

LCEEB  AND  UPPER  BOUND  VECTORS  FOR  DECISION  VAF. 
CEOECTIVE  FUNCTION  VECTOR 

RIGHT-HAND  SIDE  VAIUE  VECTOR  FOR  OBJECTIVE  FNS 

NEGATIVE  AND  POSITIVE  DEVIATION  VARIABLE  VECT. 

TEFNINATION  TOLERANCE  TEST  VECTOR 

ACCELERATION  FACTOR 

STEP  SIZE  REDUCTION  FACTOR 

MINIMUM  ALLOWABLE  STFP  SIZE 

TEST  SWITCH - IF  IMPBOV=1,  X  IMFFOVES  THE 

SOLUTION;  IF  IMPBOV=0,  X  DOES  NOT  IMPROVE  THE 
SOLUTION 

COTEUT  SWITCH - IE  IFRINT=1,  RESULTS  ARE 

PRINTED  OUT  AT  EVERY  PATTERN  SEARCH  CYCLE; 

IF  IPFI NT=0 ,  CNLY  THE  FINAL  RESULTS  ARE  OUTPUT 

TEFMINATION  TEST  SWITCH - IF  ITERM* 1  ,  THE 

SC1UTION  IS  TESTED  FCR  TERMINATION;  IF  ITEFM=0 
ONLY  A  COMPARISON  TEST  IS  PERFORMED 
CCUNIER  FOP  IMPROVEMENTS  DURING  THE  EXPLOBA- 
TCET  MOVES 

INDICATOR  FOR  ACCELERATION  MOVES; 

IACCEL  =  1,  ACCELERATION  MOVE  WAS  SUCCESSFUL 
I  ACCEL  =-1,  ACCELERATION  MOVE  WAS  UNSUCCESSFUL 
IACCEL  =  0,  ACCELERATION  MOVE  NCT  ATTEMPTED 
A  TEMECPARY  SWITCHING  INDICATOR  FOR  PATTERN 
AND  RIDGE  MOVES 


•••DATA  INPUT  GUIDE*** 


non  n  on 
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C 

C - CARD  1  FORMAT 

C  NCBJ  15 

C  .  NERICB  15 

C  NVAR  15 

C  NMAX  15 

C  IPRI NT  15 

c - CARD  2 - INITIAL  POINT 

C  X  (F)  ,  F=1#NVAF  8 F 10. 0 

C - CARD  3 - EERTUREATICN  STEP  SIZES 

C  EPS  (K)  ,K=1,KVAR  8710.0 

C - CARD  4 - LOWER  ECONLS 

C  LBD(K),  K=  1  #  NT AB  8F10.0 

C - CARD  5 - UPPER  BOONES 

C  UBD(K),  F=1,NVAR  8710.0 

C - CARD  6 - EIGHT-HANE  SIDES 

C  RHS(I),  I=1,NCEJ  8710.0 

C - CARD  7 - TOLERANCES  FOP  TERMINATION 

C  FFSY  (J) ,J=1,NPRICR  8710.0 

C - CARD  8 

C  ACCEL - ACCELEEATICK  FACTOR  710.0 

C  REDUCE - STEP  SIZE  FEDUCTION 

C  FACTOR  F 10. 0 

C  DELTA - STEP  SIZE  TEST  FACTOR  F10.0 

C 

c 

c  ***SGEFCUTINFS*** 

C  SUBROUTINES  OPJFUN  AND  ACHFUN  MUST  BE  USER  SUPPLIED 

C  FCR  EACH  NEW  NIGF  PFCELEH 

C  SUBROUTINE  OBJEUN  CONTAINS  THE  OBJECTIVE  FUNCTIONS 

C  SUBROUTINE  ACHFUN  CCKTAINS  THE  ACHIEVEMENT 

C  FUNCTIONS,  ONE  FCR  EACH  PRIORITY  IEVEI 

C 
C 

C  ***H AIN  EFCGRAfl*** 

C 

REAL*«  IFS,LBD, N 

INTEGER*2  Ifl EROV, I PRINT ,ITERM ,IPFCG ,1 ACCEL 
COMMCN/CC  MHO 1/NC8J,NFEICR,NVARr NMAX, NCYCLE 
COMMCH/CCMH02/X (2500) , XBASE  (2500) 

COMMCN/CC  KM  03/ A  (10)  ,AEIST  (10) 

CQHMCN/CCMM04/ACC  EL,RECUCI,DELTAfIRPRCV,IPRINT,ITEPM 
COMHCN/CCaaOS/LEE  (2500) ,UBD  (2500)  rEPS  (2500)  ,IES  (2500) 
COMECN/CCHE06/CBJ  (2500) ,EHS  (2500) 

COHMCN/CCSM07/N  (2500) , P  (2500) 

COHMCN/CCMMO 8/EPSY  (10) 

COMMCN/CC KB09/SA VI 

EEGIN  THE  PROGRAM 
CALL  DATAIN 

ALL  NECESSARY  DATA  IS  INPUT  AND  PRINTED  CUT 
CALL  HJAIG 

EEGIN  THE  PATTERN  S1ABCH  ALGCPITHM 
ALL  FUFTHEB  KCRK  IS  COMPLETED  IN  TBE  SUBROUTINES 
CF  THE  CODE 
STOP 
END 
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SUBROUTINE  DATAIK 

**** ********************** * *********** *************** ****** ***** ****** 

SUBROUTINE  DAT AIN  BEADS  IN  AIL  THE  NECESSARY  DATA 
AND  PRINTS  IT  CUT  AS  A  MEANS  CE  CHECKING  IT. 


INTEGFB*2  IMPBOV,IPRINl,ITERM , IPRCG ,IACCEL 
REA1*4  IPS, LED 

COMKON/COSKO 1/NCBJ, NPBIOB,NVAB,NMAX, NCYCIE 
COMMCN/CCEH02/X (2500) , XE ASE (2500) 

COEMCN/CCKM04/ACCEL, RE IUCE , DELTA / IMPRCV ,IPRINT,ITEBM 
COMMCN/CCEM05/LBD (2500) ,UED (2500) ,EPS (2500) ,IPS (2500) 
CCMHCN/CCMH06/OEJ (2500) ,RHS (2500) 

COMMCN/CCEH08/EPSY ( 10) 

C  READ  AIL  THE  NECESSARY  DATA 

READ  170,NOBJ,NPRIOF,NVAR,NMAX, IPRINT 
READ  175,  (X(K)  ,K=1,  KVAE) 

READ  175,  (EPS(K)  ,K=1,NVAR) 

READ  175,  (LED  (K)  ,K=1,NVAB) 

READ  175,  (UBE  (K)  ,K=1,NVAB) 

READ  175,  (RHS(I)  ,I=1,NCEJ) 

READ  175,  (EPSY  ( J)  , J=  1 ,  f  PRIOR) 

READ  175,ACCEI, REDUCE, EELTA 
C  PRINT  AN  ECHO  CHECK  CF  THE  DATA 

100  PRINT  1 80, NCFJ, NPRICF, AVAR, NH  AX , ACCEL, REDUCE, DELTA 
PRINT  181 
DO  110  K*1, NVAB 

110  PRINT  182,K,X(K)  ,K,EPS  (K)  ,K,LBD  (K)  ,K,UBD(K) 

PRINT  183 
DO  120  1*1, NCBJ 
120  PRINT  184,1, BHS (I) 

PRINT  185 
DO  130  J*  1, N EBIOR 
130  PRINT  1 86  ,J,  EPSY  ( J) 

PRINT  187 

IF  (IPRINT  .EC-  1)  GC  TC  140 

PRINT  188 

RETURN 

140  PRINT  189 
170  FORMAT  (1615) 

175  FORMAT  (8F10.0) 

180  FORMAT  ('-',25X, 'NUMBER  CF  OBJECTIVES  =' ,T57,I5,/, '  •  ,25X ,' NUMBER  0 
$F  PRIORITIES  *' ,157,15,/,'  *, 25 X ,' NUMBER  OF  DECISION  VARIABLES  =  ', 
JT57,l5,/,«  *  ,25X, ' EAXIECM  NUMBER  CF  AIIOHABLE • ,/, '  •  ,29X ,' PATTERN 
$  SEARCH  CYCLES  * * ,157 , 15,/, *  * ,25X,' ACCEL ERATICN  FACTOR  =',T59,P10. 
$ 6,/, •  * , 25X, *STEP  SIZE  REDUCTION  FACTOR  =  • ,T59,F10.6 ,/, '  ',251, 'MI 
SNIHUH  ALLOWABLE  STEP  SIZE  « »,T5 9, F 10. 6) 

181  FORMAT  l*-*, 251, *  VICTOR  CF  INITIAL  ESTIMATES  FCR  DECISION  VARIABLES 

$ - X*,/,’  * ,25X,*PEBTUIEATI0N  STEP  SIZES  FOB  THESE  VARIABLES - EPS 

J'#/.'  ' , 25X, ' LORES  AND  UPPER  BOUNDS  FCR  DECISION  VARIABLES - L3D  A 

$ ND  USD') 

182  FORMAT  ( '  » ,  1 0X,  «  X  {'  .15  , ' )  =  *  ,  F 10. 5, 7X, '  EPS  (•  ,  1 5,  •  )  =  ',F10.5, 

$7X,  'LBD  ('  ,15,*)  *  ',F10.5,7X,'0ED(*  ,15,')  *  *,F10.5) 

183  FORMAT('-»,25X, 'RIGHT-HAND  SIDE  VALUES  FCR  OBJECTIVE  FUNCTIONS - R 

$  HS ' ) 

184  FORMAT  ( *  •  ,2 51, 'R BS  ( ' , 15 , ♦)  =  »,F20.4) 

185  FORMAT  (»-',25X, 'VECTCF  CF  ACHIEVEMENT  FUNCTION  TOLERANCES - EPSY') 

186  FORMAT  ('  ' ,25X, * EPSY  ( •  ,12 , » )  =  ',F10.6) 

187  FORMAT  ('-',25X,' AIL  DATA  HAS  BEEN  INPUT - BEGIN  THE  PATTERN  SEARCH 
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S') 

188  FORMATf  * , 25X, • NCTE:  IFPIHT  =  0,  SO  ONLY  THE  FINAL  SOLUTION  WILL 
SHE  PRINIEE  COT') 

189  FORMAT (*  ',25X,'NCTE:  IFFTNT  =  1,  SC  THE  SOLOTICN  WILL  EE  PRINTED 
SOOT  AT  EVERY  PATTERS  SEARCH  CYCLE') 

RETURN 

END 

C  ********************************************************************** 


SUBROUTINE  HJALG 

C  ****  ******************************* ****** * ** ****** ********* *********** 
C 

C  SUBROUTINE  HJAIG  IS  THE  BOOKE-JEEVES  FATTEFN 

C  SEARCH  ALGCFITHK.  IT  FERFCRM5  EXPLORATORY  AND 

C  PATTERN  HCVES  TC  FIND  THE  'BEST'  SOLUTION  TO  THE 

C  NONLINEAR  GCAL  EECGEAHHING  PRCELEH. 

C 

C 

REAI*4  IPS, LED 

INTEGEB*2  IMFRCV, IFRI NT , ITERM, I PFCG , I  ACCEL 
COM H ON/C C MHO 1/NCBJ,NPRIOR,NVAR,NEAX,NCYCLE 
COMECN/CCKM02/X  (2  500)  ,  XEASE  (250  0) 

COMMON/COMM 03/ A  (10) , AEEST (1 0) 

COMMCN/CCMB04/ ACCEL, RE  LUCE,  DELTA, IMPRCV, IPRINT, ITERM 
COMHCN/CCMM05/LBD (2500) , UED (2500)  ,EPS  (2500)  ,IFS  (2500) 

COHHCN/CC MHO 9/SAVE 
CALL  ACHFUN 

SUBROUTINE  ACHFUN  FCRSS  THE  ACHIEVEMENT  VECTOR  FOR 
T BE  INITIAL  BASE  POINT.  IT  IS  TESTED  FOR  ALL  ZERO 

VALUES - IF  SC,  THE  INITIAL  PCINT  IS  OPTIMAL. 

OTHERWISE,  IT  IS  SET  AS  THE  BEST  VALUE  AND  IS  PRINTED. 

DO  100  1 , NFRIOR 

IF  (A  (J)  .NE.  0.0)  GO  TC  203 
100  CONTINUE 

ALL  PRIORITY  LEVELS  ARE  SATISFIED  AT  ZERO  VALUE. 

THE  INITIAL  FCINT  IS  CFTIMAL - TERMINATE. 

200  PRINT  290 
DO  201  J=1,NFBICB 

201  PRINT  292, J, A  (J) 

PRINT  291 
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DO  202  K=1,NVAB 

202  PRIM  280, F,X  (K) 

GO  10  295 

203  PRINT  281 

DO  204  J=1, NFEIOF 
ABEST  (J)  =  A  (J) 

204  PEIM  282, J, A  (J) 

**  INITIALIZATION  Cl  PATTERN  SEARCH  PARAMETERS  ** 

NCYCIE  COUNTS  THE  NCMEER  OF  PATTERN  SEARCHES  0SED 
I ACC El  IS  AN  I KDIC ATCE  WHICH  TEILS  IF  THE  EXPLORATORY 
HCVES  ARE  HADE  AEC0T  AN  ACCELERATION  POINT 
NCYCIE  =  0 
IACCE1  =  0 

SET  THE  IHITIAI  EASE  PCINT  AND  THE  INITIAL  STEP 
SIZE  EEEOBE  THE  EXPICR ATCB Y  MOVES 
DO  2C5  K=1,NVAR 
XBASE  (K)  =  X  <K) 

205  IPS  { K)  =  EPS  (K) 

START  THE  PATTERN  SEARCH 
**  EXPICRATCRY  MCVES  ** 

206  IF (NCYCIE  .GT.  NH  AX)  GC  TO  260 
IPRCG  =  0 

IFROG  COUNTS  TEE  KOFBEE  OF  IMPROVEMENTS  DURING  THE 
EXELCRATGRY  MCVES.  IT  IS  USEE  TO  DECIDE  IE  A  PATTERN 
HCVE  SHOULD  EE  PEEFCRKED. 

207  DO  210  K=1,NVAB 

PEP.TUR EATICNS  IN  THE  POSITIVE  DIRECTION 
TEST  STEP  SIZE  CCKECNENT.  IF  LESS  THAN  THE  MINIMUM 
ALLOWABLE  STEP  SIZE,  A  RIDGE  SEARCH  WILL  BE  PERFORMED 
IF  (IPS  f K)  .IT.  DEITA)  GC  TO  255 
X  (K)  «  X(K)  ♦  IPS  (K) 

IF  {XI K}  .IT.  IBD  (K)  )  GC  TC  208 
IF  (X  (K)  .GT.  UBC  (K)  )  GC  TO  208 
CALL  ACHFUN 

ACHFUN  FORMS  THE  ACHIEVEMENT  VECTOR  AT  THE  TEST  PCINT 
ITERS  =  0 
CALI  DECIDE 

DECIDE  COMPARES  THE  ACHIEVEMENT  VECTORS  OF  THE 
EASE  PCINT  AND  THE  TEST  POINT-  IF  X  (K)  IMPROVES  THE 
SCIUTICN ,  IT  BECOMES  A  TEMPORARY  HEAD  POINT  AND 
EXPLCR ATOEY  MOVES  AFE  CONTINUED  FRCM  IT.  IF  THEBE 
IS  NC  IMPECVEME NT ,  TEEN  EXPLORE  IN  THE  NEGATIVE  DIRECTION 
IF (IMPBOV  -EC-  0)  GC  TC  208 
IPRCG  =  IPRCG  ♦  1 
GO  TC  210 

PERTH R EATICNS  IN  THF  NEGATIVE  DIRECTION 

208  X  (K)  =  X(K)  -  2. 0*IFS  (F) 

IF  (X  (K)  -IT.  LBD(K))  GC  TC  209 
IF  (X  (K)  .GT.  UBD(K))  GC  TC  209 
CALL  ACHFUN 

ACHFUN  FORMS  THE  ACHIEVEMENT  VECTOR  AT  THE  TEST  POINT 
HERB  =  0 
CALL  DECIDE 

IF  X(K)  IMPHOVFS  THF  SOLUTION,  IT  BECOMES  THE 
TEMPORARY  HE AC  PCINT  AND  THE  NEXT  V J PIABLE  IS 
PESTUBEEC.  IF  THEBE  IS  NO  IMPROVEMENT,  XBASE  (K) 

REMAINS  THE  TEMFCEAFY  HEAD  PCINT. 

IF  (IMFROV  .EC-0)  GO  TC  209 
IPBCG  =  IPRCG  ♦  1 
GO  TC  210 
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I?  (IFRINT  -EC-  0)  GO  1C  206 
PR1 HI  287, NC YCLE 
DO  241  K-1, NVAB 
241  PRINT  280, K, X  { K ) 

PRIKI  288 
DO  245  3= 1, NFRIOR 
245  PRINT  282,3, A  (3) 

GO  TC  206 

EXPLORATION  ARCUKE  1HE  ACCELERATION  POINT 
DOES  NOT  HELD  ANY  IMPROVEMENT 

BESET  THE  EASE  ECIN1  AND  GO  TC  EXPLCF ATOBY  MOVES. 

250  IF  (IACCEL  .EC-  0)  GC  1C  255 

251  DO  252  K=1,N¥AR 

252  X  (K)  =  XBASE  (K) 

IACCEL  =  0 

GO  TC  206 

STEP  SIZE  REDUCTIONS  ECUALLEE. 

A  RESCIDTICN  RIDGE  SEARCH  WILL  BE  FERFORKEE 
255  PRINT  284 

259  CALI  RIDGE 

SUBROUTINE  RIDGE  ill!  TRY  TO  FINE  A  NEK  BASE  POINT 
IN  AN  CBLICUE  EIFEC1ICN  TO  THE  PRESENT  EXFLCRATCBY 
SEARCH  DIRECTICK.  IE  SUCCESSFUL,  THE  PATTERN 
SEARCH  Will  BEGIN  AT  THIS  NEW  POINT.  IF  UNSUCCESSFUL, 

THE  ALGORITHM  Kill  EE  TERMINATED. 

IACCEL  =  0 
GO  TC  206 

THE  MAXIMUM  NUMBER  CE  PATTERN  SEARCH  CYCLES  HAS  BEEN 
EXCEEDED.  THE  EES1  VALUES  OF  THE  VARIABLES  AND  THE 
ACHIEVEMENT  VECTOR  ABE  PRINTED  CDT ,  AND  THE  ALGCFITH M 
Kill  EE  TERMINATED. 

260  PRINT  289 

IF  (IACCEL  .LT. 1)  GC  1C  265 
EO  261  K= 1, NV  AB 

261  PRINT  280, K, X  (K) 

PRINT  288 

DO  262  3= 1 , NEFI OR 

262  PRINT  282,3, A  (3) 

GO  TC  295 

265  DO  266  K=1,NVAR 

266  PRIKI  280 ,K, XEASE  (K) 

PRINT  288 

DO  267  3=1, NPRICR 

267  PRINT  282,3, ABEST (3) 

270  GO  IC  295 

279  FORMAT(,-',2EX,,THE  EXPICRATORY  HOVE  HAS  BEEN  SUCCESSFUL',/,'  *,25 
$ X, ' A  PATTERN  HCVE  KILL  EE  PERFORMED') 

280  FORMAT  (*  '  ,  241,  •  X  ('  ,15  , ' )  =  ',F13.6) 

281  FORMAT ( *  ',2SX, 'INITIAL  ACHIEVEMENT  VECTOR  V1IUES') 

282  FORMAT  (*  '  ,  34X,  'A  (',12,')  =  *,F13.6) 

283  FOHHAT(»-*,2SX, 'TEMPORARY  HEAD  POINT  DCES  NOT  IMPROVE  THE  SOLUTION 
$',/,'  '  , 25X,  'THE  STEP  SIZES  KILL  EE  BFEUCED',/,'  ',251, 'NEB  EXPLCR 
SATOFY  MOVES  HILL  EE  PEFFCFHED  ABOUT  XBASE') 

284  FORMAT  ( •  '  ,251, ' MAXIHOE  STEP  SIZE  REDUCTIONS  EXCEEDED  FOR  THIS  EXP 
SLORATION ' ,/, '  * ,25X, • I  RIEGE  SEARCH  WILL  BE  PERFORMED') 

295  FOBHAT('-*,25X,*THE  ACCELERATION  POINT  YIELDS  AN  IMPROVED  SOLUTION 
S’,/.'  • ,25X, 'EXPLCFATCFY  MOVES  KILL  BEGIN  AT  THIS  NEK  POINT') 

286  FORMAT  (*-*,25X, 'THE  ACCELERATION  POINT  DCES  NCT  YIELD  AN  IMROVED 
S  SOLUTION  * ,/ , *  * ,251, 'SXEICBATORY  MOVES  HILL  BEGIN  AT  THIS  POINT*) 

287  FORMAT  (•-• ,251, 'THE  BEST  SOLUTION  AT  PATTERN  SEARCH  NUMBER  ',13,' 
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5  FOLLOWS • ) 

288  FORMAT  ( '  *,2EX,*TBE  ACHIEVEMENT  VICTOR  V ALOES  APE:*) 

289  FOEMAT  (*  *  ,25X,  *  MAXIMUM  NHMEEB  OF  PATTERN  SEARCH  CYCLES  EXCEEDED', 
$/,'  *,25Z,»TBE  ALGCBIT EM  IS  TERMINATED*,/,'  *,25X,'THE  BEST  SCLUTI 
SON  TC  THIS  FCINT  FCLLC IS ' ) 

290  FORMAT('— *,20X,43 ( 1 H* ) ,/,'  ', 25X, 'OPTIMAL  ACHIEVEMENT  VECTOR  VALUE 
is*,/,*  *,20X,43(1H*)) 

291  FORMAT (»-', 201,43 (IB*) ,/,*  *, 25X , 'OPTIMAL  DECISION  VARIAELE  VALUES 
S*,/,*  '  ,  20X,  43  { 1H*)  ) 

295  STOP 
END 

C  **************************  *  ****  **********  ******************  *********** 


SUBROUTINE  DECIDE 

C  ****  ****************************************  *****************  ********* 
C 

C  SUBROUTINE  DECIDE  IE  OSED  FOR  TWC  FUFFOSES: 

C  (1)  TO  COMPARE  TEE  ACHIEVEMENT  VECTOR  VALUES  OF 

C  THE  BASE  POINT  AND  A  TEST  POINT  DURING  THE  EXPLORATORY 

C  MOVES  AND  SELECT  THE  BEST  POINT  FOR  FUTURE  MOVES 

C  (2)  TO  TEST  AN  ACHIEVEMENT  VECTOR  FOR  TERMINATION 

C  CF  THE  ALGORITHM.  THIS  CAN  OCCUR  IN  TWO  METHODS: 

C  ****TER HI  NATION  El  MEANS  CF  ZERC  FOR  ALL 

C  EEICRITY  LEVEIS,  AND 

C  ****TERHINATIGK  BY  MEETING  SPECIFIED  TOLERANCES 

C 

C  THE  PARAMETER  ITEFM  IS  AN  INDICATOR  FCR  TERMINATION 

C  TESTING.  IF  ITEFM  =  0,  ONLY  THE  COMPARISON  TEST  IS 

C  PERFORMED.  IF  ITEBE  =  1,  THE  TERMINATION  TESTS  WILL 

C  ALSC  BE  PERFORMED. 

C 

c 

INTEGEB*2  IHPBOV, IFRINT,ITERH 
COHHCN/CCHMO 1/NOE J, RPR JCF,NVAR,  NM AX, NCYCLE 
COM MCN/CCMMO  2/X  (2500)  ,XEASE(2500) 

COMHCN/CCRH03/A  (1.0)  ,  AEIST  (10) 

COM MCN/CCMMO 4/ ACC EL ,FEIUC,DELTA,IMFFCV,I FRI NT, ITERM 
COMHCK/CCMMOS/EPSY (10) 

C  COMPARISON  CF  ACF.3EVE5ENT  VECTORS  FOR  EXPLCEATCRY  MOVES  ** 

600  DO  610  J* 1, NPEIOR 

IF  (A  (J)  .GT.  AB2ST  (J)  <0.0001)  GC  TC  630 
IF  (A  (J)  <0.0001  .IT.  AEIST  (J) )  GO  TC  620 
IP  (J  .EC.  NPRICB)  GC  TC  630 
610  CONTINUE 


J 
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C  THE  TEST  FCINT  IRPECVES  THE  SCIUTICK 

620  ISPBCV  =  1 
GO  TC  640 

C  THE  TEST  FCINT  DOES  NCT  IMPROVE  THE  SCLUTICN 

630  IMFPCV  =  0 

C  CHECK  THE  INDICATOR  ECU  FURTHER  TERRINATICK  TESTING 

640  IP  (ITEBB  .EQ.  0) GC  TC  670 
CONTINUE 

*♦  TEBEIN1TICN  TESTS  ** 

( 1)  T EFEINATICN  TEST  E  Y  MEANS  OE  ZERC  FOB  AIL  PRIORITIES 
650  DO  655  J=1,NFFICF 

I?  (A  (J)  -NE.  0.0)  GC  TC  660 
655  CONTINUE 

ALL  FRICHITY  LEVELS  ABE  ZERO - THE  OPTIMAL  SCLUTICN 

HAS  BEEN  FCUND - TEFRINATE  THE  ALGORITHM 

GO  TC  680 

(2)  TEREINATICN  TEST  El  TOLERANCES 
660  DO  665  J=1,NFFIOB 

IF  (ABS  (A  (J)  -  ABEST(J))  .GT.  EPSY(J))GC  TO  670 
665  CONTINUE 

ALL  TOLERANCES  ARE  SATISFIED - A  SOLUTION  HAS  EEEN 

FCUND - TEE  SEABCB  ALGORITHM  IS  COMPLETE 

GO  TC  680 

A  TOLERANCE  IS  NCT  SATISFIED - A  BETTER  SCIUTICN 

KILL  EE  SOUGHT - RETURN  TO  THE  PATTERN  SEARCH 

670  IFfIRFROV  -EC.  0)  GC  TC  695 
DO  675  J=1, NPFIGE 
675  ABESTfJ)  =  A  (J) 

GO  TC  695 
680  PRINT  690 
PRINT  691 
DO  6 £2  3~  1, NEEIOB 
682  PRINT  692, J,  A  (J) 

PRINT  693 
DO  665  K=1,NVAB 
685  PRINT  694, K, X  (K) 

STOP 

690  FOHMAT('-',25X,  'THE  SEARCH  ALGORITHM  IS  COMPLETE - ALL  TCLEF.ANCES 

S  FOR  A  SCIUTICK  ARE  SATISFIED*) 

691  FOBMAT(*-* ,20X, 43  ( 1 H * )  ,/,*  *, 25X, 'OPTIMAL  ACHIEVEMENT  VECTOR  VALUE 
SS*,/,'  * ,20X , 43  ( 1 H*) ) 

692  FOBMAT  (•  *  ,  34X  ,  •  A  J*  ,  12  ,  •)  =  ',F13.6) 

693  FOBMAT(*-',20X,43  (1H*)  ,/, *  * , 25X  ,  *OPTI RAL  DECISION  VARIAELE  VALUES 
• ,20X,43  (1H*) ) 

694  FORMAT  ('  ',341,'X  ('  ,15,'J  =  *,F13.6) 

695  RETUHN 
END 

C  ********************************************************************** 
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SUBROUTINE  R1IGE 

C  ***********************4**44****************************************** 

c 

c  SOBRCOTINE  RIDGE  IS  CALLED  WEEK  THE  PATTERN  SEARCH 

C  CAN  NC  LCNGER  EINC  EATTESN  OR  EXPLORATORY  COVES 

C  fiHICH  IHFECVE  THE  ACHIEVEMENT  VECTCF.  RIDGE 

C  EVALUATES  EXPLCEATCBY  POINTS  VHICH  ABE  IN  CELICUE 

C  DIRECTIONS  TC  THE  DSDAL  EXPLORATION  AXES.  RIDGE 

C  ATTEMPTS  TO  FIND  A  EESCLOTION  RIDGE,  IF  IT  EXISTS, 

C  AND  MOVE  THE  PATTERN  IN  THAT  DIRECTION. 

C 

C 

DIMENSION  APCS{10) ,ANEG (10) 

REA  1*4  IPS, LED 

INTEGER  *2  IMEBOV, IPRINI,ITERH 
COMMCN/CCEHO 1/NCE J, NFE3CF,NVAR, NM  AX , NCYCLE 
COMMCN/CCMMO  2/X (2500) , XBASE (2500) 

COMKCN/CCCM03/A  (1 0)  , AE ESI  (10) 

COM MCN/CCMM04/ ACCEL, BEIDCE, CELT A, IMPPCV  ,IPRINT , ITEBM 
COMMCN/CCBM05/LED (2500) ,UED (2500)  ,EPS  (2500) ,IFS (2500) 
CCMMCN/CCKMOS/SAVE 

FERTOEE  EACH  VABIAEIE  CNE  AT  A  TIME 
PERTOKE  IN  THE  POSITIVE  CIRECTICN 
DO  800  K= 1 , N V AB 
300  X (K)  =  XBASE  (K) 

801  DO  865  K= 1 , N V  AR 

X(K)  =  X  (K)  4  IPS  (K) 

IP  (X  (K)  .IT.  IBD  { K)  )  GC  TC  805 
IF  (X  (K)  .GT.  UBD  ( K)  )  GC  TC  810 
GO  TC  815 
805  X (K)  =  LBD(K) 

GO  TC  815 
810  X  (K)  =  OEC(K) 

815  CALI  ACHFON 

ACHFUN  FORMS  THE  ACHIEVEMENT  VECTOR  IN  THE  POSITIVE 
DIRECTION 
DO  820  J=  1 , NFEIOR 
820  APOS(J)  =  A(J) 

PERTORE  IS  THE  NEGATIVE  DIRECTION 
X  (K)  =  X  (K)  -  2.0  *  IPS  (K) 

IF  (X  (K)  .IT.  LBD(K))  GC  TC  825 
IF  (X  (K)  .GT.  UBD(K))  GC  TC  830 
GO  TC  835 
825  X (K)  =  LEE (K) 

GO  TC  835 
830  X  (K)  =  OED(K) 

835  CALL  ACHFON 

ACHFON  FORMS  THE  ACHIEVEMENT  VECTOR  IN  THE  NEGATIVE 
DIRECTION 
DO  840  J  =  1, N PRICE 
840  ANEG(J)  =  A  (J) 

COMPARE  APCS(J)  AND  ANEG(J).  SEIECT  THE  MINIMUM 
VIA  A  STEFHISE  FREE  EFTIVE  PROCEDORE.  THIS  DETERMINES 
A  DIRECTION  FCB  THE  PIDGE  POINT. 

DO  850  J= 1 , NERICR 

IF ( A POS ( J)  -IT.  A SEG  ( J) ) GC  TO  855 
IF (APOS  (J)  -GT.  AN EG  (J) ) GC  TO  860 
IF (J  -EQ.  NPFICR) GC  1C  e55 
850  CONTINOE 

855  X  (K)  =  XEASE(K)  ♦  IPS  (F) 
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GO  TC  661 

860  X  (K)  =  XEASE(K)  -  IPS(K) 

RETURN  X  { K)  1C  IIS  IASI  POINT  VALUE,  XBASE,  BEFORE 
THE  PEETDBEATICN  OF  1BE  NEXT  VARIABLE 

861  SAVE  =  XEASE  IK) 

XBASE  (K)  =  X  (K) 

865  X  (K)  =  SAVE 

TEE  NEB  PCINT,  XRIDGE,  HAS  BEEN  FORMED  IN  AS  CBLICOE 
DIRECTION  TO  THE  CLI  EXPLORATORY  MCVES.  ITS 
ACHIEVEMENT  VECTCB  IS  COMPUTED  AND  COMPARED  TO  THE 
BEST  VALUE  AT  PCI  NT  XEASE. 

CO  866  K= 1, N  V AR 
SAVE  =  XEASE  IK) 

XBASE  (K)  =  X  (K) 

866  X(K)  =  SAVE 
CALL  ACHFUN 
ITERE  =  0 
CALI  DECIDE 

IF  THE  RIDGE  PCINT  GIVES  NO  IMPROVEMENT,  TEE  PRESENT 
EASE  PCINT  IS  T BE  EIST  SOLUTION,  AND  THE  AIGCRITHH 
IS  TERMINATED. 

IF  THE  RIDGE  PCINT  C-IVES  IMPROVEMENT,  AN  EVEN  EETTER 
PCINT  RILL  BE  SOUGHT  ALONG  THE  RESOLUTION  RIDGE. 
IFfIMPROV  .EC.1)  GO  TC  £69 
CONTINUE 

XEASE  IS  THE  BEST  SOLUTION 
PRINT  891 
PRINT  897 
DO  £67  K=  1 , N V AR 

867  FRINT  892, K, XEASE  (K) 

PRINT  898 

CO  868  J=1,NPEICR 

868  PRINT  893,J, AEEST  (J) 

STOP 

THE  RIDGE  ECINT  GIVES  AN  IMPROVED  SOLUTION. 

TEI  TC  ACCELERATE  ALCNG  THE  RESOLUTION  BIDGE 

869  PRINT  894 

870  DO  871  K=1,NVAB 
SAVE  =  X  IK) 

X  (K)  =  X(K)  ♦  ACC  El  *  |X{K)  -  XBASE  fK)  ) 

871  XBASE  (K)  =  SAVE 
CALI  ACHFUN 
ITERM  *  0 

CALI  DECIDE 

IF  THE  ACCELERATION  ECINT  YIELDS  AN  IMPROVED  SOLUTION, 
TRY  TC  ACCELERATE  F  AETHER  ALCNG  THE  RESOLUTION  RIDGE 
IF  THERE  IS  NC  IMPECVEMENT,  THE  RIDGE  POINT  BECOMES 
THE  NEB  POINT  FCR  TEI  RENEWED  PATTERN  SEARCH. 

IF  (IEERCV  .EC.  0)  GC  TC  872 
PRINT  896 
GO  TC  870 

872  PRINT  895 
PRINT  897 

DO  873  K=  1, N V AR 

873  PRINT  892, K, XEASE  (K) 

PRINT  898 

DO  875  J=1,NEBICB 
875  PRINT  893, J, AEEST  (J) 

RESET  THE  EASE  PCINT  AND  STEP  SIZE  VECTORS  AND 
BSTUSN  TO  THE  1XPLCF ATCRY  MOVES. 
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C  NC  IBFECVEKENT  FECK 

DO  877  K=1,NVAB 
X  (X)  =  XEASE(K) 

877  IPS  (K)  =  EPS  | K) 

DO  879  J=  1,  NE BIOP 

879  i(J)  =  AEEST(J) 

GO  TC  899 

891  FOPEAT ('-',25X, 'THE  El 
SED  SOLUTION  CAS  EE  FCU 
$  BESI  SOLUTION  THAI  CA 

892  FORMAT ('  •  ,34X, • X  ( • ,15 

893  FOES  AT  ( '  •  ,  3  111,  •  A  (•  ,  12 

894  FOSE AT  ('-',25X,'THE  El 
S25X,'TRY  TC  ACCELERATE 

895  FORMAT  ('-',2SX, 'ACCELE 
SIMPHCVED  SOLUTION ' ' 
$  T* ) 

896  FOBKAT('-',25X, 'ACCELE 
S IHPECVED  SOLUTION  *  ,  /  ,  ' 

397  FOBEAI  ('-',25X, 'EECISI 

898  FOBEAT  ('-' ,25X, • ACHIEV 

899  RETURN 
END 

C  **** ********************** 


ACCELEBATICN  PC3NT 


EGE  SEARCH  IS  COMPLETE’ ,/, *  ',25X,'NO  IMPHCV 
ID',/,'  * ,25X,*THE  CUERENT  BASE  POINT  IS  THE 
l  EE  FOUND  THROUGH  THE  PATTERN  SEARCH') 

,')  =  '  , F 13. 6) 

#  ')  =  •  g  F13.  6) 

EGE  PCINT  GIVES  AN  IMPROVED  SCLUTICN' , /, •  ', 
A1CNG  THE  FESCIUTION  RIDGE') 

FATION  ALCNG  THE  RESOIOTICN  RIDGE  HELDS  NC 
• , 25X,  'THE  RIDGE  POINT  IS  THE  NEW  BASE  POI1I 

FATION  ALONG  THE  RESOLUTION  RIDGE  YIELDS  AS 
* ,25X,'TBY  TC  ACCELEEATE  FURTHER') 

CS  VARIABLES') 

EBENT  VALUES') 


4  ***************** ************* 


SUBROUTINE  DEWAR 

*******  4********* ** 4***  4 **** ********** *44 ************* 

SUBROUTINE  DEWAR  CAICULATES  THE  POSITIVE  AND 
NEGATIVE  DEVIATION  VARIAELES  FCB  EACH  OBJECTIVE 
FUNCTION  USING  THE  PREVIOUSLY  EVALUATED  OBJECTIVE 
FUNCTIONS  AND  THE  RIGHT-HAND  SIDE  VALUES. 


REAI*4  N 

COMHCN/CCKHOI/NCEJ, EEB1CF, NVAR, NK AX ,NCYCLE 
COMBCN/CCHM02/X (2500) , XEASE (2500) 

COHBCN/CCBE06/CEJ (2500) , RHS (2500) 

COMBCN/CCMM07/N (250C) ,F (2500) 

CALL  CBJFUN 
400  DO  430  I=1,NCEJ 

IF  (CBJ  (I)  -  RHS  (I) )  410  ,420,420 
410  N  (I)  =  RHS  (I)  -  CEJ  (I) 

P  (I)  =  0.0 
GO  TO  430 
420  N (I)  =  0.00 

P(I)  =  OBJ  (I)  -  RHS  (I) 

430  CONTINUE 
RETURN 
END 

*************************  *  *  ************** ** ********  ********  *********** 


! 


no  no  no 
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SOBECUTINE  ACHPUN 

**** **********************  *********  **************** ********  *********** 


REAl*«  N 

COf!flON/CCBHO  1/SIOBJ,  NPRIGB ,NVAR,  NKAX,NCTC1E 
COMHCK/CCEH02/X  (2500)  ,XEASE  (2500) 

CO8KCN/CCKS03/A  (10)  ,  AEEST  (10) 

COHMCN/CCEB07/N (2500) , E  (2500) 

CALI  EEVVAB 

ACHIEVEMENT  FUNCTIONS  ABE  ENTERED  HEBE. 

RETURN 

END 

**** ***********4  4*** ****** ******* **************************  *********** 


SUBROUTINE  CEJFOK 

******* ********** ********4  4 ** ************  **********************  ******* 

CCHflCK/CCBBO 1/NCEJ, EFF1CF,NVAR,NEAX, NCYCLE 
C03BCN/CCHH02/X (2500) , XBASE (2500) 

COHHCN/CCBE06/CEJ  (2500) , BBS  (2500) 

C  CEJECTIVE  FUNC1ICNS  ABE  ENTERED  HEBE 

RETURN 
END 


I 

I 

I 

I 

1 

J 

■» 


APPENDIX 


NLGP/MPS-RS  CODE: 


C 


VERSION  3 


Note:  Version  3  is  set  up  to  solve  the  Target  Allocation  Problem. 


I 

I 

I 

I 

] 

1 


i 

I 


i: 


c 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 
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nonlinear  gcal  FEC3EAMi!iNG--p;.Ti  ers  siafcf  code 

*****:s»«r*«*«:*4?$**fc#*********************«*«**»»«*t«»***«»«*4* 


THIS  15  1 H r  NLCr/MPS-PS  XAVl  CCDI 

NONLINEAR  GCA  1  EF.  CG  IAMBI  NG/MCtlllED  FATTRRK  SEARCH  CODE 
WITH  A  RESOLUTION  RIDGE  SEARCH  TECHNIQUE 
VERSION  3:  SUERCGTI M  UVA10E 
-  JULY  IE,  1E£0  - 


*  *  * CODE  SEECIrICATICNS*** 
2500  OBJECTIVES 
25  CO  VAFIAEIES 

10  F5 IC  E 1 1 Y  IEVELE 


*  *  *V A  r  I  A  E  L  E  S  AND  PARAMETER  £*** 


HCEJ 
HP  El  OR 
5  VAR 

n 

H C  iCLS 
X 

XBASE 

A 

AB  ESI 
AFC^,A.*EG 
EPS, IP  S 

L3D,  03 D 

CEJ 

RHS 

N,  E 

EPSY 

ACCEL 

REDUCE 

DELTA 

IKFROV 


IPFINT 


I  TERM 


IFFCG 

IACCEL 


SAVE 

JFI 


***DESCFI ETI C  N*  ** 

NUMBER  CE  OBJECTIVES 
NOSEER  OE  PRIORITY  IEVELS 
NUMEEB  CF  EECISICF  VfFIAELES 

“AXir.ua  KC.  CF  PATTERN  SEARCH  CYCLES  ALICKED 
CCCKIEF  FCF  ESTTEFN  Sr/RCH  CYCLES 
I EC1S10K  VARIAEE1  VECTCE  USED  FCF  EEFTU5ESIICS 
EASE  F  Cl  NT  VECTCF  FCF  EATTERN  SIARCF 
AC  El EV  EM  FN  T  VECTOR  AT  TEST  PCINT 
ACHIEVEMENT  VECTCF  ICF  EEST  SCIUTICN 
T  El'POR  AE  Y  AC  SIEVE  S  EFT  VECTCF  S  FCB  FIDGE  SEACH 
SCAIAFS  CF  PFFTDFI/11CF  STEP  SIZES  USED  FOR 
EXFLCR ATOF Y  MCVES 

I  CNF  £  AND  UFEEF  ECUM  VECTORS  ICR  IICISICN  VAR 
CEJECT3VE  FUNCTION  VECTOR 

RIGHT-HARE  SIEE  VA1DE  VECTCP  ECF  CEJECTIVE  FN  S 

NEGATIVE  AND  POSITIVE  DEVIATION  VAFIA5IE  VECT. 

1EEMINATICA  TCIEEANCE  TEST  VECTOR 

ACCEL  ERA T ICN  FACTOR 

SIEE  SIZE  FEEUCTICN  IACTCR 

MINIMUM  ALLOWAELE  STEP  SIZE 

TEST  SWITCH - II  IfIFCV=l,  X  IMPROVES  THE 

SOLUTION;  IF  IM?P.CV  =  C,  X  DOES  NCI  ItEFCVE  THE 
SCIUTICN 

OUTPUT  SWITCH - IF  1EF.INT=1,  RESULTS  ARE 

EFIN1EC  CUT  AT  EVEFY  EATTERN  SIARCE  CYCLE; 

IE  I?F,  IN  1  =  C,  CNIY  TFE  FINAL  FESCITS  AFE  OUTPUT 

TEFMIKATICK  TEST  SW1TCI - IF  !TE?:i  =  1,  THE 

SCLUIICK  IS  TESTED  ECR  TEF.  HI  NATION  ;  IF  ITE  F  K=0 
CMY  A  CCMFAFISCI  T  1ST  IS  PER  ICR  ME  I 
COUNTER  FOR  IMERC  VEFENTS  DURING  THE  EXPICEA- 
TCEY  MOVES 

INDICATOR  FOR  ACCEL  EF.  A  II  CN  MCVES: 

IACCEL  =  1,  ACCEIEF  AT ICN  MCV  E  WAS  SUCCESSFUL 
IACCEL  =- 1,  ACCELERATION  MOVE  NAS  UNSUCCESSFUL 
IACCEL  =  0,  ACCEIEF  IT  ICN  MOVE  WAS  NCT  ATTEKPTF 
A  TE3FCE A  BY  SWITCHING  INDICATOR  FCF  EATTEFN 
AND  FIDGF  MCVES 

INTERNAL  CCNTFCL  PA  FAME  TEF.  FCP  SUBFCOIINE 
LVAIUE - IF  0  F 1=  0  ,  IFEN  UVALUE  IS  NCT  USED; 

II  J P 3  =  1,  CVALCE  IS  USED  TC  EVALUATE  THE 
CEJECTIVE  IQ  NCI IC  NS  CUFING  THE  EXPLCEATORI 
P  EFT UF.  EA  T  ICN  S 


c 

c 

c 

c 

c 

c 

c 

c 

c 


A ..  A  ? 
211  F 


LN  cm 
3EFS 


t ATPI X  Cl  CEJEC31VE 
1 1 E I  E F. t !»C r S  EEI*  EEN 
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FINCTICK  COEFFICIENTS 
C  EJECT  IV  I  rUNCHC::  7  A! 


fCH  IEE  h~S1  PCI N 3  AND  THE  CCF.FEN2  TEST  ECINT 


INDEX  Cl  V  A  El  A  1 1 1  EI2N0  FE5TUEEEE  AT 
EXELOF.  ATOEY  SEARCH 

DIRECTIC  N  AI  STEF  SI2E  SCALAF-IES - II  POSITIV 

DIRECTION  EXPLORATORY  PCVE,  SZEE=IES;  IF  MCA 
31 VE  DIRECTION  EXF1CFAICEY  KCV  E,  Sf?S=-I?S 


C  ***D ATA  INPUT  GUIDE*** 

C 

C - CAFE  1  FORMAT 

C  liOBJ  15 

C  XPFICE  If 

C  UVAfi  15 

c  :;m  ax  if 

C  IP  HI  NT  15 

C - CA  EC  2 - INITIAL  POINT 

C  X(K)  ,  K  =  1  ,  N  VAR  8E10.0 

C- - CAEC  3 - LCW  ER  BOONES 

C  LBD  (K)  ,  K  =  1  ,11  VAR  8110.0 

C - CAEC  4 - UPPER  BOUSES 

C  U3D(K),  K=1,NVAR  8110.0 

C - CAEC  5 - E IGHT-KAN  E  SICES 

C  EH  I  (I)  ,  I  =1  ,  N  0  B  J  8110.0 

C - CAEC  6 - TOLERANCES  FOB  TERMINATION 

C  SPSY(J)  ,J  =  1,1«PBICF  3110. 0 

C - CAPE  7 

C  EPS - INITIAL  STEF  SIZE  110.0 

C  ACCEL - ACCELERATION  EACTOB  E1C.C 

C  RED  DCS - STEP  SIZE  FED  OCT  I CK 

C  EACTOB  E1C.C 

C  DELTA - S1SE  SIZE  TESI  IACTCB  110.0 

C - CAEC  8 

C  A  MAP  (K  ,  J)  ,  J  =  1  ,  NVAE  FCE  J-TH  110.0 

C  0  EJECT  IV  E  FUNCTION 

C 
C 

C  ***3tBRCtTISES*** 

C  SUERCUTINES  03 J FU N, AC KF UN,  ANC  CVALCE  NEST  BE  CSEF  SUFFLIEC 

C  FOE  EACH  K2R  NIGP  E  EC  El  EE 

C  SUERCUTINE  C3JEUN  CONTAINS  TEE  OBJECTIVE  FUNCTICNS 

C  SUEROUTME  ACHFUN  CONTAINS  THE  ICEIIVEEENT 

C  FUNCTIONS,  ONE  FCE  EACH  PRIORITY  LEVEL 

C  SUEHOUTINE  CVALUE  CONTAINS  THE  CEJICTIVE  FUSCTICN  E  E  EE  ER  ENC  ES 

C 

C 

C  ***HAIN  PROGRAM*** 

C 


RE  AL*  4  IPS,  LED,  N 

IS  1EGER  *2  IMPRCV,IPRINT,ITEEH 

CC  EMC  N/CO  M  MO  1/  NO  B  J,  NP E  10  H ,  N  VA  E  ,  NM A  1 ,  NC  I C L  E 

C3P.MOK/COMK02/X  (2  500)  ,XEASE  (2  500) 

CCCflC N/COMMO  2/A  (  10)  ,  AEEST  (  1C) 

CO  MM 0  N/CO  ;1  S  04/ ACC  EL, REDUCE  ,DE  IT  A , I R  F  ECV  ,  I F  F I NT , IT  I EM 
CO  PMC  K/COMMD  5/LHD  ( 2  50C)  ,UEC(2ECC)  ,SPS,IES 
COKJON/COMPOo/CDJ  (2500)  ,EUS  (2500) 

CC  PMC  N/CO M  MO  7/N  (2  50  0)  ,P  (  2  5CC) 

C3MUC  K/COMEO  c/EPS  Y  (10) 


n  o  r>  n  n  n 
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cc.'.;c::/cc  :.:09/sav  f 

CO.SHOS/COJIJSI  0/JPI  ,LSC£ 
CCf-rXiv/COMKl  1/LIFE  (25CC) 

CO  JS.SCK/CO  SSI  2/ AHA  F  (1  CC,100) 
CC  £3C  K/COiiM  1  3/S  EPS 


BEGIN  TEE  PSOGFAH 


C 

C 

C 

C 

C 

c 

c 


c 


c 


CALL  DA  IAIN 

ALL  NECESSARY  LATA  IS  INPUT  AND  PRINTED  ODT 
CALL  rSJALG 

BEGIN  i  EE  ?  ATT  ERN  SEARCt  ALGORITRS 
ALL  FUR  Hi  E  a  WORK  IS  CCM1ETED  I  b  TEE  SUERC0T1NES 
OE  TEE  CCCE 
SICP 
EMC 


SUEEOGTINS  EAT  AIN 


*444444* 


SUBROUTINE  D&1AIN  FEADS  IV  ILL  TEE  UCESSARY  LATA 

a:;e  feints  it  cut  as  a  xeass  cf  checking  it. 


IS  TEG  E3  *2  I<JFRCV,I?RINT,ITEF1J 
EEAL*4  IPS, LED 

CO  N.IO  N/CO  AZO  1/1-03  J,KEFICF,VVAF,NfcAX,KYCIE 
CCM1CK/COA30  2/X  (2  50  0)  , XEASE(2f  CC) 

CO  x:ic  N/CG  H  JSC  4/ ACC  El ,  FED  UCE ,EE IT  A  ,1 1 E  E  C  V  ,  1  P  F I KT  ,  IT  E  FH 
CCPJSCK/COHJSO  5/LED  (250C),0EL  (2  f  C  C )  ,  EP  £  ,  I E  S 
COJSXON/COSEOb/Caj  (2500)  ,  RHS  (  2  5 CO) 

CCfJSCK/CO ZH08/EPSY  (10) 

CCSXCX/C01U  10/JPI,LNUfS 
CO  JSHO  K/CO  il  151  2/AIi  A  F  (1  CO,  100) 

READ  ALL  TEE  NECESSARY  LATA 
READ  170,NCBU,KP310F,FVAS,N!SAX,IFFIM 
REID  175,  (X  (X)  ,K=  1,  NVAB) 

READ  175,  (IBD  JK)  ,K=1  ,SVAR) 

READ  175,  (UBC(K)  ,K=1,NVAH) 

READ  175,  (RH  S  (I)  ,1  =  1  ,KCEJ) 

READ  175,  (EPSY  (J)  ,J=  1,NFB3CE) 

HEAD  175, EPS, ACCEL, REDUCE, LILIA 
DC  100  J= 1,NCEJ 

100  READ  175,  (ASA  F  |L,  J)  ,1=1  ,KVA&) 

PRINT  AN  ECEO  CHECK  OF  THE  DATA 
105  PRINT  160 ,K03J,KP RICF,  KVA E,NK AX, EES, ACC  II, FETUCE, CELT A 
FEINT  181 
DO  110  K=1,NVAR 

110  PRINT  182, K,X  (K)  ,  R,LEE(K)  ,K,UED(K) 

PR  INI  183 
DC  120  1= 1,NCEJ 
120  PRINT  1  84  ,1 ,  RH  S  (I ) 

PRINT  185 
DO  1 3 C  J=1,SPRIOR 
130  PRINT  136,  J,  EPSY  (J) 

PR  IN  I  1&7 

IF  JIP  FI  NT  .EC.  1)  GO  TO  14C 
PRINT  188 
RETU  Eli 

170  FOSSA  3  (1615) 

140  PRINT  139 


u 


•  I 


175  FJhlii  1  (cr  1  C.  0) 

130  FORMAT  (*-•  ,25X,*  NUMBER  CF  OBJECTIVES  =  •  ,T57, 15,/ ,  •  * ,  2 5X,  JM  HEP.  C 
3F  EIICBITIES  =  *,i'57,  35,/,  •  1 , 2  5X  ,  *  N  L  KB  EF  CF  DECISION  VARIABLES  =  *, 
2157,15,/,*  •  ,25X,  •  MAXIMUM  NGEEE5  Cl  All  C»  AELE*  ,/,  *  •  ,  29X,  •  ? ATI  El- K 
SSEArfCE  CYCLES  =•  ,  T57,  1 5, /,  *  *  ,  25X, *  INI T IAL  STEP  SIZES  =',T59,r10. 
S  6,  /,  *  •  ,25X,  *ACCSIZRA1ICN  r  ACTCB  =  •  ,T59 , 1 10.  6,/,  •  ,,25X,*SIEP  SIZE 
i  F  5  EUC7  ICS  SACTCK  =' ,  T59,  f  1  0.  6  , /,  *  *  ,  25X  ,  «  MINIMUM  ALLCViAELS  STEP  SI 
sz£  =•  ,r5S,F10.€) 

181  FORMAT  {•-  1  ,251,*  V2CTC5  C?  INITIAL  cSI  1 1.  IT  FS  FC  5  DECISION  7  A3  I A  EL  ZS 
'  ,  25X,  '  LCi  EE  AND  OFF  EE  BCIMCS  ON  DECISION  VARIABLES - IE 

) 

,10X,*X  (*  ,  15,*  )  =  *,F10.5,7X,*13D(*,I5,  •)  =  *,F1C.5, 

15,')  =  • ,  1 1 C.  5) 

,25X,' RIGHT-HAND  SICE  VAIUE3  ICE  CEJ5CTIVE  FUNCTIONS - P. 


i— x»,/.» 

3D  AND  U3D 
132  rCFMATJ* 

5  7  X  ,  *  U  ED  (• 

183  FORMAT  {•- 

3  as  • ) 

184  FOR.1A1  {* 

185  FORMAT ( •  — 

186  FORMAT  (» 

187  FORMAT  (*- 
S’) 

138  FORMAT  (' 


,25X,'3HS  { •  ,15  ,  *)  =  '  ,  E 1 3.  4  ) 

,  25X,  '  V  ECICE  CE  ACHIEVEMENT  FtNCTICN  TCIZFA.NCES - EESY  •) 

,25X,  *E?SY  (  *,I2,  *)  =  *,F1C.t) 

,2 5X ,•  All  DATA  HAS  SEEN  INEUT EEGIS  THE  PATTERN  SEARCH 

,2  5X , ' NCT  E :  IF  S I  NT  =  0,  SC  C  MY  THE  FINAL  SOLUTION  WILL 


33E  FEINTED  COT*) 

189  FO  PM  A  T  (*  *  ,2  5  X  NCT2 :  IFFINT  =  1  ,  SC  1H  SOLUTION  WILL  EE  PRINTED 
3  CUT  AT  EVERY  PATTERN  SEARCE  CYCLE*) 

PS  TORS 


END 


<444444444444444+44; 


C 

C 

C 

C 

C 

c 

c 

c 


SO  ERGOT  IN  E  HJAIG 


4444444 


SUBROUTINE  HJALG  IS  THE  KCCKE- JEEVES  E  ATT  EPS 
SEAFCH  AL  GO SIT  EM.  IT  PERFORMS  EXFLCFA1CR  Y  AND 

PATTERN  MOVES  TO  FIND  THE  'EEST  '  SCIUTICN  TC  THE 
NONLINEAR  GCAL  PROGRAMMING  PROBLEM. 


C 

c 

c 

c 


100 

c 

c 

200 

201 


RE AL*4  IP  S  , LBD 

INTZGER*2  IMF5CV,  IPS  3  NT,  I1EF.K 
COMKOK/COMMOI/SCB J, NFFICS, EVAE,FtAX,  NCYCIE 
COMMON/ COM  MO  2/X  ( 2  500)  ,  XE  AS  E  ( 2  5  C C) 

CO "MO  K/CO  2 MO  3/A  (1 0)  , AEEST  (10) 

COKMCK/COMK0  4/ ACCEL,  E  1 1 UC E  ,  D El TA  ,  I MPE C  V  ,IFRI  NT  ,1  TERM 

CO  MHO  N/CO  M  MO  5/ 1 3D  J25C0)  ,CED  (2500)  ,IFS,IIS 

COMMON'/ COSMO 9/SAV  E 

COMMON/COME1  0/JFI  ,LNUS 

CC  CMC  S/COMM  1 3/S EPS 

J?I  =  0 

CALL  ACHFUS 

SUBROUTINE  ACHFON  FCRCS  THE  ACHIEVEMNT  VECTOR  FOR 
THE  INITIAL  BASE  POINT.  IT  IS  TESTED  FOR  AIL  ZEFC 

VALUES - 1?  SO,  THE  INITIAI  ECIM  IS  CETIMAL. 

OTHERWISE,  IT  IS  SET  AS  IRE  BEST  VALLE  AND  IS  FFINTED. 
DO  IOC  J=1  ,NF3I08 
IF  (A  (J)  .  N E.  0.0)  GO  TO  2C3 
CONTI  SUE 

ALL  FRICRITY  LEVELS  ABE  SATISFIED  AT  ZERO  VALUE. 

THE  INITIAL  POINT  IS  C FTI  E A I---I E F C 1 » AT F. 

PRINT  290 

DO  201  J*1,SPSICH 

PRINT  282,  J,  A  ( J) 


ononn  r>  nnnnoo  nn  n  n  n  r>  r>  r> 
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prini  291 
DC  20  2  X=1,NYAE 

202  PRIST  2S0  ,  N,  X  IS) 

GC  TO  295 

203  PRINT  281 
CC  20  4  0=  1,N PRIOR 
A3  ESI  (0)  =  A  |0) 

204  PRINT  232,0,3(0) 

**  INI1IAII2ATICH  C?  PATTERN  SHFCE  PAR/MFTEnS  ** 

NCYCI2  COUNTS  THE  S 0HEEB  OF  PATTERN  SEARCHES  DEED 
IACC5L  IS  AN  INDICATOR  fcHICH  TEILS  II  TEE  EXPLORATORY 
MOVES  ABE  MADE  AECUT  AS  ACC  EL  EF.  A  TICS  POINT 
SCYCLE  =  0 
I ACC El  =  0 

SET  THE  I  NX  HA  I  BASE  ICIM  AND  1BI  1MTIAL  STEP 
SIZE  VECTOR  BEFORE  EXPLORATORY  MOVES 
DO  2C£  K= 1  ,N VAR 

205  XBASE  (X)  =  X  (K) 

IPS  =  EPS 

START  TEE  PATTERN  SIASCE 
**  EXPLORATORY  MOVES  ** 

206  IF  (NCYCLE  .GT.  NS  AX)  GO  TO  26C 
DPI  =  1 
IP  FOG  =  0 

IPRCG  CCCNTS  THE  KU2EEE  CF  IEERCVFMMS  CURING  TEF 
EXILC2ATCEY  MOVES.  IT  IS  USED  TO  DECIDE  IF  A  EAT1EES 
MOVE  SHOULD  3£  FEEFCEEEE. 

207  DC  210  K=  1 , N V AB 

PERTUR3ATICSS  IN  THE  PCSITIVE  IJFFCTJCK 
TEST  STEP  SIZE  COMPONENT.  IF  LESS  TRAN  THE  MININGS 
ALLOViA  3  IE  STEP  SI2E,  A  FItGE  SEAFCE  V  HI  EE  PER  FORM  EE 
IF  (IPS  .LT.  EELTA)  GC  TO  2ES 
LNGM  =  K 

X (K)  =  X  (K)  +  IPS 

SEP  S  =  IPS 

IF  (X  ( X)  .  LT.  LED { X) )  GC  TO  2C6 
I?  (X  ( K)  .Gl.  G3D  (K)  )  GC  TC  208 
CALL  ACHFUN 

ACHFUN  FORMS  THE  ACHIEYEtE  NT  VFCTCF  AT  THE  TEST  FCIHT 
ITFRH  =  0 

I?  (K  .EQ.  NVAfi)  ITEBE  =  1 
C AIL  EFCIDE 

DECIDE  COMPARES  THE  ACHIEVE2EKT  VECTCFS  CF  THF 
BASE  POINT  AND  THE  TEST  POINT.  IF  X  (X)  IMPROVES  THE 
SOLUTION,  IT  BECCEES  A  TEKECFAFY  FIJI  F  Cl  NT  ANE 
EX  PLCEATOfiY  MOVES  ABE  CONTINUED  FRCE  IT.  IF  THERE 
IS  NC  I MPBCVEHEHI,  THEN  EXFICFF  IF  TFI  NEGATIVE  EIFICTION 
IF  (IBPROV  -EQ.  0)  GO  TO  2CE 
IPROG  *  I  PROG  ♦  1 
GC  TC  210 

C  PERTUB  B ATI CNS  IN  THE  NEGATIVE  DIRECT  1CK 

208  X  ( X)  =  X  (K)  -  2.  0  *IPS 
SEPS  =  -IPS 

IF  IX  (K)  .  LT.  LcD(X)  )  GO  10  209 
I?  (X  ( X)  .GT.  G3D  { K)  )  GC  TC  209 
CALL  ACHF'JS 

C  ACHFUN  Feans  THE  ACHI EVEEEM  VECTCF  /  T  THF  TESI  PCINT 

irEBH  *  0 

I?  (K  .EQ.  5VAB)  I TEBE  *  1 
CALL  EFCIDE 


non  n  non  nnnn  n  n  n  n  n  n o n  n  n  n  n 
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II  X  (X )  IMPROVES  THE  SCIUTICN,  IT  EICCMS  TEE 
TrMChAfY  fcrAE  POINT  AM  TEE  Nr  XT  V/FIAZLF  IS 
PERTURBED.  If  THERE  IS  KC  IMPRC  YEMENI ,  1c  ASE(K) 
REMAINS  THE  TELPCESRY  hEAI  ICIM. 

IF  (IKIFCV  .  E C •  0)  GO  1C  2CS 
IP50G  =  I PRC3  ♦  1 
GC  TC  213 

209  JC  ( K>  =  X  ( K)  ♦  IPS 

210  CC  NT  I  NO  E 

NCYCLE  =  NCYCLE  ♦  1 

THE  EXPLORATORY  MCVE  DC  ICC?  HAS  CCS  i  IE  LCTED  A 
SEN  POINT.  TEST  THE  ACHIEVEMENT  VECTCR 
FCf  ALL  Z  EEC  V  ALOES 
DO  211  J=1,NFRIC8 
IF  (A  ( J)  .  S  E.  0.000)  GO  TO  212 

211  CONTI  NOE 
GC  TC  200 

THE  ACHIEVE  KENT  VECTCE  IS  NCT  All  Z  F  F  C  5S 
TEST  THE  POINT  FOR  IMPROVEMENT-  IF  IT  GIVES 
IMPROVEMENT,  A  PATTERS  MCVE  IS  EEFfCIMEC.  I?  NCT 
THEN  STEP  SIZES  S ILL  EE  REDUCED  ASD  NEfc  EXPLCFATCRY 
MOVES  Sill  Si  FERJCRMEE 

212  IF  (IPFCG  .  GT.  0)  GC  TC  214 
I?  (IACCEL)  250,213  ,215 

THE  SXFLCRATJ3Y  MOVE  HAS  PEES  CSS  tCC  E  SSf  01 .  STEF 
SIZES  -ILL  SE  REDUCED  AND  NEW  EXFICI  FTC  FY  MOVES 
-ILL  BEGIN. 

213  PRINT  283 

**  REDUCE  STEP  SIZES  ** 

IPS  *  IPS  *  SEDUCE 
GC  TO  206 

THE  EXPLORATORY  MOVE  HAS  EEEN  SDCCESSEU1. 

A  EATTEFN  ilCV  E  »  ILL  BE  PEE  FOBS  EE 

214  PRIST  2  79 

**  ACCELERATION  *♦ 

JPI  =  0 

215  DC  220  K=  1 ,  N V  A R 
SAVE  =  X  ( K) 

X  ( K)  a  X(K)  ♦  ACCEL  *  (X(K)  -  X3ASE(K>) 

I?  (X  (  K)  .LT.IBD(K)  )  X{?)  *  IED(F) 

IFJX(K)  .LE.  U  ED (K) ) GO  TO  220 
X  (K)  =  UBD  (K) 

2  20  X3ASE  (K)  =  SAVE 

TEST  THE  AC  KIEV  EM  ENT  VECTOR  AT  THE  ACCELERATION  FCIS1 
CALL  ACHFUK 
ITEBM  =  1 
CALL  DECIDE 

IF  (IK i FCV  .EQ.  1)  GO  TO  225 
IACCEI  =  -1 
GC  TO  230 
225  IACCEI  =  1 
GC  TC  240 

THE  ACCELERATION  FT.  GIVES  NC  IEEFCVIEI  NT. 

EXPLORATORY  MOVES  i  ILL  EE  FEHECRMEC  AEOCT  THE 
ACCELERATION  PT.  TC  FINE  IMPROVEMENT. 

230  PRINT  286 

231  I?  (IP  FI  NT  .EQ.O)  GO  TC  206 
PRIST  267, NCYCLE 

DO  232  K= 1  ,N  VAR 
2  32  PRINT  230  ,  K,  X  E  AS  E  ( K ) 


noon  n  n  n  n  n  onn  no 
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?:<  i:;t  23a 
rc  23  5  j=  i,n price 
235  PRIST  2  82  ,  J,  A6EST  (J) 

GC  TC  206 

THE  AOCELERATI  CN  FCI1.1  GIVIS  IEEFCVIHVT. 

EXPLCEE  ABOUT  IT  fCB  IZ  ERC  V  EZ  EN  T  . 

240  PRINT  285 

IF  1 1 E  F I  NT  .EC.  0)  GO  TO  2C6 
PRINT  2  87  ,  KC  YCLE 
CO  24  1  K=1,NVAE 

241  PRINT  290,  K,  X  (K) 

PRINT  288 

DO  245  J  =  1  ,NP  RICH 
245  E8INT  232,  J,  A  ( J) 

GO  TO  2  06 

EXPLORATION  AEOOHE  IKS  ACC EL ERA lICS  POINT 
DOES  NOT  HELD  ANY  lEEFCVEfEKT 

RESET  TEE  EASE  POINT  ANE  GC  TC  EXPLCFA1CRY  SOVES 

250  I?  IIACCEL  .EQ.  0)  GC  TC  255 

251  DC  25  2  K=  1 , N V  A B 

252  X  (K)  -  XB  A  SE  (K) 

I ACCEL  =  0 

GO  TO  206 
255  PRINT  284 

259  CALL  RIDGE 

SUE  SOU  TINE  RIDGE  KILL  TRY  TO  FINE  A  NEK  BASE  POINT 
IN  AN  OHLI v 02  DIRECTION  TC  THE  IFESEM  EX  FLCSATCSY 
SEARCH  CIEECTICN.  II  SUCCESSFUL,  THE  PATTERN 
SEARCH  WILL  BEGIN  AT  THIS  NEK  ECINT.  IE  UNSUCCESSFUL, 

THE  ALGCRITEM  WILL  EE  TERMINATED. 

IACCEl  *  0 
GC  TC  206 

THE  “AXIflUE  N  0  ME  Z  E  OF  EATTIRN  SIAFCF  CYCLES  HAS  EErS 
EXCEEDEE.  THE  BEST  VALUES  OF  THE  YAEIAELES  AND  THE 
ACHIEVE  RE  KT  VECTCF  AFE  PFINTEE  CUT,  AND  TEE  ALGCP.IT  BN 
Will  BE  TERMINATED. 

260  PRINT  289 

IF  JIACCEL  .LT.  1)  GO  TC  265 
DO  261  K= 1  ,N  VA3 

261  PRINT  230  ,  K,  X  (K) 

PRINT  288 

EC  262  J=  1 , N P R IOB 

262  PRINT  282  ,  J, A  (J) 

GC  TO  295 

265  DO  266  K=1 ,N VAR 

266  PRINT  280  ,  K,  XEASS(K) 

PRINT  288 

DC  267  J=  1 , N P R IC B 

267  PRINT  282  ,0,ABEST  (J) 

270  GC  TC  295 

279  FORMA  T  | '  ,2  5  X  ,  • THE  EXELCFATCFY  ECVE  HAS  EIFN  SUCCESSFUL' ,/, •  ',25 
IX, 'A  EATTEFN  NOV E  WILL  EE  FEE FCK3ED ' ) 

280  FORMAT  (•  •  ,34X,  *X  ( '  ,15  ,  •)  =  ',113.6) 

281  FCFHATC  •  ,  2  5X, '  INI!  1 AL  AC  El  E  V  E  XE  NT  VECTCF  VALDES') 

232  FCFHATp  •  ,  24X,  •  A  ( ',  12,  ')  =  »,F13.6) 

283  FOBNA1I*-'  ,25X,'  TEMFCFABX  HEAE  PCI  IT  ECIS  NOT  IMFFCVE  IHE  SOLUTION 

3',/,'  ' , 25X, ' IEE  STEF  SIZES  WILL  EE  RED  ICED',/,'  *,25X,'fEW  EXI1CF 
3ATCBY  MOVHS  WILL  EE  EEFFC5EED  AECOT  XEASI') 

284  FCFZAT  {'  '  ,25X,'  2  AXIMUM  STEP  SIZE  FEDLC1ICNS  EXCEEDED  FCE  THIS  EXP 
SLORATICN' ,/,  •  ',  2  5X  ,  'A  FIEGE  SIARCt  ULI  BE  FERF  CE  FED'  ) 


2o5  f:r::aic-  »  ,25x,*the  accziefaticn  fcim  ii-iES  ai;  improve:  solution 

5'#//'  ',25X,  »  EXPLOSAICRY  MOVES  fcILI  BEGIN  AT  THIS  NE»  FOIST*) 

266  FORMAL  ( '-  •  , 2 5X  ,  T  HE  ACCZtEFATICK  POINT  ICES  NOT  ’Itll  A!  IMRCVES 
i SOLUTION*  , /,'  • ,25X,' EXPICFATCFY  MOVES  »UI  EEGIK  AT  TEIS  POINT') 
237  FCPMAI  ( •- •  ,  2  5X  ,  '  I  HE  E  ESI  SOLUTION  AT  PATTERN  SZAECH  NUSEES  ',13,' 
IFOLLO  kS  •) 

2  36  FC  FM AI  (*  *,25X,'TfcE  AC  E1EV  EXEN  1  VECTOR  VALUES  ARE:*) 

269  FORMA  1  ( '  ' , 2 5 X  ,  • M AXI E 0 E  NUtEEF  CE  FATIFIN  SEARCH  CXCIES  E  XC  EE  C  ED '  , 
*/,»  •  ,2  5X  ,  '  T  H  £  ALGORITHM  IS  T  ESMIN  AT  E  E  • ,/,  •  ',25X,'THE  E  E  ST  SCLLtTI 
5  CN  rc  THIS  PCINT  FOLLCSiS') 

290  FORMAT  ,20X,43  (IE*)  ,/,  •  *  ,25X,  '  CFT  IS  AL  ACHIEVEMENT  VECTOR  VALUE 


290  FORMAT  (•-  •  ,2  OX, 43  (1 H  *)  ' 

SS',/,'  • ,  20X  ,  43  (1  E* )  ) 

231  FORMAT  (•-  •  ,  20X,  43  ( 1H  *)  ,/,  • 
'  ,2  OX,  43  |1H  *)  ) 

295  SIC? 

END 

*************  ***************** 

SU230  LTIN  E  DEWAR 


*, 25X, 'OPTIMAI  DECISION  VARIABLE  VAIUES 


SUERCUTINE  EZVVAR  CALCOIATES  THE  POSITIVE  AND 
NEGATIVE  DEVIATION  VA5IAEIES  ECF  EACE  CEJECTIVE 
FONCriOK  USING  THE  PREVIOUSLY  EVAICATED  CEJECTIVE 
FUNCTIONS  AND  THE  RIGHT-HAM  SHE  VAIUES. 


REAL* 4  N 

CO  MMON/COMMOI/NCBJ,  NFBICF,1VAF,SEAX,  NCYCIE 
CCEMCN/COMMO  2/X  (2  500)  ,XEASE  (25CC) 

CO  MMO  N/CC  HM06/C3 J  (2500)  ,SHS  (2500) 

CC  EMC  N/COMM0  7/N  (2500)  ,P(  25C0) 

CO  KHO  N/CG  1  Ml  0/JFI  ,  LNCE 
IF  (JFI.  5Q.  0)  GO  10  4CC 
CALL  LVALUE 
GC  TC  405 
400  CALL  CBJFUK 
405  DC  430  1=  1 , NO E J 

IF  (03J  (I)  -  RHS  (I))  410  ,420,420 

4  10  N  (I)  =  BHS  (I)  -  OEJ  (I) 

P(I)  =0.0 
GC  TO  430 
420  N  (I)  =  0.  00 

?  (I)  =  OBJ  (I)  -  RHS  (I) 

4  30  CONTINUE 
RETURN 


SUBROUTINE  DECIDE 


*********** 


SUEBCUTINE  DECIDE  IS  USED  FOB  HC  PCFPCSES: 

(1)  TO  COMPARE  THE  ACHIEVEMENT  VECTCE  VALDES  CE 
THE  EASE  PCINT  AND  A  TEST  PCINT  DURING  THE  EXPICFATCFI 
MOVES  AND  SELECT  THE  EEST  E  Cl  NT  ECF  IUTUEE  MOVES 

(2)  TO  TEST  AN  ACHIEVEMENT  VECTCE  FCE  1ZFEI NATION 
0?  THE  AL3CFITHB.  THIS  CAN  CCCUE  IN  TWC  METHODS: 

**** TERMINATION  EY  MEANS  CF  ZEEC  FCR  ALL 
PRICRIII  LEVELS,  AND 

♦•♦♦TERMINATION  EY  MEETING  SPEC1FIFD  TCLEE A NCES 
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:  =  5  pa:;a..e:er  iieed  is  ::;dic;ic  =  :-c?  isreilatics 
testing.  i?  iiife  =  c,  c:;iy  ike  ccrFAii sc:;  test  is 

?Zi  ?  C=  .ISC  .  IF  I  1 E  5  £  =  1,  SHE  TiFfiK/TICS  TESI3  WILL 
ALSC  BE  ?  EF.  FOR  WEE. 


IN  TEG  EH  *2  ISFaCV,IP3IM,ITEEE 

CCMC N/CGH20  1/ NOB J ,  SF5ICE,  KVA F. ,  N HA  S  ,  SC  1 CIE 

CO  FHCN/COHH02/X  (2  500)  ,XEASE  (2500) 

CCFECL/C0:'a5!0  3/A(  10)  ,  AE  EST  (1C) 

CO  5110  K/CO  3 1!  C 4/ACC  EL,  FE  BOCE  ,  E  Ell  1 ,1 2F5  CV  ,  1 F  PINT  ,  IT  E  FS 
C  C  2  H  C 1/  CC  J  KO  6/  C  5  J  (2 50 C )  ,  RES  (25CC) 

CG  WHO  K/CO  B  SO  3/EPS Y (  10) 

cc  ekc  K/co w a i  C/ o?i ,  i  s  c  c 
CCKHCK/C3HW1  1/Elf  F(  25CC) 

**  CCWPARISCS  CF  ACHIEVE2EKT  VECTCFS  FCF  EX  FLCF.  AICRY  HCVES  ** 
DC  610  J= 1,NFF10H 

I?  (A  (J)  .ST.  A32S1  (J)  +  0. 00001)  GC  1C  63C 
IF  |A  (J) +0 . 0000  1  .  LT.  AEEST(J))  GO  1C  62C 
CONTINUE 
SC  TC  6  30 

IHE  12  S 1  PCI  NT  I EEFCYES  THE  SCIUT3CK 
IHIRCV  =  1 

CHECK  IHE  INDICATOR  FCF  TEF HI K STICK  TESTING 
IF  (IT  FEM  .  GE.  1)  GO  TO  65C 
GO  TO  670 

THE  TEST  POINT  DOES  NOT  IWPEO  VE  THE  iCLCTICN 
I  SEROV  =  0 

IF  (J  F  I.  EQ-  0)  GC  TO  6<iC 

DO  63  5  J=1  ,NC3J 

CSvI(J)  =  0  E  J  (J)  -  C  IFF  (J) 

I?  (ITEHH  .EQ.  2)  GO  TC  650 
F  E IU  S  K 

**  1EE  ill  NATION  TESTS  ♦* 

{1)7  EF.  H  IN  AT  ION  TEST  ET  WEANS  OF  ZEBC  FCE  ALL  FEICBITIES 
DO  655  J=1 ,N? RIOB 
IF  (A  (J)  -  HE.  0.0)  GO  TC  66C 
CONTINUE 

ALL  FP.1CRIT7  LEVELS  ABE  ZEFC - THE  CE1I2AI  SOIITICK 

HAS  3E2N  FCUND - TERPIKATE  THE  71GCF3TEE 

GC  TC  680 

(2)  TEB  HI  KA  1ICS  TEST  El  TCIIFAKCES 
IF  (II  EFH  .  EQ.  1)  GO  TO  67C 
DO  665  J=1  ,Nrr.ICB 

IF  (AES  (A  ( J)  -  AEEST  (J)  )  .GT.  EPSY(J))GC  TC  670 
COSH  NOE 

ALL  TOLEEANCES  ABE  SAI3SEIED - A  SCllTICN  HAS  BEES 

FCCND - THE  SEARCE  ALGC  FITK  £  IS  CCEE1ETE 

GC  TC  630 

A  TOLERANCE  IS  SC T  SATISFIII - 1  EIT1IF  SCLUTICN 

WILL  BE  SOUGHT - FETUSS  TC  THE  PA  IT  EF  K  SEARCH 

IF  (IHFHOV  .EC.  0)  GC  TC  695 


A  TOLERANCE  IS 
WILL  BE  SCUGHT- 
IF  (IHFHOV  .EQ .  0) 
DC  675  J=  1,  NFF.IOB 
AS  EST  (J)  =  A  (J) 
FETUBK 
PRINT  690 
PRINT  691 

DO  682  J=1 , N  FHI CB 
FEINT  692,  J,  A  (J) 
PRINT  653 
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u  3  5  1  ,  N  '»  n  P 

C3-:  HIM  £54  ,  ?,  X  (F) 

SIC? 

69C  F  0  R.'.  A  1  ( •-  •  ,2  5  X  ,  •  THE  SEARCH  A1GCF1LHM  IS  CCMFLETI - ALL  TOLERANCES 

BFc?  a  solution  AK  1  SiHSnEC) 

o  3 1  FORMAT  (•-•  ,20X, 43  (1H*)  ,/,*  •  ,25 X ,  •  C P T  I  ? 7 L  ACHIEVEMENT  VrCIGS  VALUE 
IS*,/,*  •,  20X,  42  (lb*)  ) 

6  92  FORMAT (•  •  ,  3  <4  X  ,  •  A  ( *  ,  I  2  ,  • }  =  »,F  15-€) 

653  FORMAT  (•-*  ,20X,43  (1  H*)  ,/,•  •  ,  25X  ,  •  C  FT  I  S.  71  EPCISICR  7  AFIAELZ  VALUES 
S*,/,*  «,20X,  42  (  1H*)  ) 

654  FORMAT  { •  •  ,34X,»  X  (•  ,I5,«)  =  ■ ,  r  1  3-  6 ) 

635  RETUEK 

END 


SUESC  CIINE  EIDGE 

C  i<:tiiti****i**U^t«*»***H^t*^**:*n*»**4t*nt*'*4*^*?4UimU;U*4,i** 

C 

c 

C  SUBROUTINE  RIDGE  IS  CAIIIE  WHEN  TEE  EAITFEN  SEARCH 

C  CAN  NO  LONGER  FIS  t  PA  IT  EF.  N  CR  EXPLCRATCPY  MOVES 

C  WHICH  I  KPS CVS  1HE  ACHIEVEMENT  VECICI.  FIEGE 

C  EVALUATES  EXPL  CB  AT  CRY  POINTS  WHICH  AFE  IS  OBLIQUE 

C  DIRECTIONS  TO  THE  USUAL  EXEICfATICN  AXES.  EIIGI 

C  ATTEMPTS  TO  FINE  A  EESCLUTICS  RIDGE,  IF  IT  EXISTS, 

C  AND  .10  VE  THE  PATTERN  IK  THAT  EIFECTICN. 

C 

c 

DI  M  NSICS  APCS  ( 10)  ,  AN  EG(  1C) 

REAL *4  I P  S  ,L3D 

INTEGERS  IK P RCV ,  IPE  IN T,  I T ERR 

CO  XMO  R/CO  M  MO  1  /  SC3  J  ,  N I  El  C  F  ,  N  V  A  F  ,  R  !  A  X  ,  N  C  Y  C I  E 

CC  *  NCR/ COM  MO  2/X  ( 2  50 0}  ,  XE  AS  E  (  2  5  C  C) 

CO  MHO  N/CC  K  SO  3/A  (1  0)  ,  AEEST  (10) 

CCMMCR/CO;iMOVACCEL,REUDCE,DELTA,IMPSCV  .IPEIST  ,ITEEM 
CO  KMO  K/COMM05/LBD  (2500)  ,UEE  (2500)  , IIS, IIS 
CCMMCN/COMM09/SAV  E 
CO  KMC  R/CO  M  Ml  0/ J  FI  ,  L  N  0  M 
CC  MMCN/CCMM1  3/SZPS 

C  PE8TCS5  EACH  VARIAEIE  CNE  AT  A  TIME 

C  PEFTUSE  IN  THE  POSITIVE  DIRECTION 

DO  80  C  K= 1  ,N  VAE 

800  X  (X)  =  XBASE  (K) 

JPI  =  1 

801  DC  865  X=  1 , S V  AF 

X  (  K)  =  X  (  K )  ♦  IPS 
IF  |X  CK)  •  LT.  LED(K))  GO  TO  6C5 

IF  (X  (  K)  -GT.  U3D  ( K)  )  GC  TC  810 

GC  TC  815 
805  X  ( K)  =  L3D  (K) 

GC  TC  8  15 
810  X  (K)  =  U3D  (K) 

8  15  CALL  ACHFUN 

C  AC  nr  UN  FORKS  THE  ACHIEVEMENT  VECTCE  IN  THE  POSITIVE 

C  DIFECTICN 

DO  b2C  J= 1 ,NPRIOS 
820  APCS  (J)  =  A  ( J) 

C  PEE1UBB  IR  THE  NEGATIVE  EIFECTICN 

X  (K)  =  X  ( K)  -  2.0  *  IPS 
I?(X(K)  .LT.  L2D  ( K)  )  GC  TC  825 
15  (X  ( K)  .  GT.  0  ED  (K )  )  GC  TO  830 


no  n  nnnon  nonn  nn  non  no 
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S25 

830 

335 


840 


350 

355 

8  6  0 


8  51 


865 


8  66 


367 


868 


369 

870 


j  r  c  3  3 

Z(K)  =  1?  j  JK. ) 

GC  TC  8  35 
X  l  K)  =  'J3D  IK) 

CALL  ACHF'JN 

AC  H? 'JX  FORMS  THE  ACHIEVEtEFT  VFC1CF  11  CHE  NEGATIVE 
DIF  ECT  ICS 
DO  8  4  C  J=1  , lir SI  CR 
AN  EG (J)  =  A  ( J) 

COMPARE  A  PCS  ( J )  AFC  A5EGJJ).  SEIFC1  TEE  MNIHUK 
VIA  A  STEP*  IS  E  PS  E  EKPTIVE  Pr.OCSElRE.  TRIE  DETEF2INES 
A  DIRECTION  FCt  THE  SIEGE  ECIM. 

EC  850  J= 1 ,  N  PRIOR 

IFJAPCSJJ)  .IT.  ANEGJ0))GC  TC  855 
IF  (A PCS  (J )  .  GT.  A NFG  { J)  )  GC  TO  EEC 
1?  (J  .EQ.  l.PSICR)  GC  TC  855 
CC.MIMJE 

X  ( K)  =  XBASE  JK)  ♦  IPS 
GC  TC  861 


X  (K)  =  XB  A  EE  (K)  -  I  PE 

ciETOFN  X  { K)  TO  ITS  EASE  POINT  VALLE,  XBASE,  BEFCRE 
THE  ?E  R 1US3AII CN  CF  THE  NEXT  VAFIiElI 
S  AVE  =  XBASE  (K) 
aSASE  JK)  =  X  (K) 

X  (K)  =  SAV  I 

THE  NEK  PCIKT,  XSIEGE,  HAS  EE  E  F  ECFCIE  IN  AN  CEIIQ'JE 
Dir  ECT I CK  TC  TEE  GL £  EXELCF ATCE I  HCVES.  ITS 
ACHIEVEMENT  VECTOR  IS  CCKFU1EI  ill  CCSEAF.E £  TC  THE 
3EST  VALUE  AT  POINT  XEAEE. 

JPI  =  0 

DC  36  6  K=1,KVAE 
SAVE  =  XBASE  JK) 

>:  =  AS  £  JK)  =  X  (K) 

X(K)  =  SAVE 
CALL  ACHF'JN 
ITZRH  =  0 
CALL  IFCICE 

IF  THE  RIDGE  PCI  NT  GIVES  KC  I?.  FFCVEEfFT,  THE  PRESENT 
BASE  POINT  IS  THE  EEST  SOLUTION,  AND  THE  ALGGIITH  E 
IS  IER  .11  GATED. 

IF  THE  SIEGE  POINT  GIVES  IEPROVEflENT  ,  AN  EVEN  3ETIEE 
POINT  WILL  EE  SCUGHI  ALCFG  THE  FESCIOTICN  RIDGE. 
IFJISPFCV  .Eg.  1)  GO  TO  £65 
CONTINUE 

X  E  AS  £  IS  THE  EEST  SOLOTION 
P3 1ST  £91 
PRINT  897 
DO  867  K=1,FVAB 
FEINT  392,  K,  X  EASE  (K) 

PRINT  698 

DC  868  J= 1,NP5IOB 

PRINT  693  ,  J,A3EST  JJ) 

ST  CP 

THE  RIDGE  PCINT  GIVES  A  F  IFFSCVEC  SCIOTICN. 

TEI  TO  ACCELERATE  ALONG  THE  RESCLITICK  PIDGE 
PRINT  894 

□ACCEL  *  2.0  *  ACCEL 
DO  871  K=  1  ,N VAR 
SAVE  =  X ( K) 

X(K)  =  X  ( K )  ♦  RACCEL  ♦  (XJK)  - 


XEASE  |  K  )  ) 


non  r»  n  o  o 
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c71  XE  #SE  {.<)  =  SAVE 

CALL  ACKFUN 
llEF.il  =  0 
CALL  CECICS 

I?  IHE  AC  CZ  LES  All  C  K  FCIM  JIEIES  If  It  l  rCVED  SCLUIICS, 

T  hi  TO  ACCELERATE  EARLIER  ALONG  IHE  RESOLUTION  PIDGE 
IF  THERE  IS  NO  IEFFCVEEENT  ,  THE  F1ICI  FCIM  EECC/.ES 
THE  SEN  POINT  EOE  1F.E  1  EN'E  *  ED  PATTERN  SEA  PC  H. 

I?  (IK  E3QV  .EQ.  0)  6 C  1C  872 
PRINT  896 
GO  TO  £70 

872  FRIST  895 
PRINT  697 

DC  87  3  R=  1,:iVAE 

873  PRINT  £92  ,  F,  XBASE  (X) 

PRINT  898 

DO  875  J=  1  ,  S  P SIDE 
975  PRINT  393, J, SEES! (J) 

RESET  THE  BASE  FCI  ST  ARE  STEF  SIZE  VICTCFS  AXE 
RETURN  TO  THE  EX  PIC  5 A  TCFY  MOVES. 

NO  IMPROVEMENT  FPCE  ACCE IE  E  ATI C  F  PC1FT 
DC  877  K=  1 ,  N  V  A  fi 
977  X  '  K)  =  XBASE  (K) 

IPS  =  EPS 
DO  87  £  J=1,NFBI0R 
879  A  (J)  =  ABEST(J) 

GO  TO  699 

391  FCFHAT  ('- •  ,25X,'TEE  RIDGE  SEARCH  IS  CC.MPIEIE',/,'  ',25X,'NC  IMFCV 
i  ED  SOLUTION  CAN  BE  FOUND*,/,'  ',2SX,'TEI  CURFENT  E  AS  F  FCINT  IS  THE 
5  BEST  SOLUTION  THAT  CAS  BE  FOUND  THROUGH  IKE  PATTERN  SEARCH') 

392  FCFHATf*  •  ,34X,  *X  { ',  35,  •)  =  •  ,  F  1 3 .  €) 

893  FORMAT  (»  •  ,3 4 X ,' A  (•  ,12  ,  • )  =  *,?  13.6) 

894  FOFKAT  ('- ' ,25X,'TKE  RICGE  POINT  GIVES  AN  IKFROVED  SCIUTICK* ,/,  *  ', 
3  25  X, ' 1R  Y  TO  ACCELERATE  ALONG  THE  FISC1UT1CN  FICGE') 

895  FCFHAT  ('-' ,25X,« ACCELERATION  ALONG  THE  RESOLUTION  FICGE  YIELDS  NO 

.3 1 3 PROVED  SOLUTION  •  ',25X,'THE  FIEC-E  FCINT  IS  THE  NES  EASE  FCIN 
.ST') 

896  FORM  AT  ('-•  ,25X, 'ACCELERATION  ALONG  THE  RESOLUTION  FIDGE  YIELDS  AN 
31 '1  FRO  VZD  SCI  U1ICN • ,/  ,  '  ',25X,'TRY  TC  ACCFLEFATE  rUFTFER') 

897  FCF8AT  ,25X,' DECISION  VARIA  ELES  ' ) 

898  FORMAT  (•-  '  ,25X, 'ACHIEVEMENT  VALUES') 

899  RETURN 
END 


C 

C  SUERCUTINE  C2JEUN  EVALUATES  EACH  C5JECTIVE 

€  FUNCTION.  THIS  SCEFCU1INE  IS  FFCEIIE  DEPENDENT 

C  AN!  BUST  EE  USER  SUPPLIED  FOR  EACH  SIGP  PROBLEM. 

C 

C 

CO M30N/C0M 202/X  (2500)  ,XEASE  (2500) 

CCEKCi;/CO320  6/C3J  (25CC)  ,  RES(25CC) 

C0K3CN/C0M21  4/A,3  ,C,E  ,E 
CEJ  (1)  =  X(  1)  >X  (3) 

03  J  (2) =X (2 )+  X  (4) 

CEJ  (3  ) -29  2  0.  *X(1)  +  177C.*X(2) 

03  J  (4)  =292 C. * X  (3)  *177C.*X(4) 

CBJ(5)=X (1) -0.064 52*  |X  (5 )  *X  (7  )  *X  (9  ) )  -0.  0£  250  *  (X  (  6)  *X  (6) ) 


C***' 


c 

c 

c 

c 

c 

c 


120 

L  =  Z[t)=x  (2) -3. 553 56*  (X(  1C)  +  X  (  12)  ♦  X  { It)  )  -  C  .  C  516  2*  (>  |1  1)  ♦  J  (13)  ) 
CaJ(7)  =X'3)-C.C6696*>  (15)  -0.06452*  (X  (16)  *X  (19)  )-C.  C625C*  (X  (17)  *X  (1 
33)  ) 

C3J  (3)  =  X(4)-0.  05362*X  (  20)  -C.C5556*  (X  (2  1)  *1  (24)  J-G.C5254*  (X  (22  )  *i  [2 
-?3)  ) 

A*  40.  *  (0.  9  99  73**  (X{5)  +  X  (  15))  *  C  .  9  9  5  53  <  «  ( 1  ( 1  C)  *  X  (2  0)  ) ) 

3  =  10.  *  (0.  9 9973**  (X  (6)  ♦  X  (16  ) )  *0.999  53**  (X  1 1 1)  *X  (2  1)  )  ) 

C=53.  *  (0.  99978**  (X  (7)  +X  (17) )  *C.9SS53**  |)  (12)  *X  (2  2) ) ) 

D  =  C. 

3=0. 

03  J  (9)  =  1C0.C  -  (A*E*C*C*S) 

5ETUEK 

2ND 

SC3E011IN2  UVAIU2 

■  ******************* *******«*********$***« 

SUEFCUTINE  UVALOE  RECALCULATES  EACH  CEoECUVE  HIHCC1 
ON  NECESSARY  2  VALDAHC  NS.  IBIS  S0EFCCT1KE  IS  FFCELEK 
DEI  END  ENT.  IP  10  =  1  IS  DEFINED,  CSEE  MCST  SCFFIY  1HIS 
SUBROUTINE  FOB  EACH  MG  I  FECEIEE. 


44*«4**4«»* 


c 

CC  EMO  K/CO  MXO  1/N03J,  NPRIOB,  HVAE.SHA  >,NC  YCIE 
C3MMON/CO.1M02/X  (2500)  , XBASE  (2500) 

CC  til  C  A/ COM  MO  6/C  SJ  (25CC)  ,£ES(25CC) 

commok/commi  c/opi  ,lme 

ccr.acN/coaai  i/life(25CC) 

COKrlCN/COMMl  2/AMAE  (1  CC,100) 

CC  t -1C  N/CG  1 1 1 3/S EP S 
C0 110  K/C0 121  4/A, 3 ,C ,D ,E 
L  =  LKUM 
DO  99  9  J=1  ,NC3 J 
DIFF(J)  =  0.0 

I?  (AMA?  (L,J)  .EC-0.0)  GC  TC  999 
GO  TO  (10  1,  10  1,  10  1,  1C  1,  1C  1,1  Cl,  1C  1,1  Cl,  1C9)  ,J 
101  GO  10  1 

109  GC  TO  (999,939,999,999,  11,  12,  13,  14,15,  1  1,12,  13,1  4,15,11  ,12,13,14,1 
35,11,12,13,14,15)  ,  1 

I  EIEF  (J)=A;iAP  (L,J)  *SEPS 
GO  TO  998 

II  DIFF  (J)=S  EFS*ALOG  (AM  AP  (L,  J)  ) 

A=  100.0-B-C-D-E-CPJ  (J) 

Y9  =  ALCG  (A)  *DIf  E(J) 

DIFF (J) =E XP  (19)-A 
GC  TC  998 

12  DIFP  (G)=SEFS*AIC£  (AMAI  (I  , J)  ) 

E=  TOO  .  0 -A-C- E-E-0  EJ  (J) 

Y9=ALCG  (3)  *DI?F  (J) 

DIFF  (J)=EXF  (Y9)  -E 
GO  TO  998 

13  DIFF  (J)=S  EFS*ALOG  (AMAP  (L,  J)) 

C=  100.C-A- E-D-E-OBJ (J) 

Y9=ALCG  (C)  *DIFF(J) 

DIFF (J) =E XP  (Y9J-C 
GC  TC  998 

14  GO  TO  999 

15  GC  TC  999 

998  OB  J(J)  =03  J  (J)OIFF  (J) 

999  CC  MINUS 


I 

I 

I 

I 

I 

J 


.  —  i  ‘J  £  4» 

END 

1> f  ***»?*♦***♦ 

SrJ  ERCUI  IN  5  A  C  fc  rU  N 
************** 


<<<«*********» 


SUBROUTINE  ACHFUK  fCr.KS  THE  ACHIEVEMET  VECTCF 
FEC2  THE  FSEVICUSIY  CALCULATED  VALCEE  CF  C5JFTJS 
AND  DEWAR.  S  CB5CDTI KE  CECrUL  ErTUElS  TEE  VALUES 
CE  THE  OBJECTIVE  EUHCTICNS  AT  TEE  CCFPENT  TEST  FT. 
SUBROUTINE  CEVVAH  FE1UFJS  THE  VALUES  CE  TEE 
DEVIATION  VARIABLES  AT  THE  CURRENT  TEST  POINT. 
THIS  S  U3R  C  t  U  N  F  IS  fSCEIZE  CEEEEEEFT  AN  t  MUST  EE 
USEE  SUtFLISD  ECR  EACE  NLGF  PRGEIE*. 


REAL*4  S 

CC  CMC  E/CC  M  MO  2/X  ( 2  £0 0)  ,  XE  AS  E  { 2  f  C C) 

COMaON/CO:iM0  3/A  (1  0)  ,  AEEST  J10) 

CC  EH  CL/ CO  H  20  7/N  (  2  £00)  ,?{  2SCC) 

CALL  DEWAR 

A  ( 1)  =  F(  1)  ♦  ?  (2)  ♦  ?  <  3)  +  P  (<l)  (£)  +  3  ( 6)  *N  (')  ♦  N  (6) 

A  (2)  =  3(9} 

F-TU  EN 
END 


] 
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